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SUMMARY 


Calculations are described of the remanent magnetization and the 
remanence after demagnetization as a function of field strength for 
powders of haematite, «Fe,0,. The calculations are based on the 
assumptions that the powder particles are single domains and_ that 
their hysteretic properties are due to magnetocrystalline anisotropy 
with threefold symmetry in the basal crystal plane. The calculations 
are intended for eventual comparison with experiments now in progress. 


§ 1. INTRODUCTION 


THE oxide mineral haematite, «Fe,O;, is one of the most important 
ferromagnetic constituents of rocks. Recently experiments have been 
started in the Department of Physics of this College, by Professor P. M. S. 
Blackett, F.R.S. and Mr. G. Haigh, designed to investigate more fully 
the magnetic properties of pure haematite powders. In particular, 
preliminary experiments have been carried out and more systematic 
ones are planned on the variation with field strength of the remanent 
magnetization and of the remanence after demagnetization, both at and 
above room temperature and also at temperatures below the Morin 
transition point (250°-260°K, Morin 1950). The present paper is 
intended to supply theoretical values for eventual comparison with 
the experimental results. The calculations are based on a number of 
simplifying assumptions which it may be necessary to modify in the 
light of these results. 

The following are the relevant properties of haematite, the assumptions 
used in the present calculations being briefly indicated and discussed : 

(1) The crystal structure of haematite is rhombohedral, with an axis 
of threefold symmetry perpendicular to a basal symmetry plane; the 
arrangement of iron and oxygen ions in the lattice is described by 
Nicholls (1955). 

(2) The magnetization « depends on temperature and field strength ; 
Néel and. Pauthenet (1952, cf. also Néel 1953) have shown that in 
moderately high fields the relation for o is 


o(T', H)=o0,(T)+x(T)H. Be ay Ae 
* Communicated by the Author, 
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The mass susceptibility y(Z7’) is anisotropic below the Morin transition 
point, 250°-260°xK, being approximately zero along the trigonal axis and 
isotropic and approximately constant, about 20x 107°, in the basal 
plane. Above this temperature y(7’) is constant and equal to about 
20 10-8, both along the trigonal axis and in the basal plane, up to 
high temperatures including the Curie point, about 950°K. These 
properties indicate that the paramagnetism of haematite is actually an 
antiferromagnetism (Néel 1953), the antiferromagnetic axis being parallel 
to the trigonal axis below the Morin transition and perpendicular to it 
above. This result is also borne out by neutron diffraction studies 
(Shull et al. 1951). 
The magnetization o,(7'’) may be written as the sum of two parts, 


Op=0;,+0,. Cet ade Se eee 


Here oc; is an isotropic component having a saturation value about 
0-2 e.m.u. per unit mass (1-0 per unit volume) and retaining this value 
both along the trigonal axis and perpendicular to it over a wide tempera- 
ture range, finally vanishing at the Curie point. On the other hand a, is 
zero below the Morin transition point but rises to a value about 0-3 e.m.u. 
per unit mass (1:6 per unit volumc) at this temperature, is anisotropic 
in the basal plane and finally vanishes at the Curie point. Along the 
trigonal axis no component of magnetization corresponding to co, can be 
observed in moderate fields at any temperature. The ferromagnetism of 
haematite is closely linked to its antiferromagnetism, but its origin is 
still uncertain. Néel (1949, 1953, 1955) ascribes it to one of two principal 
causes, ferromagnetic impurities such as Fe,0, or yFe,O3, or defects in 
the lattice planes containing the iron ions. 

The hysteretic properties of haematite above the Morin transition 
point and in moderate fields may be assumed to arise from the anisotropic 
component of magnetization co, in the basal plane. Since no information 
about the nature of the easy axes of magnetization nor of the magnitude 
of the corresponding anisotropy coefficients in this plane is available, 
it is here assumed on grounds of symmetry that there are three such 
easy axes making equal angles with each other. The simplest expression 
for the anisotropy energy commensurate with the symmetry is assumed 
(see relation (2.2) below) and the anisotropy coefficient K thus introduced 
is taken to be an arbitrary parameter the value of which should eventually 
be obtained by experiment (cf. § 2). ‘ 

Measurements of the magnetization in any field, e.g. the remanent 
magnetization, give the contributions due to both the anisotropic and 
the isotropic component. In order to obtain the contribution due to the 
anisotropic component alone a very convenient method of analysing the 
experimental results is available, suggested by Mr. G. Haigh: Since the 
isotropic component of magnetization o, varies little with temperature 
near the Morin transition point, while the anisotropic component oy 18 
zero below and approximately constant near and above this temperature, 
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measurements of the magnetization on either side of the transition 
temperature should yield at once the separate contributions o, and o,. 
_ The isotropic component is here assumed to make no contribution to 
the hysteretic properties, ie. it is assumed to be reversible; the 
preliminary experiments indicate that this assumption is not strictly 
true. ' 

Although the trigonal axis of haematite is a very dificult axis of 
magnetization, so that in moderate fields the magnetization vector of 
a crystal remains parallel to the basal plane, application of high fields 
inclined to the basal plane must cause some rotations of the magnetization 
vector towards the trigonal axis. Estimates of the anisotropy field 
corresponding to magnetization along the trigonal axis are provided by 
the results of ferromagnetic resonance measurements due to Anderson 
et al. (1954) who derive the value 3 x 104 oersted. The present experiments 
have been carried out in much more moderate fields so that it may be 
assumed as an approximation that the magnetization vector is always 
strictly parallel to the basal plane. 

To sum up, the hysteretic properties are here assumed to arise solely 
from the anisotropic component of magnetization in the basal plane, the 
isotropic component in that plane being taken to be infinitely ‘ soft ’ 
and the anisotropic component parallel to the trigonal axis infinitely 
‘hard’. 

(3) Although the mean size of the powder particles used in the 
experiments has not yet been determined, it may be assumed that 
they are single domains so that the magnetization change is one of 
rotation. No accurate estimate of the critical size of haematite powders 
has yet been made since the value of the anisotropy coefficient AK in the 
basal plane is not known, but the critical size, which is proportional to 
K1//[,, is probably large (a very rough estimate on the basis of the 
approximate K value deduced in § 2 is ly) owing to the small value 
of J,, the saturation magnetization of the anisotropic component, 
1-6 gauss. If the ferromagnetism of haematite is due to impurities 
or lattice defects it may be that each powder particle consists of smaller 
ferromagnetic islands embedded in a non-magnetic matrix (Néel 1955). 
These islands are then even more likely to be single domains. 


(4) Since the hysteretic properties of haematite are controlled by 
magnetocrystalline anisotropy, the saturation magnetization being very 
small, it is justifiable to neglect the effect of particle interaction on these 
properties (cf. Wohlfarth 1955 a). 


(5) In an assembly of powder particles it may be assumed that 
the orientations both of the trigonal axes and of the easy axes of 
magnetization in the basal planes are random. 

In §2 calculations are given of the curve relating the remanent 
magnetization I, with a reduced field strength and in $3 of the curve 
relating the remanence after demagnetization / with field strength. 
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§ 2. REMANENT MAGNETIZATION 


In the experiments a field of strength H is applied to the initially 
demagnetized powder and the remanent magnetization Ip measured 
after removal of the field. The curve obtained below gives the contri- 
bution to the remanence of the anisotropic component of magnetization 
in the basal plane, obtainable from the measured values as discussed 
ial Ale 
Figure | shows the angle notation, ~ being the angle between H and 
the trigonal axis, 6(0<@<7/3) that between the positive direction of 
the component of H in the basal plane and the nearest easy axis A 
and ys, that between the magnetization vector 1, and A. 


Fig. 1 


trigonal axis 


easy xis 
A 


Angle notation. 


The total energy is the sum of the energy in the external field Ey 
and the anisotropy energy #4, where 
Hy=—HI, sin « cos (0+y,). youre err 4 6) 


Also the lowest non-vanishing term in the general expression for the 
anisotropy energy, which alone is considered here, is of the form 


E,4=K sin? (3,/2), ee eo ee 

K being the anisotropy constant. For equilibrium 
Asin x sin (0+-,)+4 sin 34,=0, eee ont Me 
A sin x cos (0-+-4,)+3 cos 3. >0, Ne OR 


where h=HI 3K. Oe ke com eee 
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The rotations of the magnetization vector with changing field strength 
become discontinuous when the inequality in (2.4) becomes an equality, 
the critical value of 4. denoted by i,(6). being given by the relation 


3 tan (0+ %)=tan 35. Suet ate et a) 


The corresponding reduced critical field strength, denoted by f(A), is 
obtained from (2.6) and (2.3) or (2.4) with an equality sign. If 


x= —h,(6) sin «>0, 
then 
6=tan~* (4t)—4 tan-! 7, ] 
where ; (2.7) 
j=tan Bfy= {(9—42°)|(4a2—1)}”, | 
With 
O<0en/3, OSyoen/6 and 3321. 


In the completely demagnetized state the magnetization vectors for 
one-third the particles with orientation angle @ are along the easy axis A, 
for one-third along B and for one-third along C, the angle between any 
two easy axes being 27/3; see fig. 2 showing the basal plane. 


a aa 
field 
component 


basal plane AS 


Angle notation for the basal plane. In the demagnetized state the magnetization 
vectors J, of one-third the particles are parallel to each of the three 
easy axes A, Band C. On application of a field in the negative direction 
the magnetization vectors rotate towards the field direction, making 
angles ,, pp and 4%, with A, B and C respectively. In the state of 
remanence after saturation in the positive direction the magnetization 
vectors are parallel to A. 0<0<7/3. 


When the field component H sin « is applied, for convenience, in the 
negative direction in the basal plane (see fig. 2) the magnetization vectors 
rotate towards this direction until, when the field has reached a critical 
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value, ie. when h=HI,/3K=|h (0) |, those originally parallel to A, 
now making an angle 4,—y9() with A, and those originally parallel 
to B, now making an angle 4,—7/3—z%(9) with B, rotate discontinuously 
towards C. After the discontinuous rotations have taken place these 
two sets of magnetization vectors are parallel to each other and to the 
magnetization vectors originally parallel to C. On now reducing the 
field to zero the total magnetization vector rotates continuously back 
to GC. For particles with orientation angles 6 and « the component of 
the remanent magnetization along the direction H, /;(6,%), is thus 
given by 

ERO, 4) /Lo=Jn (9, %)=9, h< | ho(?) z 


(2.8) 
jg(O, %)—=cos (7/3—8O) sina, AS | h,(A)| 

h,(@) being given by (2.7), where x=|h,(0)|sina. If the orientations 

of the particles are random then the reduced remanent magnetization 

of the assembly jp=Tp/[y is given by 


8 pal2 ra/3 - 
Cai in? Cre 
IR a sin a. da.| cos (; a) dé, 


/6 


3 pal é a Mi: 
=-). sin 3 9 sin? « da, ee 
3 rh 7 
eat A CR see Ni AG 
ah 2 sin( ath 2) ‘hae, 
where h = |h,(0) | =H1)/3K, x=hsina, and where @ is given as a 


function of x by (2.7). 

The analogous much simpler calculation for particles with two easy 
directions of magnetization (Stoner and Wohlfarth 1948) has been given 
by Wohlfarth (1955 b). In the present case the main mathematical 
difficulty lies in the evaluation of the integral in (2.9). 

Since $<a< $ it follows from (2.9) that 


Jn=9, BS 

3 
Jr==1(h), 

hg’ fa 34/3 (7 1( a 
jn= 1 5 5 —sin- <n h>8, 
where 

Stee | 
I(€)=h | 1o5 sin 3 —9 [nr—arynh da, 

2.10 


; As Pes a eee ; 
9 being given by (2.7). By using as variable of integration the Ree. 


vin 
nee 5 —3%)=tan-! {(4a?—1)/(9—4a? yes 23 2 (Sh1g) 
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it is found that 


Jr=9, h<}, | 
: 4,/2 k | 
ju= —— pail), $<h<3, 
p42 34/3 
jr= (5)+ IG —sin-¥(3/2h)—3(4l2—9)19)402}, h>8, 
where | 
4 
Jin)=| {1—k~ sin? ¢}/2 cos ¢ cos (44) dd, | 
0 
and where 
k= {(4h?—1)/8}1/2=sin Bp. 
(2.12) 


The integral J,(y) cannot be evaluated in closed form except when 
h=3/2; it is, however, possible to express cos (4f) over the range 
required in a more amenable approximate form, viz. 


where (2.13) 


cos (1¢)=A-+B cos 6+C sin? ¢, | 
A=0-8400, B=0-1600, C=0-0267, | 


which gives the required values of cos (44) correctly to at least three 
significant figures and the value of the integral to an even higher accuracy 
(for h=3/2, the exact value of J, is 0-7714, the approximate value using 
(2.13) is 0-7716). Inserting the approximation (2.13) into (2.12), jp can 
be evaluated in terms of elementary and elliptic functions. For (h—4)+0 
and h->co the following series expressions for jp are obtained : 


2k? 
jr= i {1—0-0133k2+ 0(k*)}, k= (4h2—1)/8, 
(2.14) 
= 24/2(h—$){1—3-007(h—3)+O(h—$)?}, (h—-3) > 0, 
jp=3-/3/8—0-0933h-3—0-0225h-5—O(h-’), h—> ©. 


The calculated values of j, are given in the table and graphically in 
fig. 3. 

The most striking characteristic of the jp, curve is the low field 
foot, extending up to h=}. Thereafter jp rises, linearly at first with 
slope 2\/2—2-828, and finally tends asymptotically to the saturation 
value 34/3/8=0-6495. 

The preliminary experiments indicate that the field strength for which 
jn tises appreciably has a value in the range 1 to 2x 10° oersted. From 
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this the value of K is equal to roughly 10? erg em~1!, compared with 
5x 10° erg em=! for the room temperature value of iron. 


Fig. 3 


h IR Jp h JR JD 
0-50 0-000 +0-650 1-1 0-557 — 0-424 
0-52 0-053 0-564 1-2 0-579 0-477 
0-54 0-101 0-484 1-3 0-596 0-518 
0-56 0-143 0-411 1-4 0-608 0-550 
0-58 0-182 0-342 1:5 0-618 0-572 
0-6 0-216 0-280 2-0 0-637 0-620 
0:7 0-346 +-0-032 2°5 0-643 0-635 
0-8 0-430 —0-140 3-0 0-646 0-641 
0-9 0-487 0-264 o'e) 0-650 —(-650 
1-0 0-527 —(0-355 


Dependence on reduced field strength h=HJ,/3K of the reduced remanent 
magnetization jr=Ip/I, and the reduced remanence after demag- 
netization jp=Ip/Ig. For h<$, jp=0, jp=0-650. 


§ 3. REMANENCE AFTER DEMAGNETIZATION 


In the second type of experiment the powder is reduced to the 
remanent state from complete saturation. A reverse field H is then 


applied and the remanent magnetization 7, measured after removal of 
the field. 
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In fig. 2 showing the basal plane the positive direction is that of initial 
saturation. On removing the field required for saturation the magnetiza- 
tion vectors of particles with orientation angle @ turn towards A, the 
nearest easy axis of magnetization. After now applying and then 
removing a field with component H sin x in the negative direction the 
magnetization vectors will be parallel to the easy axis C, that nearest 
to the reverse field direction, if H is larger than the critical value, 
depending on @ and « and given by (2.7). Hence the reduced remanent 
magnetization jp7=Ip/I) is given by 


See ae 
j= ~| {sin—sin G -#)| sin? « da 
TJ9 3 ( i) aie a Cosh) 


=jp —Ir: } 
Here jp has been evaluated in § 2 and 


3 h | ] 
jp Sa ae | a sin 6/(h?—2?)'/2} da, 
JD a ma /( ee ? | eee) 
where h= |A,(0) | =H, /3K, x=h sin a 
and where @ is given as a function of x by (2.7). By using the same 
transformation as that leading to (2.12) it is found that 


jp =3v/ 3/8, h<d, ] 
34/3 44/2.k | 
= ive {sin~1 (1/2h)+ (4h?—1)1/?/4h?}— ak, ja? 2/P): Loh<, : 
mee on 4/2 kb _ [r\ | 
ane we {sin-! (1/2h) + (4h2—1)/2/4h2}— Sun i (5). h>8, 
where 
4 
Jn)={ {1—k- sin? ¢}1/2 cos ¢ cos (2 as, 
0 
and where 
k= {(4h?—1)/8}1/2=sin B. j 
(3.3) 
By using the approximate expression 
COs ( 3 S\=A+B cos 6+C cos? +D sin +E sin ¢ cos 4, 
(3.4) 


where A=0-2874, B=0-3216, C=—0-1090, D=—0-2874, 
EK=0-0036, 


which gives sufficient accuracy over the range required, the integral J, 
may be evaluated and hence jp’ and, from (3.1), jp calculated. The 
values of jy are given in the table and graphically in fig. 4. 

It is seen that the reduced remanence after demagnetization retains 
its saturation value 34/3/8—0-6495 until, when h exceeds 3, 7p decreases, 
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linearly at first with negative slope 31/2=4-243, passes through zero 
when h—0-720 and finally tends asymptotically to its saturation value 
__0-6495 in the reverse direction. The preliminary experiments indicate 
that the field at which j)=0 is rather larger than that calculated and 


shown in fig. 4. 


Fig. 4 


0.0 0-5 1.0 h 5 2.0 2.5 


Reduced remanence after demagnetization. 
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ABSTRACT 


The new radioactive isotope 18*°Ta has been prepared and its decay 
characteristics determined as follows :— 


Half-life : 10-5-L0-5 min 

Maximum f-energy : 2:2. Mev 

Conversion-electron energies : <0-15 Mev 

y-ray energies : 125, 200, 300, 410, 510, 610, 730, 940 and 


possibly ~1150 kev 
The mass assignment suggested by its production by (n, p), but not 
by (y, p), reactions on tungsten was confirmed by experiments using 
tungstic acid enriched in 186W, 


§1. INTRODUCTION 
Or the neutron excess tantalum isotopes which can be made by fast 
neutron irradiation of tungsten, only 18®T'a has not previously been re- 
ported. This paper describes the production of this isotope and the 
examination of its decay characteristics. 


§2. EXPERIMENTAL TECHNIQUES 


The neutron excess tantalum isotopes were produced by irradiation of 
tungstic acid with fast neutrons generated by proton bombardment of a 
beryllium target in the Harwell cyclotron. The tantalum isotopes were 
quickly separated by coprecipitation with ferric hydroxide from the 
solution of tungstic acid in sodium hydroxide. Suitable samples for 
counting were obtainable within about seven minutes of the end of the 
irradiation. 10-minute irradiations of 5-20g of tungstic acid with 
neutrons of energies up to 50 Mev produced sources, with intensities of 
about 1 yc, for the examination of the decay characteristics as described 
below. 

§3. RESULTS 
(i) Production and Mass Assignment 


After irradiation of tungstic acid with neutrons of energies up to 
20 Mev, so that tantalum isotopes were produced mainly by (n, p) reactions 
on the tungsten, a half-life of about 10 minutes was observed in the 


- * Now at the Chemistry Dept., Imperial College, London, S.W.7. 
+ Communicated by the Author. 
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beta-decay of the tantalum sources in addition to a 9-hour half-life and 
longer-lived activities. This 10-minute half-life was also observed when 
low energy conversion electrons from any 16-minute '*°™Ta present were 
absorbed by 60 mg/cm? of aluminium. The decay of the more intense 
sources produced by irradiations with neutrons of energies up to 50 Mev 
also showed the 10-minute half-life when the particles counted had passed 
through 500 mg/cm? of aluminium, which was more than enough to absorb 
any 0-6 Mev beta-particles from 1°™Ta (Wilkinson 1950). Spectrographic 
analysis of the tungstic acid, and test irradiations of iron, showed that the 
10-minute activity could not have been due to contamination by 9-minute 
Ee: 

The 10-minute activity was not observed after irradiation of tungstic 
acid with an intense beam of x-rays, of energies up to 28 Mev, generated 
by a linear electron accelerator. Its production by (n, p), but not by 
(y, p), reactions on tungsten suggests that the 10-minute activity should be 
assigned to !*4Ta or 186Ta but 1**Ta is known to have an 8-7-hour half-life 
(Butement and Poé 1955). That the 10-minute activity must be assigned 
to 18°Ta was confirmed by experiments with tungstic acid enriched in 
186W, In one experiment equal masses of normal and enriched tungstic 
acid were irradiated by neutrons with energies of up to 20 Mev. The 
irradiations were monitored with the 12-8-hour activity induced in copper 
foils placed between the neutron source and the tungstic acid target. 
Tantalum activities were separated from the two samples under closely 
similar conditions and the decays of the beta-particles passing through 
30 mg/cm? of aluminium were followed. The ratio of the yield of the 
10-minute activity from the enriched tungstic acid to that from the normal 
tungstic acid was 6-4: 1. Since the enrichment factors of the tungsten 
isotopes of mass 180, 182, 183, 184 and 186 were 0-09, 0-02, 0-03, 0-05 and 
3-4 respectively, relative to natural tungsten, the increase in yield of the 
10-minute activity was closest to the enrichment of 18¢W, 

In the other experiment the ratio of the intensity of the 10-minute 
activity to that of the 50-minute '%°Ta activity produced by (n, pn) 
reactions on '*®°W was measured for tantalum samples prepared by 
irradiation of normal and '86W-enriched tungstic acid with neutrons of 
energies up to 50 Mey. Low energy electrons from the tantalum samples 
were absorbed by 60 mg/cm? of aluminium. In both cases the ratio 
was the same indicating that the 10-minute tantalum activity was also 
produced by fast neutron bombardment of 18*W. 

The 10-minute activity can therefore be unambiguously assigned to 
186'T'g,. 

(ii) Radiation Characteristics 

The half-life was obtained by measurement of the decay of beta- 
particles passing through 500 mg/cm? of aluminium which absorbed all 
particles from isotopes other than the 10-minute 18°Ta and the 50-minute 
5Ta. The half-life found as the average of eight determinations was 
10-5 minutes, all the results lying within 0-5 minutes of this value, 
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The aluminium absorption curve of the beta-particles from short- 
lived tantalum activities indicated a maximum beta-particle energy of 
2-3 Mev together with a low energy component of about 0:15 Mev with 
approximately the same intensity. Beta-gamma coincidence absorption 
experiments, using a Geiger counter and Nal crystal scintillation counter, 
showed that the high energy beta-disintegration is in coincidence with 
gamma-rays. While some of the low energy particles must have been 
conversion electrons from !**™Ta many were conversion electrons in coin- 
cidence with the 10-minute !*°Ta beta-disintegration. This was shown by 
the aluminium absorption of those particles entering one end-window 
Geiger counter in coincidence with beta-particles entering a similar counter 
from which conversion electrons and low energy beta-particles were 
excluded by aluminium absorbers. Spurious coincidences caused by 
scattering of a particle from one counter into the other were almost 
completely eliminated by mounting the source over a small hole in a brass 
plate placed between the two counters. More than 90°% of the coinci- 
dences were due to particles with energies up to about 0-15 Mev. There- 
fore most of the conversion electrons from the various conversion 
processes possible, in view of the gamma-ray spectrum described below, 
have energies of about 0-15 Mev or less. 


1000 
200 kev 


410key 20 key 
610kev  730kev 


100 


So 


~II50 kev 


Relative counting rate / unit energy range 


500 1000 
Energy kev 


Gamma-ray energy spectrum of 1°Ta 


The beta-particle energy spectrum was examined from 1-25 Mev 
upwards with an anthracene crystal scintillation spectrometer together 
with a 30-channel pulse analyser. The spectrometer was calibrated with 
conversion electrons following the decay of !!Sn and 87Cs and beta- 
particles from Ni (2-10 Mev) and 90Y (2:24 mev). Fermi beta-disinte- 
gration functions given by Feister (1950) were used for the Fermi plots. 
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The maximum beta-particle energy of 18°Ta was found to be 2-2. Mev as 
the average of twelve closely agreeing results. 

The gamma-ray spectrum of !8°Ta was examined with a Nal crystal 
scintillation spectrometer and a 30-channel pulse analyser. The spectrum 
was covered in several overlapping stages and at two different amplifier 
gain settings, two or three successive counts of the short-lived activity 
being taken from each sample and corrections for longer-lived activities 
being applied. The spectrum up to 300 kev was obtained from samples 
prepared by irradiation of !8*W-enriched tungstic acid to avoid interference 
from the 180 kev gamma-ray from !°°™Ta. The results were all combined 
to give the spectrum shown in the figure. No 60 kev K x-ray component 
was observed because of the very large background intensity of this energy 
from 50-minute !8°Ta. The spectrum is complex and may be interpreted 
to indicate gamma-rays, or groups of gamma-rays, of energies close to 
125, 200, 300, 410, 510, 610, 730 and 940 kev. After allowing for Compton 
absorption and for photo-electric absorption efficiency the relative 
intensities of the above gamma-rays were found to be very approximately 
2:5:10:2-5:2:4:5:45:6-5:1-5. The 125 kev component is probably 
at least partly due to the transition of that energy which has also been 
observed following electron capture by '*Re (Metzger and Hill 1951, 
Steffen 1951) and Coulomb excitation of 18®W (McClelland, Mark and 
Goodman 1954). 

The higher energy component at about 1150 kev was not clearly 
resolved and not very reproducible. It was probably due in part to 
simultaneous absorption in the crystal of two lower energy gamma-rays. 
Experiments with the source | in. from the crystal suggest, however, that 
there is a genuine gamma-ray, or group of gamma-rays, at 1100-1200 kev 
with intensity <1 relative to the 730 kev gamma-ray intensity of 6-5. 
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SUMMARY 

The linear strain hardening observed after the easy glide region in face- 
centred cubic single crystals at sufficiently low temperature is explained 
by assuming that all the Frank—Read sources of the active slip system are 
im action, and form small active slip planes bounded on all four sides by 
Cottrell barriers. The initial easy glide region should occur before all 
the barriers are formed ; the parabolic strain hardening appearing in the 
final stages would be related to a thermally activated breaking of the piled- 
up dislocations through the Cottrell barriers. 


$1. INTRODUCTION 
THE stress-strain curves observed by traction of face-centred cubic single 
crystals exhibit, after the initial easy glide stage AB (fig. 1), a linear part 
BC which has as yet received no satisfactory explanation. 


Fig. 1 


Typical stress-strain curve of a face-centred cubic single crystal 


It is usually agreed that easy glide occurs as long as the dislocations of 
the principal slip system can develop without much hindrance, and, 
at least in face-centred cubic crystals. either go out of the erystal or 
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compensate each other’s stress field in deformation bands. Also greater 
hardening rates occur when some ‘secondary’ slip on another system 
hinders the principal slip system too much. In agreement with this idea, 
easy glide is in evidence particularly under conditions such that a principal 
slip system is much more favoured than the others, the extreme case being 
that of hexagonal metals (Mott 1953, Cottrell 1953). The second stage BC 
actually occurs before the axis of traction has reached the [100]-{111] 
symmetry zone, where two slip systems are equally favoured. But this 
is usually explained by inhomogeneities in the stresses due to the disloca- 
tions developed in the principal slip system ; or, more macroscopically, by 
a lower work hardening in the secondary slip systems (Rdéhm and 
Kochendorfer 1948, Stroh 1953). 

It is now clear that the more rapid work hardening rate, which sets in 
after the easy glide stage, begins with an almost perfectly linear stage 
BC (fig. 1) in all the face-centred cubic crystals investigated and whatever 
the initial orientation of their axis of traction, if the temperature of 
deformation is sufficiently low. The metals and alloys studied so far are 
aluminium, nickel, copper, silver, gold, alpha brass and Cu, Au (c.f. Masima 
and Sachs 1928, von Goler and Sachs 1929, Sachs and Weerts 1930, 1931, 
Blewitt and Kohler 1950, Jaoul 1955, Kuhlmann-Wilsdorf and Wilsdorf 
1953, Blewitt 1953, 1955, Rosi 1954, Staubwasser 1954, Kocks 1954, 
Cottrell and Stokes 1955, Diehl, Mader and Seeger 1955). The rate of 
work hardening (do/de) is about 4 10-21, increasing slightly for axes of 
traction with orientations varying from [110] to the [100]-[111] zone ; 
but the factor of proportionality to the shear modulus p does not seem to 
vary appreciably with temperature, nor with the nature of the metal. On 
the other hand, the last stage CD (fig. 1), with negative curvature, appears 
at stresses lower for higher temperatures and, at the same temperature, 
lower for aluminium than for copper. These two effects combine so as 
to reduce the linear stage BC to practically nothing in aluminium at room 
temperature, and to produce twinning without any third stage CD in 
copper at 4°K. 

Some current theories of work hardening (Taylor 1934, Mott 1953) 
predict parabolic stress-strain curves, and cannot therefore explain the 
linear stage BC. This note suggests a possible mechanism for this 
behaviour, and for the variation of the parabolic stage CD with tempera- 
ture and with the nature of the metal. The main difference with the 
previous theories will be to assume, in the linear stage, that piled up 
groups of dislocations are fixed in number and position, and increase only 
in strength with deformation. Previous theories assumed on the contrary 
that the dislocations were dispersed (Taylor) or piled up in groups of fixed 
strength (Mott) but increasing in number with deformation. 


§2. COTTRELL BARRIERS 


Work hardening is certainly due to an accumulation of dislocations in 
the crystal. This can impede further movements of dislocations in two 
ways (c.f. Mott 1953, Cottrell 1953, Seeger 1954) : 
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1. Multiplication of the forest of screw dislocations through which the 
moving dislocations have to cut. 

2. Increase in the elastic distortion of the crystal, which hinders the 
movement of dislocations. 

The creation of jogs necessary in the first process is a thermally activated 
process, with quite a small activation energy for metals like aluminium ; 
it should lead to a work hardening strongly dependent on temperature, 
contrary to what is observed in the linear part BC. A dependence on 
temperature proportional to that of the elastic constants, as is predicted 
by the second mechanism if no recovery occurs, seems in better agreement 
with observation. 

For this second mechanism to be satisfactory, one must find a reason, 
other than cutting through the forest of screws, for dislocations to remain 
in the crystal. In single crystals, this can be achieved only by elastic 
interaction with dislocations of the same slip system and opposite sign, 
as in deformation bands, or by piling up dislocations in front of Cottrell 
barriers formed with dislocations of other slip systems. Deformation 
bands cannot harden the crystal very much, for the long range actions of 
their dislocations with opposite signs compensate each other (Mott 
1953). They form indeed in the easy glide stage. The beginning of the 
linear stage BC must therefore be connected with the appearance of Cottrell 
barriers. 


CS (101) 
C4 va 
(Ot!) 


Slip band in a (111) plane. 


A great number of such barriers must be formed from the beginning of 
the linear stage. For, as will now be shown, every Frank—Read source 
belonging to the principal slip system must be active; and every active 
source must form at least four barriers in its glide plane. 

The first conclusion is reached by noting that if, in the linear stage, 
dislocations cut through the Frank—Read net without difficulty even at 
low temperatures, they must be acted upon by stresses large enough to 
make active most of the Frank—Read sources of the principal slip system. 
This is not surprising, as the stresses applied in the linear stage are always 
much larger than 10~ wu. 

An active source 8 can, on the other hand, form barriers along the two 
[110] directions which are not parallel to its own Burgers vector (c.f. 
Appendix A). It can therefore surround itself with four barriers, forming 
the parallelogram represented in fig. 2. 


4K2 
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Now, if hardening is only due to piled-up groups in front of Cottrell 
barriers, each active source must be surrounded by such a parallelogram. 
Otherwise, loops from S could escape to the surface of the crystal, and 
produce an infinite amount of deformation under a constant stress. 


§ 3. SHoRT SLIP 

The size L of these parallelograms may be deduced from the two 
following considerations : 

1. As shown in Appendix A, the density Ng, of sources 8’ able to form 
barriers of a given direction with the sources S is roughly the same as the 
density Ng of S sources. To build all the parallelograms required, all the 
possible S’ sources must therefore be used, and each of them must build 
2N,/Ne& 2 barriers with active sources 8. That loops from a given source 
can build barriers with loops from several other sources is shown clearly 
by Jacquet’s observations on alpha brass (1954). 

2. Cottrell barriers are formed readily, as soon as two intersecting slip 
systems of suitable slip directions are active. This is again shown by 
Jacquet’s pictures. This must be due to the attraction between disloca- 
tions able to form barriers readily (c.f. Appendix A) and to the con- 
siderable amount of energy gained by forming such a barrier (usually 
several electron volts per interatomic distance). We can therefore assume 
that, from the start of the linear stage, each active source S has built up 
barriers with the nearest 8’ sources available. In these conditions, the size 
LL of the parallelograms is the average distance between S’ sources : 


L= (Ng). 


The linear stage should therefore set in when the orientation of the 
axis of traction and the hardening due to easy glide are such as to activate 
the sources 8’ of the required slip system to form barriers. It should be 
a region of short slip bands, of constant number and length, distributed at 
random through the crystal. 

For a well annealed crystal, we may assume that the possible Frank— 
Read sources build up a three-dimensional net of size 1, where the various 
slip systems are equally represented (Mott 1953). Then Ny=WNg and, 
when one of the 12 slip systems is active, Ng~1/12/%. With /=10-4 
to 10-3 em, the lengths and distances of the slip bands should be of about 
10-3 cm. 

These conclusions, deduced from an analysis of the stress-strain curves, 
are in full agreement with observations made independently on copper by 
Blewitt (1955), Blewitt examined the slip lines produced by an elonga- 
tion of a few per cent, after prior traction into the linear stage followed by 
electropolishing. The surface was thus free in these experiments from 
the slip lines produced in the easy glide stage. It exhibited all the 
characteristics predicted, except apparently for some reduction in size of 
the slip lines with increasing strain. As pointed out by Blewitt himself, 
however, the dispersion of the statistical measurements made this last 
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effect uncertain. As expected too, the slip lines in the easy glide stage 
were quite different, strong lines, much longer and fewer in number, 
appearing in what is, according to Diehl, Mader and Seeger (1955), a 
background of fine slip. 

The hardening of polygonized crystals will be considered in § 5. 


§ 4. Linear HARDENING 


In our model, hardening may be attributed to the internal stresses 
acting on the sources 8. They have three possible origins: the back 
stress o, of a piled-up group on its emitting source; a ‘relaxation’ 
a,’ of this stress, described below ; and the stress c, on the source due 
to all the other piled-up groups. — If, as assumed here, all the piled-up 
groups have similar numbers n of loops, » increasing uniformly with 
strain, these three stresses are proportional to n. This explains a linear 
hardening, for the strain varies linearly with n. 


Fig. 3 
{] 
S 


S 


Relaxation around a piled up group. 


Treating the n loops of a piled up group (fig. 2) as concentric circles, we 
have (Leibfried 1951) 


b 
oe eee 5 o/s C1) 
As the strain due to the piled up groups is 
e=nN,L*b, eT re irs ga eee) 
the first term, taken alone, would give a hardening rate 
se Sra 1 hs ok Se TE) 


ema ONGLAs 
For a well annealed crystal, where Nga is of the order of Ns, 
a ~41, which is too large by a factor about 10°. 
€ 


We attribute the main part of this discrepancy to a relaxation of 
c,, produced by the neighbouring Frank—Read sources. As pointed out 
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by Mott (1953), as soon as a piled up group has more than a few disloca- 
tions, its stress field is strong enough to set into action all the neighbouring 
sources 8” (fig. 3). These emit loops of dislocations in such a way as to 
reduce the stresses around the piled up groups. 

Practically all the Frank—Read sources can be used for this relaxation. 
In a well annealed crystal and if each source 8” contributes loops to the 
two nearest piled up groups, there will be 12-l1=11 sources 8” per piled 
up group, if one of the 12 slip systems is active. The average distance 
D over which relaxation of a piled up group occurs should therefore be of 


1 2 OV 
the order of a fraction 57] SL of the distance between S” sources, so that 


D=(I/20), if lis the size of the Frank—Read net. We may reasonably 
assume that the relaxation transforms the effect of the piled up group on 
its S source into that of a dipole of average strength nD. The back stress 
on § is thus reduced from o, to o,-+-0,'=(D/L)o,; and, as Lo1 x (12), 


Ds 1 I 
L ~ 30% (12)13 — 50° 


The contributions to the strain « from the various S”’ sources during 
relaxation probably cancel each other. Thus « is still given by (2). 
And the relaxation reduces the hardening rate by the factor D/L=2 10-?, 
of the right order of magnitude to explain the rates observed : 

do FUORI ee ty fh ee oe 
de de 2L 

The stress o, due to the other (relaxed) piled up groups has been 
neglected in the hardening rate, for the following reason. These groups 
occur in pairs of opposite sign, and therefore cancel their effects, except 
for the group nearest to the S source considered, at a distance, say. A 
(fig. 6). Thus, for a well-annealed crystal, 


nubD : 
to OPT 4 a! Sica! bapele gee otty ars (5) 


A is of the order of the distance between § sources, thus L, according 
to §3; and ay is of the order of o,+¢,', but with an equal probability 
of being positive or negative. If, therefore, o is increased by do, the 
sources for which the corresponding increase do, is positive or negative 
will be present in equal number ; they will emit respectively dn + or dn_ 
loops, so that 


bl pbD 

do=do,+ > dn = —doy+ dn 
— hOD dn,+dn_ wD a 
Ter Oy 


The stresses o, due to the other piled-up groups should not therefore 
effect the hardening rate, which is the same as (4). They may, on the 
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other hand, play an important role in stabilizing the piled up groups : 
one has to explain why the piled up groups remain in the crystal when 
the applied stress is removed. Mott (1953) proposed that piled up groups 
could be locked by Cottrell barriers formed on their rear with loops 
emitted by 8S” sources during relaxation. This process is however not 
possible with the model used here, for all the 8’ type of sources able to 
form barriers with the dislocations of the piled-up groups have been 
used up to form the parallelograms on which piling up occurs (§ 2). The 
loops emitted by the 8S” sources during relaxation interact therefore 
only elastically with the piled up groups*. They reduce their energy, 
thus making them more easily stabilized by any frictional force such as 
the interactions o, with other piled up groups. This possible stabilization 
is, however, not yet clear in detail. 


$5. HARDENING OF POLYGONIZED CRYSTALS 


The equations of the last section can be put into a more general form, 
which shows that the hardening rate should not depend very much on 
the distribution and size of the Frank—Read sources in the metal. 

Let Nz, Ng and Ng., be the densities of 8S, S’ and 8’ sources in the 
erystal (S active sources; S’ sources able to form barriers readily with 
S sources ; 8” sources other than $8). VD is, according to § 3, a fraction 
N,/2N,-. of the average distance Ng,,-'/? between S” sources. Thus 


D=Ng/2(Ny)*? 


da be 
—- = Wan? 
de ALALN 3”) 


with L=(N,.)-V3 if Ng. <Ng and N,/N,4" if Ng, Ny. The value of L 
for Ns, >N, comes from the fact that there will be Ng,/Ng sources 
S’ able to form a barrier with each source 8, and their average distance 
is N,V. The value of L for Ng.<Ng was obtained in § 3. Finally, 
if the sources 8, and 8, able to form barriers along the two directions 
[110] and [011] of fig. 2 have unequal densities, thus give unequal slip 
distances L,<JL,, the factor L4 in the hardening rate is to be replaced 
by L,3L,. For the variation of o only depends on the action of the 
piled up groups at the shorter distance L,; and L, enters into the 
equation only through the area LL, of the loops. 

These equations were obtained assuming the S, 8’ and 8” sources 
distributed fairly wniformly in the crystal. With this assumption, some 
interesting results are obtained. 

1. The hardening rate do/de depends only on the densities of S, S’ and 
8” sources, not on their exact geometrical arrangement. Thus it is the 
same whether the crystal is polygonized or well annealed, provided Ng, 
Ny, Ng are the same. 


and 


nn nnn nnn nnn EEE 


* They may form Cottrell barriers between themselves. 


1176 J. Friedel on the Linear Work 


2. The hardening rate depends only on the ratios of these densities, 
not on their absolute magnitude : it does not depend on the size of the 
Frank—Read net. 

3. In the general case where Ng-<Ng, the hardening rate does not 
depend on the density of the active S sources. Thus a (finely) polygonized 
structure where the walls of the blocks are built up with dislocations of 
the active slip system should have the same hardening rate as a well 
annealed crystal with the same densities of dislocations in the other slip 
systems, although its density N, of S sources is much larger. This is 
because the lowering of the hardening rate do,/de due to the piled up groups 
(cf. (3)) is compensated by a decrease in relaxation (o,+0,')/o,;—D/L, due 
to fewer 8’’ sources per piled up group. There are, however, many more 
slip lines than in the well annealed crystal, for the number of slip lines is 
given by Ng. 

This is possibly Kuhlmann-Wilsdorf and Wilsdorf’s ‘fine slip’, 
observed on mono- and polycrystals of high purity aluminium after 
fairly large deformations (cf. Kuhlmann-Wilsdorf, van der Merwe and 
Wilsdorf 1952). It consists of parallel slip lines, of length L from 5 to 
20 10-4 cm and with distance X from 3 to 5x 10-® cm apart. It could 
be originated by a polygonized structure where the blocks of the mosaic 
structure, of size J~Z/8,/12, are covered by an array of 1//X~100 
parallel dislocations belonging to the active slip system. High purity 
metals are known to present readily a polygonized structure of that 
type, when heat treated for recrystallization (Talbot, de Beaulieu and 
Chaudron 1953, Friedel, Boulanger and Crussard 1955). Finally the 
number » of loops deduced from the height of the slip lines in the fine 
slip is near to Lo/ub, in approximate agreement with eqn. (1). There is 
therefore, as expected, little relaxation of the stress o, of the piled up 
groups in that case : 


o/0;=(0,+0,')/o,=D/L = 1. 


4, The only case where V,,> Ng is in a polygonized structure where 
the dislocations of the walls belong to one of the two types of S’ sources. 
The equations given above predict then a very short slip distance L, in 
that direction, and as trong hardening rate, varying as (Ng.)19/8, It 
would be interesting to check this point. However when the active slip 
system is sufficiently: hardened in that way, other systems become 
probably very active and distort the picture. 

5. All other cases of polygonization should lead to hardening rates of 
the same order as for well annealed crystals, if somewhat smaller. This 
is because the density of S’’ sources is large anyway: to lead to a 
noticeable decrease in the hardening rate, the polygonization must 
increase the density of sources, in one of the 11 slip systems making up 
the S’” sources, by a factor much larger than 10. This is not a very 
likely occurrence. 
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§ 6. BREAKING UP oF CorrrELt BARRIERS 


The ‘ parabolic’ hardening CD (fig. 1) occurs at a stress og which 
varies strongly with temperature, thus indicating some thermally 
activated process. We suggest that it occurs when the stress concentration 
om front of the piled wp groups is large enough to force, with the help of 
temperature, some of its dislocations through a Cottrell barrier. This would 
extend the slipped region around a source S further than the parallelogram 
of fig. 2, and give rise to a hardening rate lower than the linear stage. 

As an experimental confirmation of this point of view, we may refer 
to Diehl, Mader and Seeger (1955), who observed on several metals 
(Al, Cu, Ni) extended and multiple slip bands appearing in the parabolic 
stage, in contrast with the fine slip of the linear stage. More systematic 
observations of the Blewitt type (1955) would certainly be welcome. 

There are of course several possible geometrical configurations for 
the piled up groups and the barriers, sc various possible mechanisms 
for breaking through a barrier. Figure 4 represents two possible 
configurations. 

Fig. 4 


cp) 


! 


] 


(110) I ai 
a b 


Two possible geometrical configurations for a group of dislocations piled up in 
front of a Cottrell barrier. 


A complete discussion will be given by Stroh (1955). We shall only 
describe here a mechanism connected with the configuration of fig. 4(a), 
which we believe to be the most favourable for breaking through. 

Let us assume that, by thermal agitation, dislocation J has recombined 
with the barrier over a certain length h. Then the whole barrier has 
recombined over this length; for the upper half dislocation J’ has a 
Burgers vector (6/34/2) [112] which is now normal to the Burgers vector 
(b4/2/3. [111] of the rest of the barrier, and is therefore no longer repelled 
by it (fig. 5). 

Now, if @ is larger than the critical length of a Frank—Read source 
under an applied stress no, that is if * 
pb h 


no >= In 


nee 6 
2a Kh bo” (6) 


*h, may actually be somewhat smaller, as will be explained by Stroh (1955). 
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with b&b and K=0:8 for the orientation of the barrier (cf. Cottrell 1953), 
the recombined portion of the barrier, with a Burgers vector (b//2) [110], 
can bend in its (001) glide plane into a loop able to expand. There is 
therefore a stable final state where the dislocations have broken through 
the barrier, and which can be reached by thermal activation.t 

The activation energy is given by the difference of energy between the 
states of figs. 5(a) and (c). It is the sum of the energy U, to recombine 
the barrier over a length h, and that U, to produce a ‘ constriction ” such 
as in fig. 5(b). According to Stroh (1954, 1955)s 


Fig. 5 


A 


(III) 
ris B 
a b c 


Example of Cottrell barrier : (@) normal state ; (b) constricted ; (c) recombined. 


,  BOPh dl’ .: 
Gy= rns —2) ae re rk 
P ee. d’ 1/2 
and U,=0-9X 10-2u0a(In 5) : =< Le Ce 
0 
b? 1 
where y= eta SSS TES 
367(1—v) yt+nobWV/2 = = ~ = — * (9) 


is the width of the barrier, vis Poisson’s ratio and y the stacking fault energy. 
Now it is reasonable to assume that the expression 


: h d’ 
X = (in ) (wm by — 1) 


is positive and not much smaller than unity. Then, using (6), (7), (8) 
and (9), one can see that the constriction energy U, is always negligible 
compared with U,. The activation energy is therefore given by Uj. 
According to (6), it should be roughly proportional to (nc)-}, thus to 
o*; for, from § 4, x is proportional to c=c,+0,': 

2L?o 


eb Di 


- Ff 
c4x 10°. ee REY 


+ The loop of recombined barrier will actually soon split again under the 
applied stress. But the new barrier, having left the (111) glide plane of the 
piled up group, will no longer retain its dislocations. 
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For tractions at constant strain rate, the parabolic hardening should start 
at a constant value of U,/k7', and therefore at a stress og such that og2T 
as a constant. Table 1, deduced from Blewitt’s measurements on copper 
1955), shows that this is approximately the case. | 


Table 1. Relation between Temperature 7’ of Traction and the Stress ore 
at which Parabolic Hardening occurs (from Blewitt 1955, on copper) 


LK oc dynes cm-2x 10-8 Cpl Oso LO 
4-2 >20 ling) 
78 6-5 3:3 
200 4 3°2 
300 3 21 


The value observed for o,?7' in these measurements agrees at least 
roughly with the preceding equations. These give 
(In h/bo)(In d'/bp—1)bD (wu \?_,, 
meemnosskeizamee\G fh ce MD 
Thus, for copper and with the values L~50D~10-3 cm of § 4, U/kT is 
of the order of 3X 10!%X/o?7. One can expect for copper values of 
X =(In h/b,)(In d'/by>—1) between 1 and 5. With o,?7'~3x 10! c.g.s. 
from table 1, this gives a Boltzmann factor exp (—U/kT') somewhat 
too large, but of the right order of magnitude. 
Finally eqns. (9) and (11) show that U does not depend very much on 
the stacking fault energy y if nob>y, thus, from (10), for 
o? > Duy/2L? = 10 py em. Pres oe (12) 
It is easy to see that this condition is fulfilled over the range of stresses 
of the linear stage BC, except at its very beginning. Thus for copper 
(y=40 erg cm”) and aluminium (y@200 erg cem~?), (12) gives o>0-15 
and 0-25 108 dynes cm? respectively ; these stresses are somewhat 
smaller than that o, at the beginning of the linear stage (fig. 1). 


a 


Table 2. Values Observed for o2T'(ub)? for Various Metals 


Metal T°K | og dynes cm-2x 10-8 | [o¢?7'/(ub)?] c.g.s. x 10-7 | 
Al 90 4 2-4 , 
(300 ~0-2 Ort 
Ni |. 300 ai S22 
Cu 300 3 1-9 
Ag 300 2 2-5 
Au 300 1 0-6 
ee ee at eg i ee ee 


According to (11), the values of U and of o°T' should therefore vary as 
(ub)® from metal to metal, except at temperatures T high enough for the 
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linear stage BC to be short (gap). Table 2, deduced from measurements 
by Sachs and Weerts (1930), Staubwasser (1954), Kocks (1954), Blewitt 
(1955), Diehl, Mader and Seeger (1955), shows indeed values of o977'/(.b)° 
of the same order of magnitude for various metals, except for room 
temperature measurements of gold and aluminium, where deviations are 
expected because the linear stage is short. The value of » used in the 
table is that observed for polycrystalline material. 

Two other processes may be more favourable in some extreme cases 
(ef. Stroh 1955) : 

1. For metals with large stacking fault energies and at low stresses, 
that is at high temperature, the barrier of fig. 4(b) could recombine over 
a certain length and slip as a whole [110] dislocation, with a (100) slip 
plane. This could explain the (100) active slip planes observed in 
aluminium at high temperature (Schmid and Boas 1936). 

2. For metals with low stacking fault energies and under high stresses, 
the free half dislocations of the barrier in fig. 4(6) may be driven away 
from the barrier, thus extending the stacking fault and giving rise 
eventually to twinning. This could provide a mechanism for the forma- 
tion of numerous stacking faults (or thin twins ?) in some alloys with 
low stacking fault energy (Barrett 1952), and the deformation twins 
apparently observed in copper after fracture (Kuhlmann-Wilsdorf and 
Wilsdorf 1953) or by traction at very low temperature (Blewitt 1955). 


§ 7. PARABOLIC HARDENING 

We do not attempt a theory of the parabolic hardening CD (fig. 1), 
and restrict ourselves to the following remarks. 

1. There is an obvious relation between the traction experiments 
discussed in this paper and the yield points observed in ‘ strain recovery ’ 
by Cottrell and Stokes (1955). These authors strained a crystal in the 
parabolic region at a temperature 7’, after a prior deformation at a lower 
temperature 7',. The stress o, required for the deformation at 7’, was 
lower than the stress o, at the end of the prior deformation at 7',, and 
dropped further during the deformation, giving a stress-strain curve 
with an initial yield point. 

This behaviour could have been expected : the piled up groups formed 
under a stress o, at temperature 7’, start to break down, at a higher 
temperature 7’,, under a lower applied stress o,. A still lower stress is 
needed to develop the loops, once they have broken through the barriers, 
hence the yield point. 

2. Once dislocations begin to break through a Cottrell barrier, the 
loops which have escaped are too big to have much chance of forming 
a continuous Cottrell barrier again. We can therefore expect them to 
extend through a large part of the crystal. Their back stress on their 
emitting sources being negligible, these are expected to go on emitting 
new loops. Finally these big loops, when gliding in a slip plane near 
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enough to another piled up group, may help it to break through its 
barrier and develop into another big slip band (fig. 6). This would 


Fig. 6 


epee sp rei hy 
Production of multiple slip bands in the parabolic stage. 


explain the few, big, extended and multiple slip bands characteristic of 
the parabolic stage (Diehl, Mader and Seeger 1955). For slip lines 
extending over about #=1mm, we expect a minimum distance h 
between the components of a multiple slip band of the order of 

1 ie 1238 
Vie 
that is about 100A for a Frank—Read net of size 1210-3 cm. This agrees 
with the distances observed in aluminium (Heidenreich and Shockley 
1948, Diehl, Mader and Seeger 1955). It is to this stage that the hardening 
mechanism by cutting through the forest of screw dislocations discussed 
by Seeger (1955) can apply. 

3. From eqn. (10), the number 7 of dislocations in the piled up groups 
for monocrystals* cannot be larger than 4 10%og/u, if o, is the critical 
stress where parabolic hardening starts at the temperature considered. 
For aluminium and copper at room temperature, the values of og given 
in table 2 correspond respectively to n=250 and 2500. 

Now, in recovery and high temperature creep, the stress dependence 
of the activation energy may be related to the value of n, if these processes 
are regulated by the climb of dislocations in front of piled up groups 
(Mott 1953). The energy of activation for recovery and steady state 
creep can then be of the formt 


U’=const.—nob?3, eee ee SRY) 


where nab? is the work done by the stress no in front of a piled up group 


h= 


* In polycrystals, larger values of m can of course occur in groups piled up in 
front of grain boundaries. 

+ This formula assumes nob? >kT. For o> 10-4 and as pb? ~5 ev in alumin- 
ium, this gives n240 at room temperature. 
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when a jog of length 6 climbs over a distance 6. Values n=100 to 200 
in aluminium can indeed be deduced from recovery measurements on 
monocrystals by Kuhlmann, Masing and Raffelsieper (1949, 1950); and 
also from an analysis of 8 (transient) creep measured by Wyatt (1953), 
if one assumes with Mott (1953) an activation energy 4U’ for that type 
of creep. These last measurements were made on polycrystals with 
large grains. It would be interesting to repeat them on monocrystals 
and with other face-centred cubic metals, such as copper and nickel. 

As, according to (10), is proportional to o, the activation energy U”’ 
should vary linearly with o?, not c. The measurements made on recovery 
and creep are, however, not accurate enough to decide which law is 
observed. 


The author wishes to thank Professor N. F. Mott for his constant 
help and Drs. A. H. Cottrell, B. Jaoul and A. N. Stroh for many 
enlightening discussions. 


A PP Ng D Ax A. 


FORMATION OF COTTRELL BARRIERS IN FACE-CENTRED CUBIC CRYSTALS 


It seems reasonable to assume that a Cottrell barrier can be formed 
readily only if it is stable and if the two component dislocations D,, D, 
attract each other at large distances. 

We shall take the two dislocations parallel to the interaction of their 
glide planes. And we shall assume their energy of interaction proportional 
to the scalar product of their Burgers vectors and to the logarithm of 
their distance. We neglect therefore the differences between screw and 
edge dislocations, and also the angular dependence of their interaction. 

We shall then show that a screw dislocation cannot form a Cottrell 
barrier readily ; and that an edge dislocation belonging to a given ‘ half’ 
slip system (i.e. a system with a Burgers vector of given sign) can form 
readily a barrier with a given orientation only with edge dislocations 
belonging to another definite half slip system. The density of S’ sources 
is therefore equal to that of S sources, as stated in § 3. 


A.1. Short Range Interactions 


Let (111) and (111) be the two glide planes. By symmetry, it is 
sufficient to consider for D, a screw* [330] and an edge [033] directions. 
We consider for D, the six possible directions of the Burgers vector. 

D, and PD, are split into half dislocations X, Y separated by ribbons 
of stacking faults. The splitting can occur in two ways, leading to 
stacking faults of different natures (of the type ABC/BCAB and 
ABC/B/ABC) and half dislocations with different core energies (Frank 
and Nicholas 1953). It is reasonable to assume the difference of their 


* Burgers vectors are given here in units of b/34/2. 
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energies large enough for one kind of splitting to occur much more 
frequently than the other. This splitting, obtained following Thompson’s 
rule (1953), is given in table A 1. 


Table A 1. Decomposition into Half Dislocations X, Y 


330 033 330 303 033 330 303 033 


x 121 12] 211 211 112 121 oe 121 


ue 21 A Wie ey 112 121 211 211 {12 


The split dislocations can have one of the four dispositions pictured in 
fig. Al. The interaction between the two nearest half dislocations becomes 
preponderant when D, and PD, are near to each other. It must therefore 
be considered first. 


y! ! 
\ X5 Y2 ' x! 


Four possible cases for the interaction of dislocations D, and Dg. 


(a) If these two half dislocations repel each other (cases denoted r, 
table A 2), there is no barrier formed. 


(6) If they attract each other, the complex dislocation thus formed may 
repel the two remaining half dislocations sufficiently to form a stable 
barrier. These cases are denoted in table A 2 by the Burgers vector of 
the complex formed. They are of course linear combinations of the — 
dislocations of the ‘ stair-rod ’ type introduced by Thompson (1953). 


(c) If the complex does not make a stable barrier with the two 
remaining half dislocations, it absorbs them to make an ordinary dis- 
location, denoted in table A 2 by its Burgers vector. This vector is 
usually in one of the initial glide planes, in which case the dislocation 
will split into half dislocations and be able to glide, thus providing no 
barrier. In the case of two screws, there is no dislocation left. Finally 
the combination of [303] with [033] leads to dislocations with Burgers 
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vector not in (111) ) nor (111); these will split immediately into the more 
stable complex [110] and two half dislocations, formed directly in one 
of their dispositions (table A 2). 


A.2. Long Range Interaction 
The splitting of the dislocations may be neglected when they are at 
large distances. Depending on the angle between their total Burgers 
vector, they repel each other (R, table A 3), attract each other (A) or - 
do not interact (N). 


Table A 3. Long Range Interaction between two Dislocations 


The following conclusions may be deduced from tables A 2 and A 3. 
1. Barriers cannot be formed with two screw dislocations. 


2. A screw and an edge dislocation can form a (not very stable) barrier 
by short range interaction. This is however very unlikely to occur, for 
the dislocations repel each other at large distances. 


3. With a given edge dislocation ([033] for instance), half of the possible 
edge dislocations in the corresponding intersecting glide plane can form 
a barrier of a given direction (here [110]) by short range interaction ; 
and these edge dislocations belong to three ‘ half’ slip systems (i.e. 
systems with a Burgers vector of given sign). But only a quarter, 
belonging to one half slip system (here [303]) are attracted at long distance 
by the given edge dislocation, and are thus likely to form a barrier. 
This barrier is of the type studied by Cottrell (1952). It has the smallest 
possible Burgers vector, thus is the most stable. 
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SUMMARY 


Changes in the flow stress of copper crystals due to changes in the 
temperature of deformation have been measured over a range of strains 
and temperatures. After correction for the effect of temperature on 
elastic constants the flow stress proves to be insensitive to temperature 
from 90°K to 180°xK, falls rapidly with increasing temperature from 180°K 
to 250°K, and is virtually constant from 250°K to 473°K (the highest 
temperature examined). The change of flow stress with temperature is 
closely proportional to the flow stress itself. These results are discussed 
briefly in terms of recent ideas of dislocation processes sensitive to tem- 
perature. The effect of increasing the temperature of deformation is to 
produce a yield drop, similar to those observed in aluminium crystals 
during work softening. 


§ 1. INTRODUCTION 


Ir has been realized in recent years that in addition to the elastic inter- 
actions between parallel dislocations (Taylor 1934) several other processes 
involving dislocations may contribute to work hardening (Mott 1952, 
1953); for example the coalescence of dislocations to form sessiles (Lomer 
1951, Cottrell 1952 a), the formation of jogs on intersecting dislocation 
lines (Heidenreich and Shockley 1948, Cottrell 1952 b, Seeger 1954 a and b), 
and the creation of vacancies and interstitials at jogs (Mott 1951, Seitz 
1952, van Bueren 1953, Seeger 1955). It is therefore important to develop 
experimental methods for separating and identifying the respective 
contributions of these processes to work hardening. One possible method 
of attack, which was chosen for the present experiments, is to determine 
the change of flow stress with temperature. This ought to provide 


* Communicated by the Authors. 
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information about the contributions of the processes according to their 
sensitivity to temperature. In these experiments single crystals of copper 
were used; the results of similar experiments on aluminium crystals 
have recently been described (Stokes and Cottrell 1954, Cottrell and 


Stokes 1955). 


§ 2. EXPERIMENTAL METHOD 


Wires of 1:5mm diameter were prepared by swaging rods of both 
spectroscopically pure (99-999°%) and oxygen-free high-conductivity 
(99-99%) copper. After pickling and coating with colloidal graphite, 
these wires were grown into single crystals 17 cm long by the fusion method 
of Andrade and Roscoe (1937), being supported horizontally during 
growth in loosely fitting and carefully cleaned silica quills. The crystals 
were cleaned in dilute nitric acid and cut into 6 cm lengths with a fine 
oxy-acetylene flame, water cooling being used to protect sections near 
the cut. Small iron loops which enclosed 1 cm lengths at the ends of the 
crystals were then attached by silver solder, thus leaving a working 
length of 4cm. The completed specimens were electropolished in a 
phosphoric acid cell and their orientations determined by the Laue back- 
reflection method (fig. 1). 


Fig. 1 


Initial Oricntations of the Specimens. Numbers 4 and 5 were made from 
oxygen-free high-conductivity copper; the remainder from spectro- 
scopically pure copper. 


Tensile experiments were made with a hard-beam machine at a strain 
rate of 4x 10-° per sec. A small load was always maintained on each 
specimen, between successive experiments, to preserve the alignment of 
the hooks and loops through which the load was applied. To maintain a 
selected temperature the mounted specimen was enclosed in a‘ thermal 
bath attached to the machine. A standard time of 15 min was allowed 
in each bath before the experiment was begun and about 2 min was 
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needed to replace one bath by another. The following temperatures 
were used :— 


Temperature °K Bath Fixed Point 
90 liquid oxygen boiling point 
196 solid carbon dioxide sublimation point 
and acetone 
233 solid and liquid calcium freezing point 
chloride solution 

273 water and ice freezing point 

iiss silicone oil controlled to +1°x 
473 silicone oil controlled to +1°K 


The general method of making experiments was to strain each crystal 
by repeated increments, each about 0-01, applied alternately at two 
temperatures. The load and elongation were recorded continuously and 
the effect of temperature on the flow stress was determined by comparing 
the flow stress of the specimen at the end of one increment with that at the 
beginning of the next increment. 


§ 3. RESULTS 


Two series of stress-elongation curves are shown in figs. 2 and 3, 
obtained from crystals of spectroscopically pure and oxygen-free high- 
conductivity copper, respectively, cycled between 90°K and 473°x. 
Figure 2 in addition shows two isothermal curves obtained from crystals 
of the same orientation as that used in the cycling experiments, taken at 
90°K and 473°K respectively. 

These isothermal curves show that at small strains more work hardening 
occurs at 473°K than at 90°K, as a result of the shorter period of easy glide 
at 473°K. At large strains, however, the crystal cannot maintain the same 
high level of flow stress at 473°K as it can at 90°K; the rate of work- 
hardening falls off at 473°K but not at 90°K. This latter effect, which 
‘has also been observed by Blewitt, Coltman and Redman (1955), seems to 
be related to the yield phenomenon which appears in the strongly worked 
crystal on straining at 473°xK after a preliminary straining at 90°K. Ina 
study of similar yield drops in aluminium crystals, Cottrell and Stokes 
(1955) showed that a cold-worked crystal softened extremely rapidly 
when strained again at a higher temperature than that of working. 
They suggested that this work-softening effect was caused by the collapse, 
due to the joint action of stress and thermal agitation, of sessile dis- 
locations created during the previous straining at the low temperature. 
This thermally activated collapse of sessile dislocations during plastic 
flow has also been suggested independently by Friedel (J. Friedel, private 
communication and Phil. Mag., 1955) to explain the observations of 
Blewitt, Coltman and Redman on the effect of temperature on the work- 
hardening of copper crystals at large strains. 

From the point of view of measuring the reversible change of flow stress 
with change in the temperature of straining, the yield drops which occur 
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after changing from low to high temperatures are an unwelcome com- 
plication since they add an irreversible change of flow stress to the 
reversible one. Changes from high to low temperature on the other hand 
do not suffer from this effect. Accordingly, all the measurements of 
changes of flow stress reported below have been obtained from transitions 
of this second type. 


— DEFORMATION ALTERNATELY AT 90°K & 473°K; SPECIMEN NP I. 


= —— DEFORMATION AT 473°K, SPECIMEN NV 2° 
----- DEFORMATION AT 90°K, SPECIMEN N® 3. 
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Stress-elongation curves on crystals of spectroscopically pure copper. 
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Stress-elongation curves on a crystal of oxygen-free high-conductivity copper 
(specimen number 4) deformed alternately at 90°K and 473°k. 


Changes in flow stress produced by changes from 473°K to 90°K made at 
various stages during the elongation of the crystal are given in the table, 
taken from the results in figs. 2 and 3. An important feature is that the 
ratio of flow stresses at these two temperatures settles down after an 
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elongation of 10-15% to a constant value that is practically the same for 
both grades of copper. This observation, which has also been made on 
aluminium (Cottrell and Stokes 1955), encourages the view the ratio of 


The Ratio of Flow Stress at 473°K to that at 90°k, Measured on Copper 
Crystals by Transitions from 473°K to 90°K at Various Elongations 


Specimen Number 1 


Per cent elongation : 10-0 12:5 14-5 Uizices A Shes 21:2 24-2 27-0 

Ratio of flow stress : 0-858 0-854 0-824 0-82 0:822 0-826 0:820 0-821 
Specimen Number 4 : 

Per cent elongation : 3°8 7:3 10°5 13-5 16-3 1:90 214 

Ratio of flow stress : 0-783 0-819 0-837 0°832 0-828 0:823 0°834 


flow stresses, measured in this way, may provide a clue to the nature of 
the elementary events involving individual dislocations ; the absolute 
magnitude of the flow stress, on the other hand, depends on the number 
of such events as well as on their nature. The ratios of flow stresses 
obtained from all the subsequent experiments were measured in this 
constant range. 


FLOW STRESS RELATIVE TO THAT AT 90° K. 


° 100 200 300 400 500 
TEMPERATURE, °K 


The reversible change of flow stress with temperature for work-hardened copper 
crystals, before (lower curve) and after (upper curve) correction for the 
change of elastic constant with temperature. The numbers of the 
specimens are shown, 


Figure 4 summarizes the results of these experiments and gives the flow 
stress at various temperatures relative to that at 90°K. All values except 
those at 195°K were obtained by direct transitions from the temperature 
concerned to 90°K. To avoid contact between liquid oxygen and acetone, 
the values for 196°K were obtained indirectly, by double transitions from 
273°K to 196°K and from 273°K to 90°K. The broken curve in fig. 4 shows 
the ratio of flow stresses directly. Since the elastic constants of the 
metal vary with temperature it is important to examine the variation of 
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o/u with tensperature, where o and y are the flow stress and shear modulus 
measured at each temperature. In the ‘corrected’ flow stress curve 
(full line) of fig. 4 the values of Késter (1948) for the variation of Young’s 
modulus with temperature have been used, since this variation is very 
similar to that of shear modulus. 

Correcting for the effect of elastic constants brings out the remarkable 
fact that the residual sensivity of the flow stress to temperature is con- 
fined almost entirely to a ‘ shelf’ from 180°K to 250°K across which the 
flow stress falls from near unity to 0-93. Outside this range of 
temperature there is very little variation with temperature; this 
constancy is particularly striking over the range from 250°K to 470°K 
where the corrected flow stress remains entirely inside the limits 0-924 
and 0-934. 


§ 4. Discussion 


The main processes suggested so far that may provide large reversible 
changes of (corrected) flow stress with temperature in a metal such as 
copper at low temperatures are (a) the intersections of dislocation lines 
gliding through a ‘ forest ’ of dislocations (Cottrell 1952 b), (6) the thermal 
creation of point defects at jogs on moving dislocations with screw com- 
ponents (Seeger 1955), and (c) the thermally activated glide of jogs on 
edge dislocations (Maddin and Cottrell 1955). A detailed interpretation 
of fig. 4 in terms of the first two of these processes has already been given 
by Seeger (1955). His theory, which we shall now outline briefly, 
depends on the differences in behaviour of edge and screw dislocations 
and on the principle that at any given temperature and rate of straining, 
the strain is produced mainly by the type of dislocation that is easier to 
move. At low temperatures (i.e. below 180°K) the edge dislocations 
move more easily than the screws because the latter must create vacancies 
(or interstitials) at jogs whenever they move. The moving edge dis- 
locations cut their way through the forest and this leads to a small 
thermal dependence of flow stress in this low-temperature range. Above 
this range it becomes easier to move screw dislocations than edges because 
thermal fluctuations can play an increasingly important part in creating 
vacancies at jogs in screw dislocations ; the steep fall of flow stress at 
200°K is attributed to this effect. Above a certain temperature the 
generation of vacancies at jogs on screw dislocations occurs so easily that 
it is no longer an appreciable source of resistance to the motion of screws. 
Accordingly, the strength remaining in this high-temperature range is 
due mainly to long-range elastic interactions of dislocations. Since 
these elastic interactions depend on temperature only through the 
elastic constants the (corrected) flow stress is constant at these 
temperatures. 

Neeger’s theory attributes the flow strength of copper in various tem- 
perature ranges to several distinct dislocation processes. Looked at 
from this point of view the results in the table are intriguing since they 


Temperature on the Flow Stress of Work-Hardened Cu Crystals 1193 


imply that over a wide range of cold-worked states the change of flow 
stress with temperature is closely proportional to the flow stress itself. 
The theory thus requires a proportionality to be maintained during work- 
hardening between the density of forest dislocations, the density of jogs 
on screw dislocations, and the density of internal stresses in the metal. 
Such a proportionality might be expected if the various distributions of 
dislocations produced in the present series of experiments differed only in 
scale, and not in pattern. Clearly, an obvious extension of these experi- 
ments (which it is hoped soon to make) is to measure the flow stress as a 
function of temperature in crystals containing radically different patterns 
of dislocations. 
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SUMMARY 

The present paper treats in some detail the various mechanisms by 
which lattice defects are formed by moving dislocations. Particular 
attention is paid to the role of thermal energy in these processes and to the 
special features of extended dislocations in metals with low stacking fault 
energies. The results are applied to the temperature dependence of flow 
stress of such metals and are found to explain the temperature dependence 
observed in Cu quantitatively. Cu is known to have a low stacking fault 
energy. 

It is possible to deduce details about the mechanism of plastic deform- 
ation, in particular of the role played by the edge and the screw parts of 
dislocation rings under various conditions. The predictions of the theory 
on the nature of the lattice defects causing the increase of electrical 
resistance during cold work are shown to be borne out by recovery 
experiments. 


$ 1. INTRODUCTION 


Ir has been known for some time that moving dislocations may generate 
lattice defects (for a review see Seitz 1952). If we confine ourselves—as 
we shall do throughout this paper—to temperatures below that of self- 
diffusion so that climbing of dislocations is unimportant, the accepted 
mechanism is that a jog in a moving screw dislocation will leave behind 
it a row either of vacancies or of interstitial atoms. Very little attention 
has so far been paid to the finer details of these processes and in particular 
to the influence of thermal energy on the mode of production of vacancies 
and other lattice defects. 

In § 1 we undertake a study of the various processes that lead to a 
generation of lattice defects by gliding dislocations. We shall investigate 
in particular the temperature dependence of these processes, giving special 
consideration to the conditions under which vacancies are generated 
thermally by jogs in extended dislocations.+ Most of the discussions in 
the present paper refer to elementary jogs (i.e. jogs connecting two pieces of 
a dislocation line in two neighbouring parallel glide planes). Sometimes 
composite jogs, which are for example formed in dislocations in the h.c.p. 
ee ee eee 


* Communicated by the Author. 
+ A more detailed description of jogs, their formation and annihilation has 
been given elsewhere (Seeger 1955 a) 
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lattice that cut through dislocations with a Burgers vector perpendicular 
to the basal plane (Seeger 1955 a), are algo discussed. Throughout our 
discussion we shall confine ourselves to unit dislocations in the esc 
b.c.c. or h.c.p. structures, thereby avoiding the complications that arise 
in the case of dislocations with large Burgers vectors. 

The application of these results to the temperature dependence of the 
flow stress in f.c.c. cubic metals is based on the author’s earlier discussions 
of the temperature dependence of plastic properties of metals crystallizing 
in close-packed structures (Seeger 1954 a, 1955 b). In these papers a 
detailed account of the temperature dependence of the flow stress of 
metals with high stacking fault energies was given. The present treat- 
ment (§ 3) of metals with low stacking fault energies (in particular noble 
metals, Ni and Pt) completes the general discussion (Seeger 1955 b) and 
supplies numerical values for some of the activation energies involved. 
Very good agreement with experiment is obtained for Cu, the only metal 
for which suitable experimental results are available. Qualitative 
agreement is reached for Ni. The implication of this is that Ni must also 
be considered to have a low stacking fault energy. As in one of the 
author’s earlier papers (Seeger 1954 a) the (not very extensive) experi- 
mental results on Ni had been tentatively interpreted on the assumption 
of a high stacking fault energy we discuss in Appendix A the present 
evidence on the stacking fault energy of Ni. 

The good agreement between theory and experiment for Cu allows us to 
draw further conclusions on the mechanism of plastic deformation at low 
temperatures of this metal, which is believed to be representative of the 
noble metals group. It turns out that at 180°K a transition occurs from a 
low temperature mechanism in which the edge dislocations move more 
easily than the screw dislocations to a mechanism in which the screws move 
more easily than the edges so that the strain rate imposed in a tensile test 
is maintained mainly by the motion of dislocations with a large screw 
component. This change in mechanism is accompanied by a change of the 
nature of the lattice defects created during deformation (e.g. vacancies, 
pairs or clusters of vacancies, and interstitial atoms). As these lattice 
defects are known to be responsible for a large portion of the change of 
electrical resistance during cold work it should be possible to check these 
predictions by recovery experiments, and the experimental data confirm 
the theoretical expectations. An annealing stage at about 100°c, which we 
believe to be due to the migration of isolated vacancies, is pronounced only 
if the specimens were cold-worked at sufficiently high temperature. In 
specimens deformed at very low temperatures, in which according to the 
theory very few isolated vacancies are generated, this annealing stage is 
found to be weak (§ 4). 


§ 2. ATHERMAL AND THERMAL GENERATION OF Lattice DEFECTS 


Thermal activation may or may not be important in the production of 
lattice defects by moving dislocations. As the effect on plastic properties, 
e.g. ‘on the temperature dependence of flow stress, will be different for 
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thermal and athermal generation of lattice defects we start by reviewing 
the different ways of creating defects by moving dislocations and discussing 
the role of thermal energy in each of them. 

For simplicity we first consider pure screw dislocations with elementary 
jogsinthem We shall later return to the general case that the dislocation 
lines have both screw and edge components. 

If the shear stress + acting in the glide plane is sufficiently large, both 
vacancies and interstitials may be generated in strings (fig. 1). 7 has to be 
large enough to overcome the line-tension of the string. The energy 
necessary to add an interstitial or a vacancy to the end of the string may be 


er eten U,=«,G0° eS Se ors 
(G shear modulus). In typical metals (e.g. Cu) «, is of the order of unity 
for interstitials and about 0-1—0-2 for vacancies. If the distance between 


neighbouring jogs is /,, we obtain from considerations of the work done by 
the shear stress during the generation of a defect (Cottrell 1952) 


aed O85 (aa fF a ore 
T2305 GUibys, 0. 2) a ee ee ee 


Equation (3 b) is also the equation for the critical stress of the Frank—Read 
source of length /, that extends between the points A and B in fig. 1. In 


or, by using (1), 


Fig. 1 


f, 


Generation of strings of lattice defects by a gliding dislocation of screw character. 
A and B are the position of the jogs in the dislocation line. 


this case the numerical factor «, in eqn. (2) is about unity and will be 
denoted by «,. 
We can now state the conditions under which (in the absence of thermal 


activation) defects are continuously added to the ends of the rows of 
defects shown in fig. 1: 


(i) The shear stress 7 must be equal to or slightly larger than the shear 
stress given by eqn. (2 a). 


(ii) The shear stress 7 must satisfy the inequality 


7 <i bl sedge ahi aes 
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There are many cases in which it is possible to satisfy the conditions (i) 
and (ii) simultaneously, e.g. the production of vacancies by elementary 
jogs in f.c.c., b.c.c., and h.c.p. lattices. In the case of composite jogs the 
situation is more difficult, but it seems safe to assume that in the simple 
cases just mentioned our conditions can still be fulfilled for jogs composed 
of two elementary jogs. 

The situation is much more critical if the direction of motion of the jogs 
is such that they will produce interstitial atoms. Since for the generation 
of one interstitial at a time «, is about equal to «,, the satisfying of our 
conditions depends on the metal under discussion. If the jogs are 
composite ones the critical shear stress «, of the Frank—Read-sources, given 
by ee AC ee meme OE ee ay 
will be reached before the shear stress is large enough to add new groups of 
interstitials to the end of the strings. The dislocation lines between A and 
B will bulge out, sweep part of the glide plane and unite as shown in fig. 2. 


Fig. 2 


Se 


Generation of a series of lattice defects due to the action of the Frank—Read 
source of length J, with anchor points A and B. 


The attraction of the edge dislocations 1 and 2 will certainly be strong 
enough to annihilate the dislocations, the annihilation products being 
interstitials. It requires, however, a more detailed investigation to decide 
whether the local heating during the annihilation process is large enough 
to allow the interstitials to diffuse away (see Seitz 1952). If the inter- 
stitials do not diffuse away the situation is the same as if the interstitials 
had been created consecutively. | 

Thermal activation is not important in either of these processes for the 
production of lattice defects, unless the temperature is high enough to 
allow the defects to diffuse away from the strings. As pairs and clusters of 
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vacancies almost certainly have a smaller activation energy of migration 
than a vacancy, a string of vacancies will presumably break up into 
groups of vacancies rather than isolated vacancies, either immediately 
after formation or during an anneal. The essential feature in both cases 
discussed above is the balancing of the line tensions of the dislocations and 
the strings of lattice defects by the shear stress + acting on the dis- 
locations. Thermal energy aids the shear stress only to the extent to which 
it lowers these line tensions, which are free energies per unit length and 
depend therefore somewhat on temperature. 

If we consider dislocation lines with an increasing ratio of the edge 
component to the screw component, the distance between neighbouring 
defects in the strings produced by the forward motion of jogs will increase 
until we can regard the defects in a string as independent of one another. 
Then in eqn. (2) U, has to be replaced by the energy U, necessary for 
creating isolated defects. We can write 


. Uyga, Gb. 2. 8 ee 
The remarks on the magnitude of «, also apply to a». 

At temperatures low enough for thermal energy to be unimportant the 
situation is much the same as with pure screw dislocations. Depending on 
whether 

ga, MP Be a 
or Oe sacs . hace be rs 


the defects will either be produced consecutively or many at a time 
through the action of the dislocation source and subsequent annihilation 
of the dislocation lines. 

If condition (6 @) is satisfied the position is quite different if thermal 
energy is sufficient to aid the production of vacancies by elementary jogs. 
This may be envisaged as occurring in two ways. 

(a) In creep it has been suggested by Mott (1955) that a jog, momen- 
tarily at rest, can move forward one atomic distance with the help of the 
stress, the vacancy (or interstitial) moving away at the same time. The 
activation energy for this will be 

W-+U,—7b7l, 
where W is the activation energy for the movement of a vacancy, so that in 
f.c.c. metals W-+-U, is that for self-diffusion. 

(b) In reasonably rapid extension, however, as considered in this paper, 
we suppose that the jogs in a screw dislocation (which after all are nothing 
but short edge dislocations) will continuously move up and down the 
dislocation line under the influence of the applied stress and the rapidly 
fluctuating internal stresses due to work-hardening. A jog can therefore 
jump forward while in motion sideways and continue its motion sideways. 


The hole thus formed does not have to diffuse away to be stable, so the 
activation energy for this process is 


U,—7b4l;. 
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This process is the one considered in this paper and we shall return to it in 
§ 3. 

We now discuss in more detail this possibility for the jog to move side- 
wards, i.e. along the dislocation line. For jogs in a pure edge dislocation, 
sideward motions of the jogs are always non-conservative and therefore 
unlikely to occur at very low temperatures. On the other hand, in a pure 
screw dislocation the sideward motion of a jog is conservative, whereas the 
forward motion is non-conservative and leads, as discussed above, to the 
production of lattice defects. This general rule is that a jog in a disloca- 
tion line with Burgers vector b will move conservatively in the direction of 
b. Ifb is not exactly parallel to the dislocation line (i.e. the dislocation 
has an edge component), even a dislocation with jogs may glide without 
the production of lattice defects, the jogs moving conservatively along the 
dotted lines shown in fig. 3. 


Fig. 3 


~ 


SS 
= conservative motion of jog <\ 
ok Tia Se motion of 
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she dislocation line 
c SS 


Conservative motions of jogs A and B in a dislocation line. 


In a tensile test an unbalanced force due to the applied shear stress and 
the stress-fields of the surrounding dislocations rearranging themselves 
during the work-hardening will in general act on a jog in a screw disloca- 
tion. If the dislocation itself is held up at some obstacle this force will 
cause the jog to move along the dislocation line. Figure 4 shows an 
example in which the area in the glide plane swept out by a dislocation 
whose mean character is that of a screw can be increased by the movement 
of a jog. But even if the screw dislocation containing the jog is straight 
(apart from the jog) work will be done by the shear stress during the motion 
of the jog. This is so because the force on a dislocation line will in general 
have a component perpendicular to the glide plane. Since by the move- 
ment of the jog the dislocation is lifted from its glide plane to a neigh- 
bouring plane this force causes to move the jog along the dislocation line:* 
In this mechanism jogs of opposite sign move in oppcsite directions. If 
a pair of jogs of opposite sign meet each other they will recombine, thereby 
reducing the number of jogs in screw dislocations. Irrespective of the 
distance travelled by a screw dislocation through the forest an average 
distance /, between dislocations corresponding to a dynamical equilibrium 


* It is assumed that the screw dislocation cannot yield by cross slip to the 
force acting on it. This assumption is satisfied in, metals with low stacking 
fault energies, to which the present considerations will be applied in § 3. 
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will be established. J, will depend on such parameters as the density of 
the forest of dislocations threading the glide plane, the velocities of both 
dislocation lines and jogs, and the overall strain rate «. As the motion of 
jogs along the screw dislocations is envisaged not to depend on thermal 
activation, temperature will influence J; only indirectly through the varia- 


tion of these parameters. 
Fig. 4 


1% 


—— 


Un slipped area 


Slipped area 


conservative motion of the jog 


The shear stress, tending to increase the slipped area in the glide plane, will 
’ cause the jog A to move conservatively to the left. 


Fig. 5 


X,  X3X5 x 


oa versus distance diagram for a dislocation line containing a jog 

Juring the re x ] The 

ve ing the movement from #, to x, the jog creates a vacancy. The 
maximum of the potential energy is reached at the position a. The 
applied shear stress is 7. < 


2 . 

Since the forces per unit length acting on the jog are of the same 
order of magnitude as those acting during the glide motion of the dislo- 
cations the usual arguments (see e.g. Nabarro 1951) can be applied as to 
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ie velocity of a jog. The velocity attained by a jog under these forces 
should be of the order one-half of the sound velocity or slightly less. 
As we shall see in the subsequent discussion, velocities of this magnitude 
are sufficient to allow the production of isolated vacancies by the motion 
of jogs. 

Figure 5 shows the diagram of energy as a function of distance for a screw 
dislocation which contains a jog. For the production of a vacancy by a 
non-conservative movement of the jog we adopt the following model : 
the vacancy is created during the forward motion of the dislocation 
line by a distance x,—x,=b. The maximum of energy during this move- 
ment is reached at a position x, given by 


Cet | Boa ho— 0/4. 2 . » . . ee (7) 


With a shear stress + acting on the dislocation the activation energy for 
motion in forward direction is 


Ces pee Si brite ee le, ee (8) 


The activation energy for the jog to jump backwards, thereby recombining 
with the vacancy, is 
O=1,67 7/4. Meee es re” (9) 


Equation (9) shows that in order to prevent a substantial fraction of the 
jogs from recombining with the vacancy 7 must not be too small. Other- 
wise eqn. (8) would have to be amended by a term taking into account 
this possibility. We shall base our discussion on fig. 9, which applies to 
Cu. The full line denoted by 3 corresponds to the process under discussion. 
The lowest stress 7, for which eqn. (8) is used in that diagram, is that at the 
temperature 7’, and is given by 


1 jes 
p= —7, === et soe 7.88 0) 
47% ee cy 
(for further details of notation see § The frequency of backward jumps 
18 
ee Li) 


which, with ,=0-86 ev and 7’,=250°K, gives 
y=0-05 v9. Peer 11) 


vp Will be of the order of Debye’s frequency. Therefore the velocity of the 
jog along the dislocation must be larger than 0-05 of the sound velocity if 
the probability of jumping backward and recombining with the vacancy 
is to be neglected. On the basis of the discussion given above we believe 
that this condition is satisfied in tensile tests of typical metals. 

If the conditions just discussed are fulfilled a jog in a screw dislocation 
will carry out a zig-zag motion, consisting mainly of a rapid movement 
along the dislocation line, from time to time interrupted by thermally 
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activated steps perpendicular to the Burgers vector, in which isolated 
vacancies are created, and possibly reversals of the direction of motion. 
This is shown schematically in fig. 6. The jog’s motion will end in a 
recombination with a jog of opposite sign. 


Fig. 6 


motion of dislocation 
t line 


Zig-zag motion of a jog Ain ascrew dislocation The dotted lines represent the 
conservative movements of the jog, the circles represent the lattice 
defects created in non-conservative motions. 


If in the process just discussed a complete separation of the defect and 
the dislocation line occurs the activation energy of the process is given by 
the energy U, of the free defect minus the work done by the applied shear 
stress. This has been made use of in eqn. (8). If the dislocation involved 
has exactly screw character it is, however, likely that the defect will 
remain attached to the dislocation line immediately after its creation. 
The activation energy will be lowered by the ‘ binding energy’ U, and 
will be given (the influence of the applied stress 7 being neglected) by 


G.30,—U,. ~. ae 


Having in mind the application of these considerations to the flow 
stress of close-packed metals, we consider in some detail the case of a 
f.c.c. or h.c.p. structure in which there is an appreciable separation of 
complete dislocations in close-packed glide planes into partial dislocations. 
At the position of a jog the partial dislocations of the extended dislocation 
have been brought together, forming a constriction (Stroh 1954, Seeger 
1954 ¢, Schoeck and Seeger 1955). If the jog and therefore also the 
constriction move away from a defect which has just been generated, the 
defect will remain attached to one of the partial dislocations as shown in 
ee 

The energy U', by which an interstitial or a vacancy is bound to a dis- 
location in a f.c.c. metal consists in general of two contributions, namely 
an interaction energy L’,, due to the electric interaction between the lattice 
defect, and an energy term U7,;, which arises from the elastic interaction of 
the lattice defect. 
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An edge dislocation in a metal containing ‘ free ’ electrons will form an 
electrical line dipole (Landauer 1951, Cottrell, Hunter and Nabarro 
1953), the compression side being positively charged and the dilatation 
side being negatively charged. A lattice defect in a metal such as copper 
can be looked upon as a negative elementary charge which is screened by 
the electrons in the Fermi gas (Fumi 1954). Due to the spatial variation 
of the electric field around a dislocation a vacancy on the compression side 
of the dislocation will be repelled if sufficiently far away from the centre of 
the dislocation. This results from the fact that the attraction of the 
negatively charged centre of the vacancy is more than overbalanced by the 
repulsion of the positively charged screen. If, however, the vacancy is so 
close to the dislocation line that the screen overlaps appreciably into the 
dilated region around an edge dislocation, a net attraction results. At a 
distance of about 1 4 from the dislocation centre, where a maximum of the 
elastic interaction is to be expected, the two opposing effects discussed 


Fig. 7 
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b dislocation line 
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A jog (A) in an extended dislocation having created a lattice defect moves 
away from the defect (D) which remains attached to one of the partial 
dislocations. 


above cancel, so that for the present problem the electrical interaction 
between dislocations and vacancies may be neglected. It should be 
pointed out, however, that these results are based on a classical calculation 
of the interaction between the charge distributions around a vacancy and 
around a dislocation and that a more refined quantum mechanical 
treatment might change some of the details. For impurity atoms that do 
not cause a lattice distortion the elastic interaction can be neglected in 
comparison with the electric interaction. In such a case the foregoing 
considerations apply fully, whereas in the case of vacancies in Cu and 
probably in all noble metals only the elastic interaction needs to be 
considered. 


4M2 
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In f.c.c. structures the lattice distortions due to both interstitials and 
vacancies show cubic symmetry. We may therefore to a good approxi- 
mation use the formula (Bilby 1950) 

a= 4G0r,” SU ee 

to calculate the elastic energy of interaction between a dislocation and a 
lattice defect whose position with respect to the dislocation line (edge- 
dislocation of dislocation strength 6) is fixed by the polar coordinates 
R,¢. 7) is the atomic radius, and 6 relates 7) to the apparent radius 

r=(1-+-8) ry 2) eee 
of the defect. From the Born—Mayer potential proposed by Huntington 
(1953) for the interaction of the ion cores of Cut, a value for 6 of 0-0175 is 
easily deduced (Fumi 1955). The modulus of b sin ¢/R will be assumed to 
be equal to unity. A change in this term by a factor of $ would change the 
tinal result only by 5°,,. For want of better information we use the same 
numerical values also for Ag, Au, Ni and Pt. The interaction energies 
Uy, are given in the table. 


Data referring to the Generation of Lattice Vacancies at Jogs in 
Extended Dislocations in various f.c.c. Metals 


Metal Cu Ag Au Ni Pt 

Energy U, of a free 

vacancy (ev) 0-90 0-8 0-67 1-6 1-2 
Shear modulus 

G (kg/mm?) 4400 2700 2800 7600 6600 
Atomic radius 7'9(A) 1-275 1-444 1-442 1-245 1-378 
Uas=U, (ev) 0-04 0-03(5)| — 0-03(6)| —-0-06(3) 0-07(5) 
U,=U,—U;, 0-86 0-7(6) 0-63 1-5 iT 
SV/ke 3-4 (2-5) 1 (4) (1) 
T,‘*) (theoretical, °K) 364 (330) 290 (420) (530) 


Values in brackets less certain than others. 
For significance of 79‘) and comparison with experiment see § 3. 


The sources of the energies U, of formation of free vacancies used in 
the table are stated in Appendix B. The activation entropies S) for 
the formation of free vacancies for Cu and Au are deduced from experi- 
ment (see § 3), those for Ag, Ni and Pt are estimates. 


§ 3. TEMPERATURE-DEPENDENCE OF FLow Stress tn CLOsE-PackED 
METALS 
In f.c.c. and h.c.p. metals the main processes which contribute to the 
temperature-dependence of the flow stress below the temperature of self- 
diffusion are the intersection of dislocation lines and the production of 
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lattice defects with the aid of thermal energy. The mechanism of the latter 
process has been discussed in § 2. The Peierls force which would also give 
rise to a temperature dependence of the flow stress is thought to be 
unimportant for f.c.c. and h.c.p. metals with the possible exception of 
extremely low temperatures. We shall take into account the following 
three processes, which will be denoted by superscripts (1,) (2) or (3h: 
(1) Cutting of edge dislocations through the forest of dislocations threading 
the glide plane, (2) cutting of screw dislocations through this forest, 
(3) thermally activated generation of vacancies by jogs in screw disloca- 
tions. They have this in common, that the activation energy of the process 
depends linearly on the shear stress 7, on the dislocation.* The reason for 
this simply is that during the activation the dislocation is pushed forward 
through an ‘ activation-distance ’ d’ and the work done during activation 
is proportional to d’7,._ In this context 7, is to be understood as the shear 
stress effective on the dislocation line. It is given by the difference of the 
applied shear stress t and the shear stress tg due to other obstacles to the 
dislocation movement such as the elastic stress field of neighbouring 
dislocations. stg depends on temperature only indirectly through the 
shear modulus G. 

If only thermally activated processes of one of the two types just 
mentioned occur the strain rate is given by 


é=NAbvy exp Go ais “a ee ena 5) 
(Seeger 1954 a, 1954 d). WN is the density of sites at which the activation 
takes place, A is the area swept out by the dislocation during and after 
the activation process, vy is a frequency of the order of magnitude 
of Debye’s frequency, Gy is the free energy of activation at zero stress, 
and v is the ‘ activation volume’. For a dislocation of strength 6 the 
activation volume is given by 


(a0 A ae eee ee 82 (1.6) 


is the distance between the sites at which the dislocation lines are being 
held up. The activation distance d’ has already been introduced. If 
an isolated vacancy (or an interstitial) is created d’ is about equal to the 
interatomic distance b. The same is true if a dislocation cuts through an 
unextended dislocation. If the intersected dislocation is an extended one 
d’ will depend on the characters and the orientations of the dislocations 
involved. d’ may be larger or smaller than the separation 27) of the half- 
dislocations in an extended dislocation (Seeger 1955 b). In Cu, Ag, Au 
and Ni, 27 is about 10 b, and it seems reasonable to expect also d’~10 6 in 
these metals. 


* For process (3) this holds only if the probability of jumps in backward 
direction discussed in connection with eqn. (9) is negligible. We therefore 
expect deviations from the law of eqn. (17) near to and above 7"). 
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Equation (15) neglects the work done by the shear stress acting on the 
intersected dislocation during the activation process. As this would bring 
in at best a factor of 2 in v, and in many cases far less than that, we shall 
neglect this contribution to the activation energy. 

The application of the foregoing discussion to the temperature dependence 
of the flow stress of Al, Mg, Zn and Cd, which are all believed to have high 
stacking fault energies and therefore only a small separation of the partial 
dislocations, is discussed elsewhere (Seeger 1954 a, Seeger 1955b). We 
shall deal here particularly with Cu, Ag, Au, Ni and Pt which are all 
thought to have a small stacking fault energy and therefore a wide 
separation of partial dislocations. If eqn. (15) is solved for the shear 
stress as a function of strain rate and temperature we obtain 

={"6 + [H,—T (St+tkIn(NAby/e))]/v T<T a 

T= 2 ee ees 
TE f Maes hes 
T =H, /[S+h In (NAby,/€)]. < ¢5. ap See 
In 'eqns. (11) and (12) we have split the free energy of activation Gp, 
according to the equation 


with 


G,=H,—TS, ot Geter 
into the energy of activation H, and the entropy of activation S. 


Fig. 8 


To T 


Schematic representation of the temperature dependence of the flow stress S 
according to eqn. (17). 


The plot of 7 (eqn. (17)) versus 7’ is shown in fig. 8. It is characterized 
by three quantities : the high temperature flow stress t,, the temperature 
7’), and the flow stress at absolute zero tg-+7,. Here 

Teall j/ UL ee at a Lely, ee 
is the ‘ activation stress’, i.e. the shear stress that has to act on the 
dislocation line to push it through the obstacle without the aid of thermal 
energy. 
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We shall now compare the temperature dependence for the three therm- 
ally activated processes under discussion. As Cu is the only member of 
the group mentioned above for which theoretical estimates of all the 
quantities determining the temperature dependence of 7 are available, we 
shall confine our present discussion to this metal. 

Tq will be the same for all three processes if the dislocated state of the 
metal is the same. In testing the theory this has to be secured by using a 
proper experimental arrangement. /) is given by the distance [, of 
jogs in a screw dislocation. It will be of the same order of magnitude as 
Ui) —/), which are equal to the distance /’, of those dislocations in the 
forest that are effective in forming obstacles for the gliding dislocations. 
In experiments carried out after large pre-strains (Cottrell 1954, Wyatt 
1953) the forest consists mainly of extended dislocations. For Cu there 
is both experimental and theoretical evidence that the forest in an 
unstrained crystal is made up mainly of extended dislocations, too (see 
Seeger 1954). We expect therefore d‘=d"” to be about ten times as 
large as d®), as discussed above. 

For the processes (1) and (2) the free energy of activation under zero 
stress, Gp, will be taken to be equal to Hy. There may be small entropies 
of activation S™ or S®, but it is even hard to assess their signs. H, is 
made up of two contributions, one from the formation of a constriction 
(and possibly also a jog) in the intersected dislocation and the other from 
the formation of a constriction in the gliding dislocation. Both contri- 
butions depend on whether the dislocation passing through the forest is an 
edge or a screw dislocation. H,‘!) is just twice the energy of a constriction 
in an extended edge dislocation, which was calculated for Cu to be 3-9 ev 
(Schoeck and Seeger 1955).* H,'?) is the sum of the energy of a jog in an 
extended edge dislocation (4:1 ev) and the energy of a constriction in a 
screw dislocation (0:84 ev). We thus obtain 


Hy, 2P=49 ev. a 


In estimating G,®) we shall take into account both H)®) and S®. In 
§ 2 we have found H,®)=0-86 ev. S® takes into account the, thermal 
expansion of the crystal and the increased vibrational entropy of the lattice 
in the immediate environment of the vacancy. It can be found from the 
measurements of Meechan and Eggleston (1954) by using the results of 
Jongenburger (1953) and Abelés (1953) on the electrical resistance of 
vacancies in Cu. In this way we obtain S®—3-4k. (In applying these 
results to the present problem the assumption is made that S‘®) is the same 
a ee a Sn ee rr i i a A er 

* The forest is envisaged to consist of both edge and screw dislocations (and 
also intermediate types which will be neglected). Since the energy of a con- 
striction in an extended edge dislocation is 4-5 times that of an extended screw 
dislocation (Schoeck and Seeger 1955) edge dislocations form the most effective 
obstacles in the forest and determine the activation energies. It is, however, 
the screw dislocations that are responsible for the formation of jogs in the dis- 
location rings passing through the forest. 
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whether a vacancy is attached to a dislocation or not. This assumption 
is probably justified within the accuracy of the present calculation.) 

The numerical value of the logarithm in eqn. (17) and eqn. (18) varies 
somewhat with the experimental conditions under which the flow stress is 
determined. It depends directly on the strain rate and indirectly 
(through N and A) on the amount of pre-strain and on the substructure of 
the crystal. The quantities under the logarithm will be about the same in 
the three cases considered here with the exception of A, the area swept out 
by a dislocation line after each activation process. A‘ and A®) will be 
the same and larger than A“) by a factor of the order 10-10. For slow 
strain rates, of the order €-=10~4 sec~1, the logarithmic terms will give a 
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Temperature dependence of the flow stress of Cu. Only the full parts of the 
straight lines are observable experimentally. 


numerical factor of about 32 for processes (1) and (2) and of about 24 for 
process (3). These numerical values can either be deduced by choosing 
plausible values for NV, A and v), or by a comparison with experimental 
results on other metals, notably Zn and Cd (Seeger 1954 d). 
Inserting the numerical values the equation for 7',°) for Cu reads 
Byte yi 
i 3-4-1 24 
This temperature is given in the table, together with the corresponding 
temperatures for the metals treated here. The comparison with experi- 
ment will be given below. 
For Cu 7y') is about 2600°K and 7, is about 1600°K. These temper- 
atures, being higher than the melting point, cannot be observed directly. 


ev=k x 364°xK. on ae Deen ra} 
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The same will be true for the other metals with low stacking fault energies 
discussed in this paper. In a more indirect way the predictions of this 
paper can be tested by a comparison with the dependence of flow stress on. 
temperature and strain rate. 

At a given temperature the strain rate imposed in a tensile test is main- 
tained by those processes which require the lowest stress for the movement 
of the required number of dislocations. Figure 9 shows the 77’ relation 
for the three processes considered in this paper. Since the edge and the 
screw parts of a dislocation ring move independently of each other 
the shear stress that will actually be observed is given by the full line in 
fig. 9 (Seeger 1955 b). As the processes (2) and (3) are prerequisite to 
each other, the one that requires the higher shear stress will determine the 
shear stress necessary for the movement of the screw dislocations. If l, is 
of the same order as or smaller than /,’ the straight lines corresponding to 
processes (1) and (3) will intersect and the 7—7' relation will show two bends. 
The reason for this can be seen from eqns. (16) and (20): H,!Y/d’ is 
slightly smaller than H,/d’®), therefore 7, will be smaller than 7, 
provided /, is not larger than J)’. Figure 9 corresponds to 1, 1)'/2. 

Adams and Cottrell (1955; see also Cottrell 1954) found experimentally 
a temperature dependence of the flow stress of pre-strained Cu single 
crystals of the type shown in fig. 9. The temperatures of the two bends 
were 7’,=180°K and 7’,=250°K. According to the theory presented here 
the temperature dependence of the flow stress outside the temperature 
interval 7’, to 7’, is very small only and comparable with that of the shear 
modulus. It will therefore be difficult indeed to separate the direct 
temperature dependence of the flow stress due to the thermally activated 
processes (1) and (2) from the indirect temperature dependence caused by 
the temperature variations of 7,, H,') and H,). These in turn are 
caused by the temperature dependence of the elastic constants and the 
stacking fault energy. 

It is very much easier to check the prediction of thermally activated 
processes contributing to the temperature dependence of the flow stress by 
tests on the change in the strain rate. The experiments of Carreker and 
Hibbard (1953) on polycrystalline copper give indeed a strain-rate depen- 
dence of the flow stress of the right magnitude both below 200°k and above 
270°K. 

Unfortunately it is not possible to compare the calculated temperatures 
7,2) with experiment, since this temperature is not observable directly. 
Tf, however, the temperatures 7’, and 7’, are known from experiment and 
the ratio H,‘/H,\) is known from theory, it is possible to find 7).* The 
evaluation of Cottrell’s data gives 7')°)=350°K, which compares very well 
with the theoretical value of 364°x. This close agreement is a strong 


* The ratio of H) and H,®) is given to a rather high accuracy by the 
considerations of §2, although the values for these quantities themselves may 
be less accurate due to the uncertainty of the exact value of the stacking fault 

energy. The evaluation of 7’, is best done graphically as shown in fig. 9. 
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indication of the correctness of the theory presented here. Further 
evidence will be presented in § 4. 


$4. Toe Mecuanism or Piastic DEFORMATION OF NOBLE METALS 


In § 3 we have seen that a straightforward application of the theoretical 
considerations of § 2 on the generation of vacancies leads to an under- 
standing of the temperature dependence of the flow stress of copper single 
crystals. We shall now discuss briefly the implications on the mechanism 
of plastic deformation of noble metals at temperatures below the temper- 
ature of self-diffusion. We shall base our arguments on the results on Cu, 
as there is no reason to suppose that the essential features should differ 
within the noble metals group. We shall leave out Ni as we do not know 
at present whether a second bend at the temperature 7’, occurs at low 
temperatures. 

The mechanism below 7', (in copper 7';=180°K) is as described by 
Seeger (1954 b). The edge dislocations move at lower stresses than the 
screw dislocations. The form of the dislocation rings spreading across 
the glide planes shows therefore a marked elongation in the direction of 
the Burgers vector. Thermal energy is not sufficient to make thermal 
production of lattice defects an easy process. As the stress is raised 
during an elongation test the athermal generation of lattice defects as 
discussed in § 2 will become possible. In elongation tests carried out well 
below 7’, we therefore expect the defects generated by moving dislocations 
to be present mainly in the form of interstitial atoms and pairs or clusters 
of vacancies. There will be a rather small fraction of isolated vacancies 
only generated at points where /, is exceptionally large. 

Above 7', the situation isreversed. The thermal production of vacancies, 
and therefore the spreading of the screw dislocations, is easier than the 
jog-formation in edge dislocations. The edges will be held up and the 
shape of the dislocation rings will be elongated perpendicular to the 
Burgers vector. 

In this temperature range many isolated vacancies are created. If the 
temperature is low enough not to allow appreciable migration of vacancies 
they are prevented from condensing into clusters or annihilating at 
dislocation lines. The number of interstitial atoms generated during 
deformation is of the same order of magnitude as that of vacancies. The 
reason for this is that jogs creating interstitials are about as frequent as 
jogs creating vacancies. A screw dislocation stopped by its jogs will move 
forward preferably at places were vacancies will be created. It will 
bulge out and generate interstitials by the mechanisms described in § 2. 

It can be seen from eqns. (20) and (21) that the stress necessary to move 
extended edge dislocations through the forest is considerably larger than 
that necessary to move extended screw dislocations. This cutting of 
screws through the dislocation forest is the rate determining process in Cu 
deformed at room temperature. According to Cottrell (1953) the 
activation energy of such a process can be found from increment tests in 
logarithmic creep as follows. 
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In our present notation the equation for the activation energy H, at 
zero stress reads 
a ee 


ice ter ee) 
Here 7 is the applied shear stress which is assumed to be lar ge compared 
with the stress increment dr. ¢ is the time and » a constant with the 
dimension of a frequency. dInvt denotes the shift in the t-scale of the 
creep curve e=<«(¢) to which the raising of the applied stress by dr is equiva- 
lent. From a series of measurements on polycrystalline Cu at room 
temperature by Wyatt (1953—in particular table 4) it follows from 
eqn. (23) that Ho) =4-4 ev with an accuracy of about 10°, in reasonably 
good agreement with the theoretical value H,'==4-8 ey (eqn. 21).* It 
seems as if the stacking fault energy y,,=40 erg/em?, on which the latter 
value is based (Seeger and Schoeck 1953) were too low by about 10%. 

Cottrell (1954) has observed that for larger pre-strains the only effect 
on the 7-7’ curve of changes in the pre-strain is a variation of the r-scale. 
The shape of this curve is determined by the magnitude of the lengths 
L,, lp’ and J)—/) is a mean distance of the dislocations in all glide systems and 
determines 7, (Seeger 1954 b). For larger pre-strains, where slip in more 
than one glide system occurs, J) will be proportional to /,’, which measures 
the distance of dislocations threading the glide planes. As pointed out 
earlier, for constant rate of strain /; will be proportional to /,’. At 
larger pre-strains all three lengths involved are proportional to each other 
and decrease with increasing strain. This is believed to be the explanation 
of Cottrell’s observation. 

Another way of comparing the predictions of the theory with experi- 
ments is by investigating the magnitude and the annealing behaviour of 
the additional electrical resistance due to deformation (for surveys see 
Broom 1954 and van Bueren 1955). The increase of electrical resistance 
anneals out in four stages (see Appendix B), stages IT, ITI, [V corresponding 
to the disappearance of point defects created during cold work. Stage V 
constitutes the additional resistance that does not anneal out below the 
temperature of recrystallization (in Cu) and is attributed to the scattering 
of electrons by dislocations. 

It is remarkable that after deformation stage IV, corresponding to the 
migration of vacancies, is observed. This shows that Mott’s creep 
mechanism (1955) cannot be the only one capable of creating isolated 
vacancies, but that there must be at least one more mechanism. This is 
so because Mott’s mechanism requires a temperature which allows the 
vacancies to anneal out as soon as they are created. As far as the author 
is aware the mechanism employed in this paper is the only one so far 
proposed for the generation of single vacancies at low temperatures. 

Whereas stage IV is quite pronounced after deformation at room temper- 
ature it is, compared with stages IT and IIT, weak after deformation at low 


be The strain-rate in Wyatt’s increment tests was of the same order of magni- 
tude as that usually applied i in tensile test. The comparison of creep data and 
tensile test appears to be justified. 
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temperatures and increases somewhat with temperature (preliminary 
experiments by van Bueren and Jongenburger—private communication). 
This is in qualitative agreement with the present theory which predicts 
the generation of a small number of single vacancies at low temperatures. 
Anothereffect borne out by these experiments is that stage V, corresponding 
to the resistance of the dislocation lines, increases with increasing temper- 
ature. This is at first sight rather surprising, as the coefficient of work- 
hardening is practically independent of temperature in the investigated 
range of temperature and deformation (Blewitt, Coltman and Redman 
1955), so no change in the number of dislocations is expected. The 
experimental results find a natural explanation, however, in terms of the 
changes occurring in the shape of the dislocation rings with increasing 
temperature of deformation. At very low temperatures the dislocation 
rings, being oblonged in the direction of the Burgers vector, consist 
mainly of screw dislocations, and these are known to have a smaller 
electrical resistance than edge dislocations (Hunter and Nabarro 1953). 
The form of the dislocation rings changes with increasing temperature. 
consisting at room temperature mainly of edge dislocations. The 
electrical resistance due to the dislocations created during cold work 
increases therefore with temperature until recovery effects reduce the 
total number of dislocations appreciably. 

It appears from the preceding discussions that for Cu the theory agrees 
not only with the known facts about mechanical properties but also with 
the observations on the electrical resistance and its annealing behaviour. 
In the latter case, however, more complete data are required for a 
quantitative comparison between theory and observations. 
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APPENDIX A 


Evidence for the Stacking Fault Energy of Nickel and other Transition 
Metals 

As pointed out elsewhere (Seeger 1955 b) the stacking fault energy y of a 

f.c.c. or a h.c.p. metal is closely related to its electronic structure. Quali- 

tative conclusions on the magnitude of the stacking fault energy may be 


Lattice Defects by Moving Dislocations 1213 


drawn from a study of the factors governing cohesion of these metals. 
It is, however, hard to obtain quantitative results on a purely theoretical 
basis as the stacking fault energy per atom does not exceed a few per cent 
of the binding energy of the metal and is often of the order of only 10-3 
of the binding energy. It has been possible to group the close-packed 
metals into metals with high stacking fault energies and metals with low 
stacking fault energies (for details see Seeger 1955 b). The situation is 
straightforward with regard to metals without partly filled d-shells. 
Multivalent metals with f.c.c. or h.c.p. structure belong to the first group 
(Pb possibly being an exception) ; monovalent metals, in particular the 
noble metals, belong to the second group. 

In the case of transition metals the situation is more complex as both 
d- and s-electrons have to be taken into account. The author (Seeger 
1954 e) has recently undertaken a systematic study of the energetically 
favourable structures of the transition metals (‘ energetically favourable ’ 
means in this context that the cohesive energy of the structure differs by 
not more than, say, 3x10-% from the binding energy of the stable 
structure). Which of the simple metallic structures are energetically 
favourable depends on the number of d-electrons per atom. There are 
some cases in which neither the f.c.c. nor the h.c.p. structure is energetic- 
ally favourable (V, Nb, Ta, Cr, Mo, W), some others in which the h.c.p. 
but not the f.c.c. structure is energetically favourable (Ti, Zr, Hf) and yet 
others in which both the f.c.c. and the h.c.p. structures are energetically 
favourable (e.g. Fe, Co, Ni, Ru, Rh, Pd, Os, Ir, Pt). If both the close- 
packed structures are energetically favourable, the stacking fault energy is 
expected to be low ; if only one of them is favourable, the stacking fault 
energy will be high. 

The known facts on allotropy and stacking fault energies of transition 
metals are in agreement with these considerations. A further test would 
possibly be provided by a knowledge of the main glide plane of the 
hexagonal metals Ru and Os. These metals have subnormal axial 
ratios (Ru: c/a=1-5824; Os: c/a=1-5790) which are of the same order 
as that of Ti (c/a=1-5873). Tiis known not to show the basal plane as the 
main glide plane. This is not surprising, as in such a subnormal structure 
the basal plane is less densely packed than the {1010} and the {1011 }- 
planes, and (because of the high stacking fault energy of Ti) there is no 
preference for dislocations in the basal plane. In Ru and Os, however, 
complete dislocations lying in the basal plane should be able to dissociate 
freely into partial dislocations and to lower their elastic energy. This 
energy gain could perhaps be large enough to put the basal plane back in 
its role as the main glide plane of the hexagonal metals. 

Turning now to Ni this transition metal should have a low stacking 
fault energy. There is no reason why Ni, being situated in the periodic 
table between Co and Cu, both of which have definitely low stacking fault 
energies, should have a high one. It appears, however, from the experi- 
mental data so far available that the stacking fault energy yy; is somewhat 
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larger than that of Co (yg,=10 erg/em? at room temperature) and of 
Cu (y¥qy=40 erg/em?) (for evidence see below). 

The experimental evidence for a low stacking fault energy of Ni is as 
follows : 

(a) The stacking fault energy of f.c.c. metals is closely related to the 
boundary energy of coherent {111} twins and about twice as large as the 
latter. The smaller the ratio of the twin energy to the average grain 
boundary energy is the more frequent will be the occurrence of recrystal- 
lization twins. The fact that annealing twins are common in Ni and that 
their frequent occurrence makes it virtually impossible to grow single 
crystals of Ni by the strain-anneal method (P. Haasen, private communi- 
cation ; E. Kneller, private communication) is therefore an indication that 
yy; is comparatively low. The same situation with respect to the recry- 
stallization twins holds for Cu. 

(b) It was recently possible to establish experimentally the correctness 
of the suggestion by Koehler (1952) (see also Leibfried and Haasen 
(1954), Seeger (1954 f), Haasen and Leibfried (1954), Seeger (1955 b)) that 
the lamellae within the slip bands of f.c.c. metals originate from glide in 
one plane by thermally activated screw dislocations (Diehl, Mader and 
Seeger 1955). The frequency of thermally activated cross-slip depends on 
the stacking fault energy through the width of extended screw dislocations 
and also on temperature and applied stress (Schoeck and Seeger 1954). 
In Al, which has a high stacking fault energy, slip lamellae occur at room 
temperature already at stresses as low as 10-4 G (G=shear modulus). For 
both Cu and Ni, however, an applied shear stress of about 10~% G is neces- 
sary for the formation of slip bands. As the equivalent temperature of Ni 
is somewhat higher than that of Cu, this experimental result shows that 
the width of extended dislocations in Ni must be slightly smaller than that 
inCu. It appears at present that a stacking fault energy yy;=80 erg/em? 
might be a reasonable estimate. . 

Other experimental data on Ni are on (i) the temperature dependence 
of the critical shear stress, (ii) activation energy of high temperature steady 
state creep, and (iii) critical shear stress of the binary alloys Cu—Ni. 

(i) The temperature dependence of the critical shear stress has already 
been discussed in § 3. It was shown there that it fits qualitatively with 
the present picture of a low stacking fault energy. An interpretation of 
the data on the basis of a high stacking fault energy (the low temperature 
part of the dashed curve in fig. 10 being due to jog-formation in edge 
dislocations) would require that ‘extended’ dislocations are as narrow 
in Nias they are in Al, and this in turn would require that yxy; is at least 
of the order of 400 erg/em?. From the point of view of the electronic 
structure of Ni this seems to be unreasonably large. 

(ii) The activation energy of high temperature creep is found to be 
equal to that of self-diffusion. It is not possible to decide theoretically 
whether this is due to a very small or to a large splitting of complete 
dislocations into partials (Seeger 1955 a) 
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(iii) A recent investigation (unpublished) showed that the analysis of 
the concentration dependence of the critical shear stress of Cu-Ni alloys 
by Suzuki (1952) is not quite satisfactory. The resistance to flow due to 
short range order (Fisher 1954) must not be neglected. Taking into 
account this contribution to the critical shear stress of the alloys, the twin 
energy of Ni will come out smaller than given by Suzuki (1952). No 
final conclusions can be drawn as yet. 


Fig. 10 
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Temperature dependence of the critical shear stress of Ni-single crystals after 
Andrade and Henderson (1951). ©: values given by Andrade and 
Henderson ; e : extrapolated from the easy glide region according to 
Seeger (1954 a). The dotted lines show possible interpolations between 
the measurements. 


On the basis of the above evidence it is concluded that Ni has a low 
stacking fault energy of the order of 80erg/cm?. Since the electronic 
structure of Ptis very similar to that of Ni we believe Pt, the other transition 
metal treated in this paper, to have a low stacking fault energy, too. 


APPENDIX B 
Energies of Lattice Defects in the Noble Metals and in the Transition Metals 
Ni and Pt 

In this paper formation and migration energies of lattice defects have 
been used for the calculation of the temperature 7')°) and for an identifica- 
tion of lattice defects created during plastic deformation. The values 
employed for the noble metals have been given elsewhere (Seeger 1955 c). 
They are based on the classification of annealing stages observed after 
cold work, radiation damage and quenching as given by van Bueren (1955) 
and are in agreement with the energies of formation of single vacancies as 


1216 Alfred Seeger on the Generation of 


calculated by Fumi (1955). Following the Delft school, van Bueren 
(1955), and others, the annealing stage observed in Cu between 130°K 
and 200°K (so-called stage IT) has been attributed to the migration of pairs 
and groups of vacancies. Stages III and IV (activation energies in Cu 
of 0:7 ev and 1-2ev) are ascribed to the migration of interstitials and 
vacancies, respectively (following Brinkman, Dixon and Meechan 1954). 

The only experimental data on noble metals known to the author which 
do not fit into this scheme are those of Kauffman and Koehler (1955), 
which were obtained by quenching and annealing experiments on Au wires. 
These authors give (1:28-0:03) ev for the formation and (0-68-+-0-03) ey 
for the migration of vacancies in Au. The reasons for leaving out these 
data are as follows: (i) Lazarew and Ovcharenko (1955) performed 
essentially the same experiments as Kauffman and Koehler and found an 
energy of formation of vacancies in Au of 0-78 ev. This is in reasonable 
agreement with the value of 0-67 ev obtained by Meechan and Eggleston 
(1954), which is used in this paper. The activation energy of migration 
obtained by Lazarew and Ovcharenko from the annealing kinetics above 
room temperature is 0-52ev. This value, which agrees with an earlier 
value of Kauffman and Koehler (1952) of (0-40-14) ev, falls within the 
activation energies of stage II and must therefore be attributed to the 
migration of double vacancies. This interpretation is substantiated by 
the observation that migration and formation energy as found by Lazarew 
and Ovcharenko do not add up to the activation energy of self-diffusion of 
1-96 ev. (11) The calculations of Fumi (1955) suggest that the energy of 
formation of a vacancy should be smaller in Au than in Cu. For the 
latter an activation energy of 0-9 ev appears to be well established. (iii) 
It appears reasonable to assume that under non-equilibrium conditions 
a vacancy cannot carry out more than 10® jumps until it is either trapped 
at a dislocation line or combines with another vacancy to form a double 
vacancy. At temperatures where double vacancies are still very mobile 
they will anneal out even more rapidly than single vacancies. On this 
basis it can be shown that the temperatures from which the quenching 
occurred in the experiments under discussion it was not possible to retain 
the majority of vacancies. 

For the energy of formation of lattice vacancies in nickel we use the value 
1-6 ev obtained by Nicholas (1955). In quenching experiments on Pt 
wires Lazarew and Ovcharenko (1955) found 1-2 ey for the energy of 
formation of vacancies. From the annealing kinetics after quenching 
these authors found an activation energy of migration of 1-1 ev. This 
agrees well with the migration energies of 1:20 ev and 1-19 ey 
observed by Dugdale (1952) after cold-working and irradiating Pt wires. 
Unfortunately the activation energy of self-diffusion of Pt does not seem 
to have been measured. It has been pointed out elsewhere (Seeger 1955 a) 
that in high-temperature creep by the climbing-mechanism the activation 
energy Q, for steady state creep cannot be smaller than the activation 
energy @, of self-diffusion. For metals with low stacking-fault energies 
(among them Pt—see Appendix A) Q, should be equal to Q,- Dorn 
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(1954) gives for Pt Q,=2-4ev. It appears therefore that it is correct to 
attribute to both the formation and the migration of vacancies in Pt 
activation energies of 1-2 ev. 
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ABSTRACT 


In this paper elementary lattice kinetic theory is applied to the problem 
of impurity diffusion by the vacancy mechanism, special attention being 
paid to the case where there are appreciable attractions between impurity 
atoms and vacancies. We make use of the approximate methods of 
association theory, in which we treat impurity—vacancy pairs at small 
distances of separation (‘ complexes’) on a different footing from pairs 
at farther separations. Detailed formulae are presented for (i) the 
diffusion coefficient of impurities in face-centred cubic metals, for which 
the association model should be very good (Johnson mechanism), and 
(ii) the diffusion coefficient of divalent ions in polar crystals of the NaCl 
type. One important aspect of the results on divalent ions in polar 
crystals is the dependence of the diffusion coefficient D(c) on impurity 
concentration, c, when there is appreciable association. At low concen- 
trations D(c) is proportional to c but at higher concentrations D increases. 
more and more slowly finally tending to a saturation value. In the 
absence of appreciable association D(c) remains proportional to c at all 
concentrations and temperatures. The form of our results shows that 
experimental determinations of impurity diffusion coefficients are capable 
of throwing a very direct light on the conditions existing within many 
of the systems of interest in the field of colour centres (e.g. NaCl and KCl 
containing alkaline earth chlorides). 


§ 1. INTRODUCTION 
Ir is now widely appreciated that many of the interesting properties of 
crystalline solids are only understandable in terms of lattice imperfections 
or departures from perfect lattice periodicity. (For a general account of 
this subject see for example Seitz 1952.) Foreign atoms constitute one 
important type of imperfection and in insulating polar crystals, impurity 
ions with valency differing from that of the ions of the host lattice are 
of especial interest for work connected with colour centres (Seitz 1954) 
and luminescence centres (Runciman 1954). The particular interest in 
multivalent ions comes from the necessity for compensating the nef 
charge that they carry when embedded in the crystal lattice. Thus, to 
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take one of the simplest examples, in NaCl containing substitutional 
Cd?~ ions, the Cd?+ ions each carry a net charge +e and this is compen- 
sated by the creation of an equal number of Na+ ion vacancies each of 
which has a net charge —e (é@ is the magnitude of the electronic charge). 
Thus by the addition of a number of Cd?+ ions much greater than the 
number of intrinsic defects we can obtain a system containing a known 
constant number of Na* vacancies. This allows properties of the vacancy, 
such as its mobility, to be studied with far less ambiguity than is possible 
otherwise. A general complication, however, is the tendency of the Cd2+ 
ions (net charge +e) and the Na+ vacancies (net charge —e) to form 
stable pairs or ‘ complexes ’ as a result of the attractions existing between 
opposite charges. Following on earlier work by many others we have 
previously examined the effects of this association into complexes on the 
electrical and dielectric properties of these impure salts (Lidiard 1954, 
1955a). The present paper on impurity diffusion represents a continua- 
tion of the investigations into the transport properties of these systems. 
Some of our principal results have already been given in a shorter 
publication (Lidiard 1955 b). In all these calculations it is assumed 
that at all the temperatures of interest the impurity remains atomically 
dispersed and does not coagulate to form colloids. The coagulation of 
these multivalent impurities has not been much studied so far. Exceptions 
are provided by the work of Haven (1950) on MgF, in LiF, Ziickler (1949, 
quoted by Seitz (1954), p. 21) on CdCl, in KCl and Miyake and Suzuki 
(1954) on CaCl, in NaCl. 

For cases such as Cd?+ ions in NaCl, where the impurity and the 
compensating vacancy carry equal but opposite net charges we only need 
consider aggregates formed from one impurity ion and one vacancy. It 
is then a good approximation to neglect the interactions of these (neutral) 
pairs or ‘complexes’ amongst themselves or with the unassociated 
impurity ions and vacancies. In the original theory of association in 
these crystals as given by Stasiw and Teltow (1947), the mutual inter- 
actions of the unassociated impurity ions and vacancies were also ignored. 
With the help of these approximations the association—dissociation 
reaction can be described by the mass-action law. Let us adhere to the 
example of divalent ions in 1:1 salts. Let the molar concentration of 
impurity be c, let the fraction of impurity which has associated be p 
and let the concentration of thermal defects be negligible (<c). Then 
if we only consider association into pairs at the nearest neighbour 
separation the mass-action law gives us (Lidiard 1954) 


p : : 

(i—p exp (AG /RT). init Sce bs (led) 
Here z, is the number of distinct orientations of the complex ; w= 12 in 
a NaCl type lattice, where either sub-lattice is face-centred cubic, and 
296 in a CsCl type lattice where either sub-lattice is simple cubic. The 
quantity AG, is the Gibbs free energy of association, i.e. the work gained 
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under conditions of constant pressure and constant temperature in 
bringing a vacancy from a particular distant position to a particular 
nearest neighbouring position of an impurity ion. R is Boltzmann's 
constant. 

Now it is possible that the binding energy at the next nearest neighbour 
and even further separations may also be appreciable and we should then 
expect the bound vacancy to spend a significant fraction of its time at 
these remoter separations. We therefore extend the idea of association 
so as to include complexes in excited states—as we may regard these 
extended pairs. The calculated thermodynamic properties are rather 
insensitive to the exact separation at which one draws the distinction 
between associated and unassociated vacancies (see, for example, a 
calculation by Fuoss (1934) for the analogous problem in hquid electrolyte 
solutions). We therefore only need to include excited states of the complex 
if their energy above the ground state is comparable to or smaller 
than kT’. With this extended definition of association, eqn. (1.1) becomes 


) 
pt eH ex (AG,/RT) Ji, 5p AP 2) 
where z; is the number of orientations of a complex in the 7th excited 
state with energy 4G,)— 4G, above that of the ground configuration. 
The factor Xz;exp(4G,/kT) is the partition function for a single 
complex. ; 
It is assumed that an electric field will alter the relative populations 
of the various states and orientations of the (dipolar) complexes but that 
the electrolytic current is contributed solely by the unassociated vacancies. 
In a previous paper (Lidiard 1955 a) we extended this idea to alternating 
current phenomena and calculated the response of the complexes and of 
the free vacancies separately, without regard to the details of the 
dynamical equilibrium existing between them. This assumption can be 
supported in the followmg way. In the definition of association we 
include excited states until they individually no longer make an 
appreciable contribution to the partition sum in (1.2). The details of 
the dynamic equilibrium between the unassociated vacancies and the 
complexes then only concern one or two of the most highly excited states. 
Since these states do not contribute appreciably to the properties of the 
complexes, the effects upon the association—dissociation reaction of those 
agencies which upset thermodynamic equilibrium (e.g. electric fields, 
concentration gradients, etc.) will also be negligible. This is the point 
of view which has been taken in association theory up to now. In § 2 
we shall adapt this point of view to a discussion of the problem of impurity 
diffusion in NaCl type lattices, with the main emphasis on divalent ions 
diffusing through 1:1 salts. At this stage the discussion can be made 
more general than the Stasiw—Teltow (1947) model, and eqn. (1.2) can 
be replaced by an equation which allows for the Coulomb interactions 
among the unassociated ions and vacancies (see, e.g. Lidiard 1954). 
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This approach may be seen to apply only to systems containing tightly 
bound complexes, at least one of whose configurations is of such low 
energy that the mean lifetime of the complex is long compared with the 
mean times governing the response of the complexes to electric fields, 
etc. Now in the hypothetical limiting case of only random association 
the ‘ associated ’ vacancies will contribute to the electrical conductivity 
as easily as if they were ‘ unassociated’. Similarly for weakly bound 
complexes we should expect an electric field to lead to a preferential 
dissociation, and hence to a contribution to the current proportional to 
the number of complexes present. The theory which we have discussed 
so far and which is applied to the diffusion problem in § 2 does not go 
over into this limit. It needs to be extended by explicit consideration 
of the association—dissociation reaction in order to include weakly bound 
complexes. There is a practical interest in this case because the system 
NaCl+-CaCl, does not show strong association effects (Bean 1952, cf. the 
system NaCl+CdCl, as studied by Etzel and Maurer (1950) and by 
Aschner 1954). In §§ 3, 4 and 5 we shall therefore attempt to discuss 
the details of the effects which the dissociation—association kinetics have 
on the observed diffusion constants. Unfortunately we have not yet 
found it practicable to do so outside the strict framework of the 
Stasiw—Teltow association model, which assumes that there is an 
attraction between an impurity ion and a vacancy only when they are 
nearest neighbours. This is obviously only a very crude approximation, 
to the real state of affairs but it is one which gives a fairly accurate 
representation of tightly bound complexes in ionic crystals if the unknown 
parameters are chosen to give the best fit with experiment. It also 
corresponds exactly with the Johnson (1939) model for the diffusion of 
substitutional impurities in metals. Accordingly the general equations 
which we derive in § 3 will first be particularized in § 4 by dropping all 
terms connected with the net charges on the impurities and vacancies 
so as to apply to metallic diffusion. In § 5 we apply the general equations 
to the diffusion of divalent ions in 1:1 salts. In all the ionic crystals, : 
we have in mind diffusion on only one of the sub-lattices and we restrict 
ourselves to NaCl type crystals where the sub-lattices are face-centred 
cubic. Accordingly our discussion of metallic diffusion applies only to 
metals with a face-centred cubic lattice. 

The principal results for systems containing tightly bound complexes 
are eqns. (2.11) and (2.12) for the impurity diffusion coefficient. The: 
form of eqn. (2.12) is actually obvious on very general grounds. For, 
when the complexes have a long lifetime, we can think of them as 
diffusing particles and the current (in suitable units) is —D0(pe)/dx. 
Referred to the gradient of the total impurity concentration this is equal’ 
to —D,(@(pc)/dc)(dc/Ax) so that the impurity diffusion coefficient has the 
form (2.12). The main point of § 2 is thus the rigorous derivation of the 
form for Dy. This has only been done for the NaCl lattice (face-centred 
cubic sub-lattices) but the extension to other lattices is straightforward. 
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D(c) for divalent ions in 1:1 salts increases proportionally to ¢ for low 
concentrations, but for larger concentrations D(c) increases more slowly, 
finally tending to a saturation value, Dp. 

The principal results for systems containing complexes of arbitrary 
lifetime are eqn. (4.4) for face-centred cubic metals, and eqn. (5.6) for 
divalent ions in 1: 1 salts of the NaCl type. The shape of the diffusion 
isotherm (5.6) is roughly the same as (2.12) whatever the lifetime of the 
complexes. On the other hand as pointed out in § 6, the ionic conductivity 
isotherms for short-lived complexes will be appreciably different from 
those appropriate to complexes of long lifetime. The hypothesis that the 
difference between the observed conductivity isotherms of NaCl+CdCl, 
(Etzel and Maurer 1950) and NaCl+CaCl, (Bean 1952) is largely due to 
a difference in lifetime (and is not because the Ca?* ions are entirely 
unassociated) could then be tested experimentally by studying the 
diffusion of Cd?+ and Ca?* ions in these systems. 


§ 2. Dirrusion oF TIGHTLY BoUND COMPLEXES 


We begin with a discussion of impurity diffusion in systems where the 
impurity—vacancy complexes are tightly bound but where there are no 
tightly bound excited states. We consider one-dimensional diffusion 
parallel to a [100] axis in a NaCl type lattice. The NaCl lattice, being 
cubic, displays isotropic diffusion properties so that our choice of a 
[100] direction, although simplifying the calculation, does not restrict 
the results. Figure | shows the diffusion direction and the twelve nearest 
neighbour positions for a vacancy associated with the impurity atom x. 
These twelve positions can be divided as shown into three groups (a), 
(b) and (c) according to the projection of the complex on the diffusion 
direction (x-axis). The allowed orientations are shown more schematically 
in fig. 2. We shall refer the orientation of the complex to the position 
of the impurity ion. Let an(x) be the number of associated impurity 
ions per unit area of plane x ; here a is the shortest anion—cation separation 
in a NaCl type lattice. Of these associated ions let an,(a), an,(v) and 
an,(x) have their vacancies in the forward (a), middle (b) and backward (c) 
directions respectively. Then 


n(x)=n,(%)+n,(x)+ (2). 

Let w, be the probability per unit time that an associated vacancy will 
jump from one attached position to another particular attached position 
(one of four possibilities). A complex with a-orientation can jump 
backwards in this way to either of two b-orientations (and to two other 
a-orientations) and a complex with c-orientation can jump forwards to 
either of two b-orientations. Conversely every complex with a b-orientation 
can jump backwards to one of two c-orientations or forwards to one of 
two a-orientations. However, direct changes from a to c and c to a are 
also possible by an exchange of places between the impurity ion and the 
vacancy. Let this process occur with a probability w, per unit time. 
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On the basis of these definitions we see that the equations for the time 
rates of change of n,, n, and n, are, in the absence of any electrical effects, 
On,(2) } 


At = 2W4N,(X)+ Won,(X+a)—wyn,(x)—2w,n, (x), 
on x 
SN) ayo ynq(v)+ 20 nq( x) —2reyn le) —2w yn (0), nica) 
On (x 

Me =2w4n,(z)+ wen, (tx—a)—won,(x)—2wn, (x). 


Fig. 1 


X = IMPURITY ION 
O = NEAREST NEIGHBOUR 


The twelve nearest neighbours to an impurity ion (x) in a face-centred cubic 
lattice. If the impurity is part of a complex one of these twelve sites 
will be vacant. When the diffusion direction defines the 2-axis as 
shown then a, 6 and ¢ represent the planes x=a, x=0 and x=—a, 
The impurity lies at the centre of x=0. ) 


Schematic diagram showing those orientations and positions of the complexes 
which must be considered in setting up the equations of §§ 2, 3. 


As -explained in §1, when the complexes are tightly bound we 
can pursue our calculations without regard for the details of the 
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association—dissociation reaction. This does not necessarily mean that 
the presence of unassociated vacancies will not be felt in the equations 
for the motions of the complexes. Thus on the vacancy diffusion model 
the unassociated impurity ions are unable to move. Their immobility 
and the mobility of the oppositely charged unassociated vacancies then 
give rise to a Nernst diffusion potential which prevents the diffusion of 
the vacancies (see, e.g. Denbigh 1951, especially § 5.3). In fact there 
arises an electric field just strong enough to oppose the effect of the 
concentration gradient and reduce the free vacancy current to zero. We 
must therefore modify eqns. (2.1) by inserting new probabilities appro- 
priate to the case where an electric field E exists along the diffusion 
direction. It may be shown that the new probabilities are obtained 
from w, and w, by multiplying by an appropriate factor of the form 
exp (qr. E/k7’), where q is the charge on the ion which jumps (e for 
parent cations, 2e for impurity cations, etc.) and r is the vector joining 
the initial equilibrium position of this ion to its position at the saddle 
point for the. transition. (See, for example, Frohlich 1949, especially 
Chap. IIT and Wert and Zener 1949, 1950.) For fields of the magnitude 
arising experimentally we can generally replace this exponential factor 
by 1+(qr.E/RT). Abbreviating aeH/RT by A we obtain in place of 
eqns. (2.1) 


Mal) — mo} d)n,(2)+-w,(1—A)n,(a+-a) 


—w,(1+A)n,(%)—2w,(1+4A)n, (x), 


ral) ono, (1+ $)nglt) +20 (1— 3X)n (x) —4w n,(x), ies ano) 
On, (x) 


aT =2w,(1+4A)n,(x)+w,(1+A)n,(a—a) 
—w,(1—A)n,(x)— 2w,(1—4$A)n,(z), 


for diffusion in the presence of a field E along the a-axis. The value of A, 
although it will not be needed explicitly, is given by eqn. (5.3) when k,=0 
and m, is the concentration of unassociated vacancies. We shall now 
suppose that we are only interested in small concentration gradients so 
that, to the first order of small quantities, eqns. (2.2) become 


Ong 
Bp 2 ilt—n Nag) +w(n, — nq) +a, ~? — 2AMigg (Wy +2), o (2:3) 
on 
Sp = 2wrltha—2ny+N,), Seep oir ents he Ee 
ON- 


Bp elo — Me) + Wa(Ma— N.)— —aw, nt + 2dtqo(w Wy+w,). . ' (2.5) 


In the terms in \ we have neglected the differences between Na, Ny, and n, 
and their (equal) values at equilibrium, »,,. These differences are only 
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of the order of the concentration gradient, and when multiplied by A, 
which is also proportional to the concentration gradient, they give only 
second order terms—which we have agreed to drop. 

An equation for the current density j,(~) can be obtained by enumerating 
the numbers of impurity ions entering and leaving plane x and then 
counting each as contributing half an atom to the current. In this way 
we find, 


J lte)= Jaf w9(1-+ A) (ae) —wa(1—A)n,(@-+a) 
—W(1—A)n,(%)+-w(1-+A)n,(~—a)}. 
or to the first order of small quantities, 


aw. ON, Aw, On 
4 (x) =aw,(n,—n,)— —- —? —— ee ZAGW Nig. . (2.6) 


We now wish to derive an equation of the form 
j,(x)=—D,[en(x)/dx] (Fick’s law), 


so that we can deduce an expression for the diffusion coefficient Dy. 
To do this we notice that the time derivatives dn,/dt ete. occurring in 
eqns. (2.3) to (2.5) are an order of magnitude smaller than any of the 
other terms appearing in those equations. (The general diffusion equation 
shows that dn/dt is, with an appropriate conversion factor, of the same 
magnitude as 0?n/dx? which we have agreed to neglect by comparison 
with the other terms in eqns. (2.3) to (2.5).) We shall therefore set the 
time derivatives 0n,/dt etc. equal to zero. Equation (2.4) then gives 
immediately 


N(x) =3[N4(x)+n,(@)]. ei es (2.7) 
Substitution into (2.3) and (2.5) then yields, 
aw,(dn,/0x)=aw,(0n,/0x)=(Wy+We)(NM,—N,+2ANg). . (2.8) 


The left-hand equation of this pair requires that 
N,—N,=a constant, K. ene es ren 229) 


Using (2.6), (2.7), (2.8) and (2.9) it is then a simple matter to show that 
j (x)= — [a7 w.2/3(Wy+W,)|On/ox. eee (10) 


The diffusion coefficient of the complexes is therefore 
Dy=a?wW./3(Wy +g). Ee CR) 


For the continued diffusion of a complex it is necessary that each 
forward jump of the impurity ion be followed by two forward jumps 
of the vacancy about the impurity ion as centre. Whichever of these 
processes is the slower is the one which limits the value of Do. Equation! 
(2.11) shows that the mean time for a translation of the complex is the 
sum of the mean times associated with the three necessary jumps, Le. 


wo, (2w,)-1 and (2w,)?. 
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Although we are here mainly interested in the diffusion of multivalent 
impurities in NaCl type polar lattices, the assumption of a polar lattice 
has not really been made in deriving (2.11). We did indeed introduce 
the Nernst diffusion field, but \ was eliminated from the final equations 
without having been assigned its special value. The above calculation 
of D, therefore applies equally well to the diffusion of substitutional 
impurities in face-centred cubic metals when the Johnson (1939) 
mechanism is operative. 

The coefficient D, which we have calculated is not of course the diffusion 
constant for the impurity but only of the complexes. Thus the current 
in (2.10) is the total impurity current, but the concentration gradient 
refers only to the complexes. When a local equilibrium is maintained 
between the associated and the unassociated impurity atoms, we can 


D (c)/ Dp — 
po SP 
+ na oa 


id 
G 


} 0.2 0.4 perc 

ee 0.6 0.8 0 
The diffusion coefficient of divalent impurity ions in NaCl type crystals plotted 
as a function of the fractional molar impurity concentration, for three 
different temperatures. The full curves are drawn for a model in which 
the interactions are coulombic at all distances, as treated in an earlier 
paper (Lidiard 1954). The broken curves are for the same model as 
calculated with the Stasiw—Teltow approximation. The numbers on the 
curves represent absolute temperature divided by 7’), where k7', is the 

binding energy of a nearest neighbour complex. — : 


obtain the impurity diffusion constant, D(c)—such as we would measure 
experimentally—by multiplying Dy by the rate of change of » with total 
impurity concentration. We get 


D(c)=D,ld(pe)/de] Grea ope 
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where p is the degree of association. For impurities in metals p is 
independent of ¢ (see § 4) and hence D is also independent of c. For 
divalent ions in ionic crystals the dependence of D on c can be obtained 
by using the Stasiw—Teltow theory (eqn. (1.1)). Details of the predicted 
dependence are given in Lidiard (1955 b). However, it is also possible 
to calculate D(c) by using the more involved equation, connecting p 
with c, which follows when we allow for the Coulomb interactions existing 
among the unassociated ions and vacancies. Figure 3 shows diffusion 
isotherms for a simplified model in which the interactions are coulombic 
at all distances (a) as treated by Stasiw and Teltow, and (b) as treated 
previously by us (Lidiard 1954). We observe that there are no 
qualitative differences between the two cases. 

Before concluding this section we note a further possible extension by 
inclusion of excited configurations of the complexes. The appropriate 
model for this case is one where the energies of the excited configurations 
are only slightly higher than those of the ground configuration 
(AG,—4G,;< kT), but where the association energy of the ground 
configuration AG, is much greater than RT’. This case represents a 
fairly obvious extension of the simpler equations which we have derived 
above and will therefore not be considered explicitly here. 


$3. GENERAL EQUATIONS FOR STASTW—TELTOW AssocIATION MODEL 


We shall now proceed with our calculations without the restriction to 
tightly bound complexes, i.e. we shall no longer assume that w, and wy, 
are large compared with the probabilities for those vacancy jumps which 
lead to dissociation of the complex. The kinetic equations now become 
much more complicated, for we must take explicit account of the effects 
of the concentration gradient and the associated Nernst diffusion field 
on the association—dissociation reaction. As explained in the introduction 
these complications have so far restricted our considerations to a crude 
model in which all pairs at greater separations than the nearest neighbour 
distance are considered to be non-interacting. Only impurity—vacancy 
pairs at the nearest neighbour separation are considered to be associated. 

Firstly let us consider the geometry of the dissociation jumps which 
can be made by an associated vacancy. (We are considering motion on a 
face-centred cubic (sub-)lattice only.) If, starting from the position of 
the vacancy, we consider in turn all the twelve jump vectors of the 
type [110] we get twelve new possible positions for the vacancy. One 
of these corresponds to changing places with the impurity and four 
correspond to jumps into new nearest neighbour positions. This leaves 
seven jumps leading to dissociation. If the vacancy was originally in 
an a-position (figs. 1 and 2) we find that four of these seven dissociation 
jumps have increased the vacancy x-coordinate by a, two have left it 
unaltered and one has decreased it by a. If the vacancy was originally 
in a b-position the result is different ; two Jumps will increase x by a, 
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three leave it unaltered and two decrease it by a. The case of a vacancy 
in a c-position is analogous to case a; four jumps decrease « by a, two 
leave it unaltered and one increases x by a. 

Although these final dissociated positions are not all at the same 
distance from the impurity ion, we shall make the simplifying assumption 
that jumps into any one of them from an associated position are equally 
probable; we let the field free jump probability into a particular 
dissociated position be k, per unit time. We let the probability for 
the inverse jump be k,. Jumps of the vacancy which do not involve 
movement to or from an associated position will be supposed governed 
by a probability wy. In the presence of an electric field all the jump 
probabilities are altered by appropriate factors of the form 


exp (qr. E/RT) =1+ (qr. E/RT) 
(see § 2). 

We are now in a position to write down the kinetic equations for this 
system. We let m, and m, be the numbers of unassociated impurity ions 
and unassociated vacancies per unit volume respectively. In terms of 
the abbreviation A for aeH/kT we have 


On, (%)/Ot= 2w,(1—4A)n,(x)+w(1—A)n,(a+a)—w,(1+A)n,(x) 
= 910, fy) 2 ain CG er 
—k,(1+4A)n,(x)+4m,(x) {4h (1+ 4A)[m, (a+ 2a)]/N 
+ [2legn,(e-+-a)|/N + by(1—BA)frn 2), pus cea: 


On, (a)/Ot= 20, (14-3A)n,(x)+ 2w,(1—3A)n,(w)— 401 ,(x)— 2ky(1—3A)n, (2) 
—3kn,(%)—2k,(1+ dA)n,(x) + 4m,(x) {2k,(1+4A)[m,(a+a)]/N, 
+[3k m,(x)]/N,+2k,(1—4A)[m,(x—a)]/N,}, ~~. Gare 


On,(x)/ Ot=2w,(1+4A)n,(~)+w ye: r)n,(x—a) —w,(1—A)n,(x) 


—2w,(1—JA)n,(x) —4hy(1 +B vel) —Bheyn, (a We (1—$A)n,(2) 
+4m,(x){4k, ae — 2a)]/N,.+[2kym,(x—a)]/N, 
+k,(1-+4A)[m,(x)]/N,, }. ae ol ke eee 


In these equations V, is the number of cation (or anion) sites per unit 
volume ; hence m,(x)/N, is the probability of finding an unassociated 
vacancy at a particular site in lattice plane x. The last term in eqn. (3.1) 
without the factor 4 is the rate of formation of complexes with a 
particular a-orientation, by recombination of free vacancies with free 
impurity atoms. The factor 4 is inserted since there are four distinct 
orientations in group (a). Similarly for groups (6) and (c) (eqns. (3.2) 
and (3.3)). The equation for the rate of change of free impurity 
concentration is likewise, 
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Gm (x)/Bt=4le4(1—4A)nq(@) + 2eym4(@) +-hey(1-+ BA) @) + 2ky(1—4A)n, (2) 
++ Bleqn p(x) + 2k (1+ Fayre (ew) + 4h (1+-4A)n,() + 2kyn,(x) 
+-ky(1—$A)n,(x)—4m,(w) {4k9(1-+4A)[m, (a+ 2a)]/N, 
+[2kym,(%-+-a@)]/N +h (1—4A)[m,(x)]/N, 
+ Be(1+4A)fm,(e+a)]/N,+[3kym,(x)]/N, 
+2,(1—4)[on,(7—a) IN, +-44e9(1—4A)[m,(w—2a) |, 
+2k[m,(x—a)]/N + (1+4A)[m,(a)N,}. . . . (34) 


We could write down a similar equation for the rate of change of m, 
by the process of vacancy diffusion and by the dissociation and formation 
of complexes. But such an equation would have no particular value 
since additional sources and sinks of free vacancies exist in the crystal 
in the form of dislocations. These ensure that the concentration of 
vacancies is everywhere such as to satisfy some additional condition 
which we have not yet included. In metals (where of course A=0) this 
condition is that the density of free vacancies is everywhere equal to 
the equilibrium value 


m,(%)=exp (—g,/RT')=a constant 


where g,, is the Gibbs free energy of formation (see e.g. Seitz 1950). In 
polar crystals containing multivalent impurity ions the corresponding 
condition is the requirement of electrical neutrality everywhere. For 
divalent ions in 1:1 salts this means m,(x)=m,(z). The appropriate 
condition will be inserted later when we consider, particular cases. 

Before proceeding it may be noted that eqns. (3.1) to (3.4) neglect 
the existence of slight interference effects in the jumping probabilities. 
Thus the probability for a jump of a vacancy to a particular position is 
really modified by the possibility that there may be a second vacancy 
or complex already in that position. The equations for the complex 
re-orientations likewise assume that these motions are not hindered by 
the presence of other species of the mixture. These effects, being by 
assumption random, will introduce corrections to wo, w,, etc. of the order 
of the molar concentrations of the various species and are therefore 
negligible in the systems we are considering. 

Following the procedure adopted in § 2 we shall now put eqns. (3.1) to 
(3.4) into differential form by using Taylor’s expansion and retaining only 
quantities of the first order in the concentration gradient. As remarked 
previously it is consistent with this to put the time derivatives equal to 
zero (small departures only from steady state conditions). Furthermore 
in the terms proportional to A we set ,, 2», etc., equal to their equilibrium 
values (denoted by suffix 0). With these simplifications we eventually 
get from eqns. (3.1) to (3.4) the new equations 
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O= —n,(2w + Wet Thy) +2wyn,+ wn, +aw,(An,/0x) 
+ 28k,(mym,/N ,.)+40ak(m,/N ,)(m,/0x) 


+ N{1N49(— 20 — 2W2+3k,/2)4 6hemigmMy/N}, - + (3-9) 
O=2w n,—N, (40,4 Thy) + 2wyn,+ 28k(mym,/N.), . . (3-8) 


O= Won, +2wyn,—N,(2w,+W2+ 7h,)—aw,(On,/02) 


+ 28k,(mm,/N ,)—40ak,(m,/N ,)(0m,/dx) 


+A{Mqo(2W y+ 2w2—3k4/2)— 6kq(MipMy/N+)},- - - (3-7) 
O=—12k,(mm,/N,)+k (ng tn,tn,). - - - «. « «. (3.8) 


This last equation shows that the mass action constant (cf. eqn. (1.1)) 
for the association reaction remains unaltered provided that the departures 
from thermodynamic equilibrium are small. Addition of (3.5), (3.6) and 
(3.7) with use of (3.8) gives 


On,,/Ox=On,/Ox, rr ee eS): 
Le: 
N,—N,=~a constant, K. f) at eh coe enh OD 
Also, (3.6) taken with (3.8) yields 
N,=(N_+n,)/2. ere cma ee 


Before finally leaving the general equations and proceeding to a 
discussion of particular cases we must obtain the general expressions 
for the impurity ion and the vacancy current densities, j); and j,, respec- 
tively. For the discussion of impurity diffusion in the absence of 
electrostatic effects we only need j;. But when we are considering say 
the diffusion of Cd?+ ions through NaCl we need both j; and j,, since 
the diffusion potential gradient (i.e. A) is determined by the condition 
that electrical neutrality be maintained locally, ie. 


j i—Jv= 9. 


The calculation of j; is quite straightforward and the result has been 
given already (§ 2). Using (3.9) and (3.10), eqn. (2.6) becomes 


j =aw,K —a*w,(dn,/0x)+ 2dAawysn,o- Ar TEES 


On the other hand the evaluation of j,, is rather tedious and is therefore 
given in an Appendix. The result to the first order of small quantities is 


Jy=Ka(—2w,—w.+3k,)—4wya?(Om,/dx)+a?(dn,/0x)(w y—20k,) 
+ 80a7k o(m,,/N (Om 0x) — Aan 4w,+2wy+ 14k,)—4adwym,. (3.13) 


In arriving at eqn. (3.13) we have made use of eqns. (3.9) to (3.11). 
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$4. Dirrusion oF Susstrrutionan Impuritias tx Face-Cenrrep 
Cusic METALS 
In this case there is no diffusion potential and ) is zero. Furthermore 
the concentration of free vacancies is everywhere maintained at. its 
equilibrium value and dm,/dx is also zero. Equation (3.5) with (3.8), 
(3.10) and (3.11) then gives 
On, On, On, K(w,+w.t7k,/2) 


Ox Ox Ox AW, 


(4.1) 


Equation (3.12) for the impurity current then becomes 

),;=—akK (w,+7k,/2). et a ee 42) 
To find the diffusion coefficient it only remains to find the gradient of 
impurity concentration in terms of K. This can be accomplished by 
combining (3.8) and (4.1). Finally, elimination of AK from (4.2) and 
comparison with Fick’s law gives for the impurity diffusion constant 


FE LHL eee (AS) 
3 (wy+W2+ Tk,/2) 

where p=(1+k,N,/12k,m,)~1, is the degree of association in this model, 
i.e. the ratio of the number of impurity complexes to the total number 
of impurity ions. The degree of association and hence D, are both 
independent of impurity concentration; p does however depend on 
temperature. If, as seems most likely, the free energy of association 
AG=kT In (k,/k,) is less than g,, the free energy of formation of an 
isolated vacancy, then as 7-0, p+0 as 12 exp (4G—g)/kT. We shall 
discuss two special cases of (4.3) : 


(i) Tightly Bound Complexes 
Here k,<wy,, Wy, i.e. the complexes have a lifetime [(74,)~1] which is 


long compared with the mean times associated with the w, and w, jumps. 
For this case the diffusion constant becomes 


i.e. the result (2.11) essentially. 


(ii) Loosely Bound Complexes 

_ There are several possibilities here but it is instructive to discuss the 
limit corresponding to an impurity of the same chemical species, i.e. a 
radio-isotope as ordinarily used in self-diffusion measurements. ww, Wp, 
k, and k, are all equal to wy, the free vacancy jumping probability, so 

that for m,/N,<1, eqn. (4.3) gives 

9 (4a7?wym, 

= — | —— ek tera eee (oS 
oe a rea et 
i.e. nine-elevenths of the self-diffusion coefficient as usually calculated for 
tracers. However, as pointed out by Bardeen and Herring (1952), the 
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simple calculation of the tracer diffusion coefficient neglects the fact that 
successive jumps of a given tracer atom are correlated with one another. 
Thus if a tracer atom has just jumped into a vacancy tts next jump is 
not entirely random since it is possible that it will move back into the 
position it had just left. The factor m,/N, which expresses the assumption 
of complete randomness in the tracer’s jumps must therefore be reduced 
somewhat. By a rigorous random-walk calculation one can derive* the 
factor 0-80 for a face-centred cubic lattice, compared with our factor 0-82. 
The slight difference arises because the present calculation is also 
incomplete in its handling of the correlation between jumps. We have 
dealt with nearest neighbour correlations by our use of a rudimentary 
pair distribution function (numbers of complexes) but we have not taken 
account of the statistical correlations as they concern pairs at second 
and farther neighbour separations. Nevertheless it is gratifying that the 
present calculation is so accurate when applied to the problem of tracer 
diffusion since the method is designed for systems where the correlation 
effects are largely due to an attraction between the impurity and the 
vacancy. 


$5. Dirrusion or DivaLeNntT ImMpurtiry Ions IN 4 1:1 Saur 
witH NaCl STRUCTURE 
In this case electrostatic effects play a large part in determining the 
nature of the general solution. Thus we shall solve eqns. (3.5) ete. 
by using the condition of electroneutrality. Since the vacancy and the 
divalent ion carry equal and opposite net charges this condition requires 
m;=m, everywhere and also j,=),. Use of this last equation will enable 
us to fix the value of A. 
The relevant equations are obtained as follows. Firstly eqns. (3.5) and 
(3.8) (or alternatively (3.7) and (3.8)) yield 
K(wy+wet 7hy/2)=a(wyt5ky)(On,/0x)+ANgo(—2w,—2w,+3k,), . (5.1) 
where A is defined by (3.10). Also, using (3.8) to simplify (3.13) for 7, 
and then setting j,—=)j; gives 
A (2w,+ 2w,—3k,)=a(2w,—10k,)(0n,/0x)—4wya(dm,/0x) 
—A{Nqo(4w,+4w.+14k,)+4womy}. . (5.2) 
The third equation which we need is eqn. (3.12) for j;, namely 
Ji = dW, K — aw, (0Ng/0x)+ 2AAWgNao- eee Sal 2} 
In principle we can now eliminate K and X from these three equations 
and express the impurity current in terms of dn,/dx and dm,/dx; both 


* Bardeen and Herring (1952) give the value 0-90 for this factor. However, 
as pointed out by LeClaire (private communication to the author), eqn. (A.2) 
of their paper is incorrect. This should be replaced by 

f=1+20,P,+265P,?+. . 
When this correction is made the correct result 0-80 ‘is obtained. This has 
been realized also by Haven (1955). 
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of which may be related to the gradient of impurity concentration by (3.8). 
Elimination of A from (5.1) and (5.2) gives an equation for A, i.e. for the 
diffusion potential 


A= — ®i5ak,(0n,/0x)(2w,+w,+2k,) 


+ aw (Om,/0x)(2w,+2w.+7k,)}, . . (5.3) 
where 


D= (20k Ny9(W 1+ w.+h,)+ wom, (2w,+2w,+7k,)}-. (5.4) 
We can now also find K and hence the current density j,, 
Ji=— {20k 2 49(Wy +1) + wom, (2w,—3k,) }a2w,G(dn,/dx) 
— 20k yn, 9Wpt?w,(Om,,/0x). eae ee eect ie (525) 
We can now make the following substitutions and get the diffusion 


constant D(c) as a function of the degree of association p, i.e. implicitly 
as a function of c. We have 


Np =F . pe, 
m,—WN_.c(1—p), 
One TSO 1 a(N_.c) d(pc) 
= - — {4 ele == = 
Ox 3 Ox (Ng tN, I N.) 2 Ox de ’ 
om, <a(N.c) a 
Also, using the fact that on the Stasiw—Teltow association theory (eqn. 
(3.8) 


d(pe) 2p 


de  — (1+p)’ 
we arrive at the following expression for the diffusion constant 
aw, { 2p Xp(W3) 
D(oj= *( ) oe eee tO. 0 
3 Ip Xp(Wit Ws) ( 


where 
Xp(W1)= Wo 2+ 7Th,) — p {2wyw,— k(20w,/3—Twy)— 20k ,"/3 }. 
We may note several features of this result. 
(1) In the second part of the expression for D(c) the denominator is 
the same function of (w,+-w,) as the numerator is of w,. Hence if w,<w, 
then. irrespective of the value of /,, 


_ @w, 2p \_ aw, a(pe) aie 
D{(6\= 3 (So oe Ne Sia) 


as is clear on elementary grounds (Lidiard 1955 b). Result (5.7) shows 
that the simple form for the diffusion constant appropriate to tightly 
bound complexes (§ 2) will also hold for complexes with any lifetime 
provided only that w,<w,. Expressed in words this condition states 
that the mean time for a direct impurity—vacancy exchange should be 
long compared to the mean time for jumps of the vacancy about the 
impurity ion as centre. We might expect this to be the case for impurity 
ions of large radius. 
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(2) Under conditions of low association energy or high temperatures 
or small impurity concentrations, the degree of association will be much 
less than unity. Equation (5.6) then gives 

D(c)= a0 (w+ Thy/2) 2p, 
3(wW + W2+ 7k,/2) 
which may be compared with (4.4). The additional factor 2 inj the 
present case comes from the fact that the degree of association in ionic 
crystals is proportional to ¢ at low concentrations whereas in metals, it 
is independent of c; hence d(pc)/de in metals is equal to p but in ionic 
crystals it is 2p. 


(5.8) 


§ 6. CONCLUSION 

In the first part of the present paper (§ 2) we have considered diffusion 
in systems where the complexes are tightly bound, i.e. systems for which 
we could proceed without explicit consideration of the association— 
dissociation reaction. We only studied the case where the ground state 
of the complex was of much lower energy than any of the excited states, 
but we indicated how this limitation could be removed. From previously 
published work (Etzel and Maurer 1950 and Teltow 1949) it seems likely 
that the calculations for tightly bound complexes are applicable to 
systems such as NaCl+CdCl, and AgBr+CdBr,. If this is the case 
then as pointed out previously (Lidiard 1955 b) measurements of the 
diffusion of Cd?+ ions in these salts are capable of providing information 
about the degree of association with a directness which is not possible 
with ionic conductivity measurements. It should perhaps be pointed out 
however, that direct measurements of the extent of the association are 
now also possible by application of nuclear resonance techniques (Cohen 
and Reif 1955). 

While the system NaCl+CdCl, seems to afford a good example of 
tightly bound complexes, the ionic conductivity measurements made by 
Bean (1952) on NaCl+CaCl, indicate the absence of any appreciable 
association in this system. Partly for this reason we have therefore 
made an examination of the effects of short lifetime on the predicted 
properties of the simplified model described in §3. In general the 
predicted diffusion isotherm (eqn. (5.6)) is more complicated than is the 
case for tightly bound complexes. However in the special case where 
the jump time for the associated impurity ion (w,) is long compared with 
the jump time for motion of the vacancy from one associated position 
to another (w,)—but without restriction on the lifetime of the complexes 
—then the diffusion isotherm again takes on a simple form (5.7), easy 
to relate to experiment. 

In this paper we have not explicitly considered effects connected with 
external applied fields. However it is not difficult to see from eqns. (3.5), 
(3.8), (3.12) and (3.13) that the d.c. electrical conductivity is 

4a%e?w, 40a7e7k,(wy+w,t+k,) 


psoas RM See OP ee eA ee 
o= pp SUP) Ser, tobik ey tee oe Bal 
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Only the first term of (6.1) is given by the Stasiw—Teltow theory (see 
Lidiard 1954) which essentially applies to systems where the complexes 
have long lifetimes (small £,). When studying NaCl-+CaCl,, Bean (1952) 
inferred the absence of association from the linearity of the o vs. c isotherms. 
The idea here is that the (1—p) factor introduces a curvature into the 
isotherms. However, we now see that if k, is appreciable the presence 
of the second term in (6.1) may easily be sufficient to bring this curvature 
back to zero without the necessity for negligible association. Measure- 
ments of the diffusion of Ca?+ ions through NaCl-+-CaCl, should enable 
us to determine whether /, is indeed large for this system. Negligible 
association, such as Bean inferred, would give a diffusion constant 
closely proportional to the concentration of Ca?+ ions. On the other 
hand the effects of association on the diffusion isotherms should be 
clearly visible whatever the value of k, (eqn. (5.6) or possibly (5.7)). 
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A:P RB EONt DAR: 


In this appendix we shall obtain the full expression for j,, the total 
vacancy current density. The expression which we obtain may then 
be simplified to apply to conditions where the departures from equilibrium 
are small, and will then be found to give the formula (3.13) already 
employed. 

We proceed in the same way as in the evaluation of j; in § 2. That is, 
we shall enumerate the numbers of vacancies entering and leaving plane x 
and count each as contributing half an atom to the current. If the 
vacancy leaves in the forward direction it contributes +} to j,, while if 
it leaves in the backward direction it contributes —} to j,. Similarly 
for vacancies which arrive at plane x. Movements of vacancies within 
plane a may, of course, be ignored. 

In the following we need the numbers of complexes, ions etc., in 
particular lattice planes. From our definitions these are a times the 
corresponding quantities n,, %,, %., m,; ete. There is therefore a factor a 
in j,. There is also a factor } since each jump in or out of plane x 
contributes either +4 or —} to the current. We therefore remove the 
factor 4a at the beginning and assert that 7, is $a multiplied by the sum 
of the terms which we now enumerate. 

Firstly consider the terms which come from the re-orientation jumps 
of the complexes (w, and w, jumps). From fig. 2 it is clear that they 
contribute to 2,/a an amount 

ie ue ee ere 
ae ra ae a) 2w (1 + 3A)n,(%—a) 
—2w,(1+4A)n,(a@+a)+ 2w,(1—4A)n,(x—a) 
4-w,(1—d)n,(a-+a)-+2w(1—}A)n,(«+4). 

Next consider the terms coming from dissociation of complexes (k,). 
We recall that of the seven ways in which an a-complex can dissociate, 
four lead to an increase of a in the vacancy’s 2-coordinate, two leave it 
unaltered and one decreases x by a. For a b-complex two of the 
dissociation jumps increase the vacancy x by a, three leave it unaltered 
and two decrease it by a. Lastly for a c-complex only one dissociation 
jump increases the vacancy a-coordinate by a whereas two leave it 
unaltered and four decrease it by a. We then see that dissociation jumps 
contribute to 2),/a an amount 

—k,(1+4A)n,(v)+4hk,(1—4A)n,(a—a) 
—k,(1+4A)n,(a—a)+4k,(1—4A)n,(a—2a) 
— 2k (1+ 3A)n,(@+-a) + 2k (1—3A)n, (a) 
2k,(1+-4A)n,(a)+ 2k,(1—4A)n,(a—a) 

— i 1+-3A)n,(v-+ 2a)—4k,(1+4A)n,(a+a) 
+hy(1—}A)n,(a+a)+-hy(1— 3am, (2). 
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Thirdly we find the contribution to 2j,/a coming from jumps of free 
vacancies into associated positions. Free vacancies at 2 can contribute 
both by jumping forwards and backwards. We get 


M,(%) X 164,(1—$A)m (a+ 2a)/N,. (c at v+2a) 

M(x) xX 8k(1—4A)m,(x-+a)/N (6 at x+a) 

mM,(%)x 4h(1—4A)m,(x)/N, (a at x) 
mlx) x Akeg(1+-4A)m,(x)/N , (c at x) 
—™M,(%)X 8k(1+4A)m(a—a)/N, (6 at x—a) 
—M(%) X 16k,(1-+-4A)m (x —2a)/N (a at x—2a). 


Free vacancies at (x—a) can only contribute through jumping forwards, 
giving 


M,(x—a) X 16k,(1—4A)m (x+a)/N, (c at x+a) 
M(%#—a)X 8k(1—sA)m,(x)/N., (b at x) 
Mm,(x—a)xX 4k(1—dA)m(x—a)/N, (a at x—a). 


In the same way free vacancies at (~-+-a) can only contribute by jumping 
backwards. They give contributions to 2),/a of 


—m,(x+a) x 16k,(1+4A)m (~—a)/N (a at x—a) 
—m,(x+a)x 8k,(1+4A)m,(a)/N, (b at x) 
—m,(«+a)x 4h(1+4A)m,(a+a)/N, (c at x+a). 


The last terms to be enumerated are those which come from the motion 
of free vacancies from one unassociated position to another. They are 


4A)m,(x)—4w,(1+ 3A)m,(x) 
4w,(1—4A)m,(x—a) —4w (1+ A)m,(x+a). 


The sum of all these terms which we have written down is 2j,/a. In 
asserting this we are omitting any interference terms. Thus the number 
of free vacancies jumping to new unassociated positions should strictly 
take account of the fact that not all final positions are equally likely, 
owing to the presence of the other ions and vacancies. However, this 
fact will introduce a correction to w) which is only of the order of the 
impurity concentration and such corrections may therefore be neglected. 
To the first order of small quantities, 7, may then be shown to reduce to 
eqn. (3.13). 
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ABSTRACT 


The 27-day recurrence tendency of the cosmic ray intensity has been 
studied for the nucleon, meson and electron components of the cosmic 
radiation and the amplitude of the variation is found to be greater for the 
nucleon component, furthermore the changes in cosmic ray intensity 
seem to precede those of magnetic activity by 4 or 5 days. The 27-day 
period corresponds to the time of rotation of the sun and this would seem 
to indicate that the cosmic ray intensity changes are associated with solar 
disturbances which, after a delay of several days, produce magnetic 
storms. There is some evidence to show that the recurrence tendency 
is due to decreases in intensity which occur at 27-day intervals. 


$1. InrTRopUCTION 


Ir is now well established that the cosmic ray intensity shows a marked 
27-day recurrence tendency, similar to that shown by geomagnetic storms 
of moderate intensity. The continuous observation of intensity varia- 
tions using ionization chambers and counter telescopes is limited to 
primary particles with energies above about 10!eyv. It is also possible 
to study the time variations in a lower energy region. The cosmic ray 
neutron flux at sea-level is largely produced by latitude sensitive primaries, 
ie. by primary particles of energy between 10° and 10!ey. Simpson 
et al. (1953) have shown that the intensity variations in this energy region 
can be studied by means of neutron monitors. 

From June 1952 to June 1954 continuous measurements of the cosmic 


ray intensity have been made at Manchester using a neutron pile monitor, 


and two arrays of Geiger—Miiller counters inclined at 45° to the vertical 
in the North and South directions. It has therefore been possible to 


' compare directly the intensity variations found in the nucleon component 


of the cosmic radiation with those found in the ionizing component. In 
this paper the experimental results on the 27-day recurrence tendency 
in cosmic ray intensity are reported. 


* On leave from National Physical Laboratory, Pretoria, South Africa. 
+ Now at University of Ceylon, Ceylon. 
t{Communicated by Professor P. M. S. Blackett, F.R.S. 
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§ 2. THe Neutron DerectTor 


The details and dimensions of the neutron pile monitor are shown in 
fig. 1. The entire pile was housed in a thermally insulated hut, and the 
temperature was kept constant by thermostatically controlled heaters. 
The counters used were of the copper-in-glass cathode type with a central 
anode wire of tungsten. The effective length of each counter was 50 cm 
and the diameter 4 cm. They were filled with BF, gas enriched with the 
B-10 isotope at 40cm Hg pressure, and operated in the proportional 
region. Three counters were connected in parallel and directly coupled 
to a high-gain linear amplifier. Output pulses from the amplifier were 
selected by a discriminator circuit which also provided the correct pulse 
shape for the scaling circuit. The pulses were scaled by 16 before being 
fed to a drive unit which operated electro-mechanical registers. These 
registers, together with a clock, an aneroid barometer and a mains 
voltmeter were mounted on a panel, which was automatically photo- 
graphed every 15 minutes. 


Fig. 1 
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The apparatus was designed to contain two trays of counters operating 
colaterally, but the experimental results described in this paper were 
obtained with only one set of counters installed. The counting rate was 
about 2000 per hour. 

» To check its overall performance, the neutron monitor was calibrated 
periodically by placing a (Ra, Be) source at a fixed position near the 
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centre of the pile. By subtracting the number of neutrons generated in 
the pile by the cosmic radiation from the number recorded with the 
source in position, the instrumental efficiency for the detection of neutrons 
could be determined. From June 1952 to March 1953 the efficiency 
gradually increased by about 4% of its original value. From April 
1953 to June 1954 the change has been less rapid. 


§ 3. EXPERIMENTAL RESULTS 


(a) The Influence of Meteorological Variables on the Neutron Intensity 
Measured at Sea-level 


If the number of neutrons locally produced in a neutron pile, were 
related to the primary radiation only by secondary nucleons, the neutron 
intensity would depend only on the atmospheric pressure, because of the 
well-known mass absorption effect. The propagation of the nucleon 
component through the atmosphere is, however, more complicated than 
this, since secondary mesons can give rise to further high energy nucleon 
production in the atmosphere, or star production within the pile. The 
presence of these unstable links could therefore result in a dependence of 
the intensity of the nucleon component at sea-level on atmospheric 
temperature as well as pressure. 


Table | 
Group 800 mb 500 mb 200 mb 100 mb 
1 — 0-08 —0-01 —0-02 —0-10 
2 — 0-21 0-0 +0-24 + 0-05 
3 — 0:26 —0-29 —0-29 —0-23 
4 + 0-07 +0-01 +0-23 —0-10 
Mean —0-12 —0-07 -+0-04 —0-10 


To determine whether the neutrons had any appreciable temperature 
effect, a partial correlation analysis was carried out for a total of 118 
days during the period June 1952 to October 1952. The variables used 
were the daily neutron intensity N, the mean atmospheric pressure B, 
and the mean height H of different pressure levels in the atmosphere. The 
upper air data were taken from the radio-sonde measurements of the 
Liverpool Meteorological Station. The partial correlation coefficients 
wu, Obtained from the analysis are shown for four monthly groups of 
data in table 1. 

For all the pressure levels considered ryq , had a very small value, and 
never reached the 0-1 level of significance. It was therefore assumed that 
the neutron intensity depends only on the atmospheric pressure, and 
the daily neutron intensity was correlated with the mean atmospheric 
pressure for twelve monthly periods from June 1952 to June 1953. The 
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mean value of the correlation coefficient was — 0-97, and the weighted mean 
of the barometric coefficient was 
Byp=—10-4+0-1% per cm Hg. 
The error given is the standard deviation calculated from the scatter of 
the individual values about the mean. The absorption mean free path 
in air for the nucleon component as calculated from the barometric 
coefficient is 
L=130+1 ¢ cm-, 


(b) The Relative Amplitudes of the Recurrence Tendency Measured by 
Different Detectors 


Using barometric coefficients of —10-4°%, per em Hg for the daily 
neutron intensity, and —2-64°% per cm Hg for the ionizing intensity, 
the data obtained at Manchester during the period June 1952 to June 
1954 were corrected for fluctuations in atmospheric pressure. The daily 
intensity totals were then plotted as percentage deviations from the 
mean intensity. As an example of the results obtained, fig. 2 shows the 
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Day-to-day changes in cosmic ray intensity from December 4, 1953 to 
March 18, 1954. 


day-to-day changes in the cosmic ray intensity for the period December 4, 
1953 to March 18, 1954. Local time midday-to-midday totals were 
used. The variations in the ionizing intensity, for which the mean of the 
two directional telescopes were used, were multiplied by a factor of 2-3. 
The reason for this will be discussed later. The standard deviations of the 
daily intensity totals due to random fluctuations were ~0-5%. for the 
neutron and ~0-1% for the ionizing intensity. 

Because of the rather large statistical uncertainty of the daily neutron 
intensity, it is not possible to be definite about the day on which a maximum 
or a minimum in the neutron intensity occurred. Nevertheless, it is 


obvious that : 
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(1) Although the neutron detector responds to low energy primaries, 
and the ionization detector to relatively higher energy primaries, there is 
a fairly good agreement between the intensity changes measured by the 
two types of detectors. Perfect agreement cannot be expected since no 
temperature correction has been applied to the data for the ionizing 
component. 

(2) The intensity records for the period under consideration show a 
remarkable recurrence of maxima and minima with each series having a 
period of approximately 27-28 days. For example, the maximum 
occurring on December 9, 1953 is followed by other maxima near January 7, 
February 1 and February 28, 1954. 

(3) Occasionally more than one series of 27-day variations may be 
present. Thus the smaller maximum occurring near December 27, 1953 
is followed by fairly prominent maxima near January 23, February 18 
and March 15, 1954. 

Fig. 3 
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Recurrence diagrams for K-figures, neutron intensity and ionization intensity 
for period June 20, 1952 to July 20, 1953. 


This recurrence tendency could be seen in all the cosmic ray intensity 
records, and a rigorous analysis was therefore carried out to determine the 
period and amplitude. The superposed-epoch method developed by 
Chree (1913) was used in this analysis. Days of maximum and minimum 
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intensity were selected at the rate of about one per month. Plotting the 
average intensity as a function of the day-number, showed a recurrence 
tendency in the maxima and minima of both the neutron and ionizing 
intensities. For the purpose of smoothing, the recurrence diagram 


obtained by using minima as zero days, was subtracted from that obtained 
by using maxima as zero days. 


Fig. 4 
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Recurrence diagrams for neutron and ionizing intensities. eeageyee Ons 
from the mean of the ionizing intensity have been multiplied by 2:3. 
Data used covered the period June 20, 1952 to July 20, 1953. 


The difference curves obtained from a Chree analysis of the daily neutron 
and ionizing intensities for the period June 20, 1952 to July 20, ag are 
shown plotted in fig. 3. The difference curve obtained from a Chree 
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analysis of the K-magnetic figures at Abinger is also shown plotted in 
fig. 3. The same zero days have been used for all these diagrams. It 
is obvious that both the neutron and ionizing components exhibit a 
27-28 day recurrence tendency. 

To determine the relative amplitudes of the 27-day recurrence tendency 
for the neutron and the ionizing components the percentage deviation 
from the mean of the ionizing component for each day-number was 
plotted against the corresponding data for the neutrons. These points, 
for which the mean deviations over three consecutive days were used, 
are shown plotted in fig. 4. The correlation coefficient between the two 
sets of data is 0-82. The slopes of the two lines drawn in fig. 4 represent 
the limits of the regression coefficient. These limits, which result if 
infinite weight is given to one or the other of the two co-ordinates, are 
1:83 and 2-73. Equal weights were assigned to the two co-ordinates, 
and a value 2-30 obtained for the regression coefficient by the method of 
least squares. To check this the Chree diagrams, upon which the above 
analysis was based, have been plotted again in fig. 5 with the percentage 
deviations from the mean of the ionizing intensity multiplied by a factor 
of 2-3. There is a very good agreement between the sets of results. 
The 27-day neutron variations are therefore about twice as large as the 
variations in the ionizing component. 

In addition some earlier data obtained by Elhot and Dolbear (1950) on 
the hard component, and by Dawton (1952) on the hard and soft compo- 
nents have been analysed using the superposed epoch method. Elliot 
and Dolbear made their measurements with two directional counter tele- 
scopes which were shielded by 35 cm of lead. Dawton’s measurements 
were made by using two similar counter arrays, one containing no 
absorber and the other containing 10 cm of lead. The unshielded array 
measured the total flux, the shielded array the intensity of the hard com- 
ponent, and the difference in the rates gave a measure of the soft com- 
ponent intensity. Figure 6 (a) shows the Chree diagrams obtained from 
Elliot and Dolbear’s measurements on the hard component. The inten- 
sity records obtained between January 5, 1949 and May 3, 1950 were 
used in this analysis. There is a good similarity in the points obtained 
for the North and South directional records, the correlation coefficient 
being 0-82. The approximately 27-day recurrence tendency is quite ap- 
parent, and seems to persist on the average for at least four solar rotations, 
with an amplitude of about 0-2%,. The points plotted in fig. 6 (b) were 
obtained from Chree analyses of Dawton’s measurements on the hard and 
soft components. The diagram drawn in this figure is the Chree diagram 
for the total cosmic ray intensity. The above analyses covered the period 
May 9, 1950 to April 30, 1951. There are two points worthy of note : 


(1) There is a good agreement between the total 7’, hard H, and soft S 
intensities, and the correlation coefficients were found to be 


rp, 4=0-97, 1p g=0-91, ry s=0-85. 
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These values all reached the 0-01 level of significance. The limits of the 
regression coefficients are 


‘Fig. 6 
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Recurrence diagrams for cosmic ray intensity in the North and South directions. 
The measurements have been made using telescopes shielded by 35 cm 
lead. Data used covered period from January 5, 1949 to May 3, 1950. 
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Recurrence diagram for the total cosmic ray flux. The points shown plotted 
were obtained from Chree analyses of the hard and soft component 
intensities. Data was from May 9, 1950 to April 30, 1951. 


With equal weights assigned to the different sets of data, the method of 
least squares gave the following values for the regression coefficients : 


Br a~l 9, Br, g7~w9-87, Bu, sw 0-87. 


In view of the large statistical errors involved in the measurements of the 
soft component intensity, however, it is doubtful whether any significance 
can be attached to the departures from unity of the regression coefficients, 
and it therefore seems very likely that the 27-day recurrence tendency 
has the same amplitude for the total, hard and soft intensities. 

(2) There are two distinct series of peaks present. The first series per- 
sists for three solar rotations with peaks of diminishing amplitude at 28, 


1246 I. J. van Heerden and T. Thambyahpillai on the 


56 and 84 days. The second series has a small peak at 46 days, followed 
by three definite peaks at 73, 100 and 127 days (corresponding to the 
dotted lines). 

From the above it is obvious that the neutron intensity as well as the 
hard and soft component intensities show a 27-day recurrence tendency. 
and that the ampitude thereof is the same for the total cosmic ray flux, 
the hard and soft components, and approximately one-half that of the 
neutron component. 


(c) The 27-day Recurrence Tendency and Geomagnetic Field Variations 


Moderate geomagnetic disturbances are known to possess a pronounced 
27-day recurrence tendency. This recurrence tendency is clearly shown 
in the Chree diagram for the A-figures shown in fig. 3. It is also obvious 
from fig. 3 that there is a strong negative correlation between the Chree 
diagram for the A-figures, and those for the neutron and ionizing inten- 
sities. It is known that decreases in intensity of up to 8°, can occur at 
the time of intense geomagnetic storms and the relationship shown in 
fig. 3 would seem to be of the same kind, viz. an intensity decrease at 
the time of a magnetic disturbance. However, there appears to be 
a difference in detail. Inspection of the A-figure and cosmic ray intensity 
plots in fig. 3 shows that the two phenomena are not quite inversely 
related, i.e. the time of minimum cosmic ray intensity is not quite the 
same as for maximum /-figure. 

In order to determine more precisely the relationship between the 27-day 
recurrence tendency of the cosmic ray intensity and that for the geo- 
magnetic activity, the following procedure was used. The correlation 
coefficients between the neutron and ionizing intensities and A-figure 
diagrams shown in fig. 3 were evaluated, as the magnetic data was moved 
by an integral number of days in the positive or negative direction relative 
to the cosmic ray data. For instance, for a separation interval of —1 day, 
the neutron intensity on the zero day, Do, i.e. No, was correlated with the 
magnetic activity on day D_,,, ie. K,,, N,, with K,, and soon. Fora 
separation interval of ++ 1 day the neutron intensity on day D,,,i.e. V_, was 
correlated with the magnetic activity on day Do, i.e. Kg, N,. with K ,, and 
so on. Figure 7 shows a plot of the correlation coefficients obtained for 
different separation intervals in the case of the neutron intensity, whereas 
fig. 8 shows the corresponding data for the ionizing intensity. Maximum 
correlation occurs in both cases for a separation interval of about —5 days, 
indicating that the changes in cosmic ray intensity precede those of 
magnetic activity by about 5 days. 

Forbush (1954) has also studied the relation between cosmic ray de- 
creases and various measures of geomagnetic activity, and has shown that 
these intensity decreases occur during magnetic disturbances. On the 
other hand, Simpson (1954) and Kane (1955) have put forward some evi- 
dence indicating that the maxima of the cosmic ray intensity precede the 
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Fig. 7 
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Correlation coefficients between neutron intensity and K-figures for different 
separation intervals. 


Fig. 8 
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Correlation coefficients between ionizing intensity and K-figures for different 
separation intervals. 
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K, maxima by a few days, and they conclude that the most probable 
time for the occurrence of a geomagnetic disturbance is near the beginning 
of a cosmic ray intensity decrease. The high degree of negative corre- 
lation between the cosmic ray intensity and the geomagnetic activity in 
the above results does not necessarily mean that the cosmic ray intensity 
changes are causally related to the geomagnetic field disturbances. The 
two effects may have a common cause. Simpson (1954) has shown that 
the day-to-day variations in the neutron intensity at latitudes above the 
knee of the intensity vs. latitude curve were very similar to those observed 
at low latitudes, and he therefore concluded that the 27-day recurrence 
tendency of the neutron intensity and of the highly correlated ionizing 
intensity is not produced by geomagnetic field variations. The 27-day 
recurrence corresponds in time to the rotation of the solar equatorial lati- 
tudes. The present results therefore suggest that the cosmic ray intensity 
changes are associated with solar disturbances, which, after a delay of 
several days, produce magnetic storms. 

From the point of view of interpretation, it is important to know 
whether the 27-day variation is due to (a) screening off of part of the 
primary radiation, (b) the production of additional cosmic ray particles, 
or (c) a modulation of the primary radiation. 

Suppose the effect is due to the screening off of part of the primary 
‘vadiation. In fig. 9 the normal intensity, assuming an isotropic distri- 
bution, would be represented by /,; the radiation screened off is repre- 
sented by the shaded area. Under these circumstances, the mean inten- 
sity over a specific 27-day period, J,,, would be higher when the peak-to- 
peak amplitude was smaller, and vice versa. If on the other hand the 
effect is due to the production of additional cosmic ray particles, 7, may 
be taken as the normal level for isotropic distribution, and the mean inten- 
sity, /,,, would be lower when the peak-to-peak amplitude was smaller. 
and vice versa. In the case of a modulation effect, /,, would represent 
the normal intensity, and remain constant and quite independent of the 
peak-to-peak amplitudes. 


Schematic representation of the 27-day variation. 


The neutron intensity has a negligible atmospheric temperature effect 
and is therefore very suitable for investigating these three possibilities. 
Unfortunately, the day-to-day changes in the neutron intensity were to 
some extent obscured by the progressive increase in counting rate with 
time, due to changes in, the efficiency of the detecting system. A further 
correction for this gradual increase in efficiency of the neutron detector 
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was therefore applied to the barometer-corrected daily intensity totals. 
A procedure was then adopted which gives a measure of the peak-to-peak 
amplitudes of separate 28- -day periods. This measure, D, is the difference 
in the mean intensity of a 7-day interval centred around a day of maxi- 
mum intensity and the mean intensity of a similar interval centred around 
the day of minimum intensity following the maximum. The procedure 
for computing D is shown in fig. 10 (a). The next step was to calculate 
the mean intensity of each separate 28-day period, each period being 
taken as from the 7th day before the maximum. The results are given in 
table 2. 


Table 2 
Mean of 5 Mean of Mean of 
No. on ae of 7-day Date ol 7-day | D-value| 28-day 
graph maximum obs minimum : : : a 
interval interval interval 
1952 : 1952: 

1 20 June 3033 3 July 2946 87 2990 

2 17 July 3039 3 Aug. 2975 64 3012 

3 14 Aug. 3027 31 Aug. 2979 48 3003 

4 18 Sep. 3043 27 Sep. 2990 53 3018 

5 11 Oct. 3027 16 Oct. 2995 32 3012 

6 24 Oct. 3021 1 Nov. 2983 38 3009 

7 18 Novy. 3044 27 Nov. 2961 83 3011 

8 13 Dec. 3058 28 Dec. 2974. 84 3003 

1953 : 1958 : 

2 Feb. 3020 17 Feb. 2966 54 2992 

8S Mar. 3049 16 Mar. 3022 ail 3033 

3 Apr. 3052 15 Apr. 2974 78 3011 

18 June 3054 7 July 2985 69 3017 

16 July 3038 28 July PSS) 39 3014 

7 Sep. 3046 25 Sep. 2993 53 3024 


These results are shown graphically in fig. 10 (b) together with the line of 
best fit given by the method of least squares ; it is assumed that the rela- 
tionship between the D-value and the mean intensity is linear. The mean 
intensity for the whole period under consideration was calculated. The 
D-values were then expressed as percentages of and the mean intensities as 
percentage deviations from the mean of the whole period. The slope of 
the line given by the method of least squares is —0-868. From the large 
negative slope, it is obvious that for smaller peak-to-peak amplitudes of the 
28-day variation, the mean neutron intensity for separate 28-day periods 
is higher and vice versa, which would seem to indicate that the recurrence 
tendency is due to a screening off of part of the primary cosmic radiation. 
This result suggests that the 27-day recurrence tendency is caused by 
decreases in the intensity which occur at 27-day intervals, and which 
supports the view of Neher and Forbush (1952) that the 27 -day variation 
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and the large cosmic ray storms are very similar. It therefore seems 
quite likely that both phenomena can be produced by the same mechanisn. 


Fig. 10 me 
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Value of ) plotted as a function of the mean intensity of each 28-day period. 


§ 4. CONCLUSION 
The present results confirm the existence of a marked 27-day recurrence 
tendency for both the neutron and ionizing intensities, and it is found 
that the amplitude of the 27-day variation is more or less the same for the 
total cosmic ray flux, the hard and soft components, and approximately 
one-half that of the neutron intensity. Some evidence is then given 


27-day Recurrence Tendency of Cosmic Ray Intensity 1251 


that the 27-day recurrence is caused by decreases in the cosmic ray inten- 
sity which occur at 27-day intervals. It can therefore be the same kind 
of phenomenon as the large cosmic radiation decreases which sometimes 
occur at the time of intense geomagnetic storms. The 27-day recurrence 
tendency of the cosmic ray intensity, however, precedes that of magnetic 
activity by about 5 days, and it is tentatively suggested that the 27-day 
intensity changes are associated with solar disturbances, which after a 
delay of a few days, produce geomagnetic field disturbances. 
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SUMMARY 

The purpose of this paper is to bring to notice, and discuss the conse- 
quences of, the phenomenon of ambipolar diffusion of electrons and 
holes down a temperature gradient in a semiconductor. A formula for 
the magnitude of the ambipolar flux is derived, and is applied to the 
electronic contribution to the thermal conductivity and to the Nernst 
Effect. The anomalous thermal conductivity of indium antimonide, 
reported by Busch and Schneider, cannot be accounted for in terms of 
the energy transported by this ambipolar flow. 


$1. LyTrRopUCTION 
THE purpose of this paper is to bring to notice, and discuss the conse- 
quences of, a transport phenomenon in semiconductors which has no 
counterpart in metals. If a temperature gradient is maintained in a 
semiconductor then, since the equilibrium concentrations of electrons 
and holes increase with temperature, there will be a concentration gradient 
of each in the direction of the temperature gradient. Hence diffusion 
will cause a flux of electrons, and a flux of holes, down the temperature 
gradient. (This effect is in addition to that due to the Soret Effect? and 
is in general much greater, as we shall see.) In the absence of an electric 
field these fluxes cannot be expected to be equal, so there will be a net 
electric current. If the specimen is left ‘ open circuit ’, an electric field 
such as to reduce the current to zero will appear. The two fluxes will 
then be equal, one having been increased and the other decreased by 
the field, so that the final effect of the temperature gradient is an 
ambipolar flow of current carriers down the gradient. This flow is of 
particular interest because it contributes extra terms, which can be 
large compared with the ‘ conventional’ ones arising from the electron 
and hole bands considered separately, to the magnitudes of the electronic 
thermal conductivity and the Nernst Effect. The latter phenomena are 
discussed in the following sections. In the present section a formula for 


* Communicated by the Author. 
} The Soret Effect is the drift due to a temperature gradient, given by the 
term in D;? of eqn. (2). 
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the magnitude of the ambipolar flux, for conditions as general as is 
convenient, is derived. The conditions assumed are as follows : 


(a) The specimen is homogeneous and a single crystal. 
(6) The crystal has cubic symmetry. 


(c) The temperature gradient, and its divergence, are not so large as 
to cause appreciable deviations of the carrier densities and distribution 
functions from their equilibrium values: the transport equations remain 
linear. 

(d) Only a single conduction-band minimum energy (that of the lowest 
conduction band), and a single valence-band maximum energy (that of 
the highest valence band) are involved in the transport phenomena. 
(That is to say the energy distances to the next maximum and to the 
next minimum are ezther large compared with or small compared with k7’.) 
We do not require the conduction (valence) band involved to have a 
unique energy minimum (maximum) in pseudomomentum space; so 
that for example germanium and silicon, which have a multiple energy 
minimum in their lowest conduction bands and a degenerate maximum 
(two bands touching) in their highest valence bands, are covered by the 
theory. 


(e) ‘ Impurity band ’ effects (Hung and Gliessmann 1950) are negligible. 


(f) At the band maxima and minima involved, the energy may be 
taken to depend quadratically on the pseudomomentum p, 


e(P)—e(Po) = 4 | p—Po /?, eee 4) 
over an energy range at least of order k7’. In (1) A is independent of 
|p—pyo|; but we do not require it to be independent of the direction 
of p—p,, nor that (1) be an analytic function of p (the ‘ warped ’ valence 
bands of Ge and Si (Dresselhaus et al. 1954) are covered). 


(g) The concentrations of electrons and holes are low enough for 
Boltzmann statistics to be applied. 

Let the number-densities of carriers be 7; and the mobilities, diffusion 
coefficients and Soret-effect coefficients be .;, D; and D,” respectively, 
where the subscript labels the band: i=1 for electrons and i=2 for 
holes. With a temperature gradient grad 7’ and electric field E, the 
drift velocities u, are given by 


—u,=[D; d(log n,)/dT7+D,"] grad T +E. ee ea) 
If we substitute the Einstein relation 
ID ye AU Pati, scab mate (3) 
and set 
D2 = 7, Dil, (4) 


in (2), it becomes 
—u,=p,{(k/e)[T d(log n,)/d7'+-y,] grad T+E}. . (5) 
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For an ‘ open circuit ’ crystal we have equal fluxes : 
,U,;=NoU, =N, say.. . » + . - « (6) 


The value of E required, in (5), to satisfy (6) gives the Thomson 
then results in 


coefficient. Eliminating E in (5) 
Ne Zabol (; “)| 7 d( ae es og mt syste | grad (7) 
Mf +N fg 


Now, on assumptions (c), (d), (f) and (gy) above, we have 
n4n,=const. 7% exp (— A/kT), reer ish 


where 4 is the forbidden energy gap between conduction and valence 
bands. Then, writing for the conductivities 


o=0,+0», o ;=EN jf ;; RTE eee) os Ft k 
(7) becomes (Price 1954) 


N= —22(5)[ ap tdtovty, | ema 7. ae 10) 


We may expect the first term of the square bracket of (10) to be often 
the dominant one: at room temperature for Ge and Si the first term 
is about 30 (somewhat smaller for many intermetallics), while for the 
conventional lattice scattering and spherical energy surfaces y, and yg, 
are each 4. The factor o,0,/o of (10) cannot exceed o/4. For ‘ intrinsic ’ 
conditions it is 


of [Ee y/o 2+ Mo/ey)- 


The essential origin of this flow process is the fact that conduction 
electrons and holes are not separately conserved in number: the hotter 
end of the crystal is a source of pairs and the colder end is a sink of 
them. An analogous phenomenon is the superfluid ‘ internal convection ’ 
in liquid helium IT, according to the Landau (1941) model, in which there 
is a flow of excitations (carrying mass, energy and entropy) relative to 
the centre-of-gravity reference frame. Another analogous phenomenon 
is the counter currents of a dissociating molecule and its constituents 
set up by a temperature gradient in a gas (Dirac 1924), 


§ 2. ELecrronic THERMAL CONDUCTION 


It is obvious that the diffusing pairs must transport oe energy 
required to create them. The energy transported is 4+-(8,+8,)k7. 
where the factors 5; have magnitudes, of order unity, which w ng aeteda 
on the relaxation processes and band structures (for lattice scattering 
and spherical energy surfaces they each have the value 2). Hence the 
total energy transport is 


| Q=N[4+(6,+6.)kT]+@Q,+@,, Pett eh 
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where the Q@; are the standard electronic thermal conduction effects 
contributed by each band considered separately. Combining (10) and (11), 
we arrive at an electronic thermal conductivity of form (Price 1954) 


E lee = OT (k/e)?( 2+ Qo) | 


§2(1.)=(¢105/0?)[( A/kT )?+-«( A/kT') +B] é (12) 


where again « and 8 depend on the detailed relaxation processes, and 
£2, is the ‘ conventional’ contribution from the Q, of (11). The factor 
Q-+ 2, replaces that in the conventional Wiedemann—Franz Law (Seitz 
1940 a) for metals: clearly in semiconductors it has in general a much 
larger value. For intrinsic germanium at 300°xK, for example, we have 


Q=0-22[292-+ 2904 B]: 


this value probably lies between 200 and 300. The anomalously large 
Wiedemann-—Franz ratio for semiconductors does not lead to electronic 
contributions to the thermal conductivity greater than in metals, because 
the electrical conductivity is so much smaller for semiconductors. Thus, 
for Ge and Si at their respective melting points one estimates a contri- 
bution to K,,, from the term of (12) proportional to A? of about 
0-3 watt units in each case. Since K,,, should increase, and the lattice 
contribution to thermal conductivity decrease, with increasing tem- 
perature, this result indicates the possibility of a minimum in the total 
thermal conductivity, K, fairly near the melting point in each case. 

Recently Busch and Schneider (1954) observed a remarkable anomaly, 
of this kind, in the variation with temperature of the thermal conductivity 
of indium antimonide. The value of K has a minimum of a third of a 
watt unit somewhat below 100°c, and has doubled again by a little above 
400°c. It can be interpreted as the sum of two terms, 


KAR (TART), ee (13) 


where K, only is appreciable below 0°c. K, varies as 1/7’, and may be 
identified with the lattice contribution. The other term, A, varies 
according to 
i K,=4-7 exp (— 0-26 ev/2k7’) watt units. tens) 


Busch and Schneider interpret AK, as electronic in origin, as is very 
reasonable considering its temperature variation (14). Because of this 
temperature variation it may be assumed, presumably, that the crystals 
yielding the published data were intrinsic in the temperature range in 
question (100°c-500°c). Since electrical conductivity data for the same 
erystals are not available, we have estimated K,/o7’ by using the 
conductivities published by Tanenbaum and Maita (1953),* in particular 
for their ‘ crystal A’. From the data in the range 60° to 364°C, we find 


K,/oT = 50(k/e)?, experimental. Pe cn (LO) 


« * The following numerical analysis is due to Mr. P. B. Linhart. 
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(The authors’ ‘ crystal B’, evidently more nearly intrinsic, would pre- 
sumably have given a somewhat larger value for (15) if its conductivity 
in the same range had been available.) The factor 50 im (15) compares 
with the value Q,—2 (see eqn. (2) of Busch and Schneider) which is 
obtained by treating each band separately (i.e. not including our diffusion 
effect) and also assuming lattice-scattering relaxation and spherical 
energy surfaces. 

To compare (15) with (12) we estimate 2 as follows: for 4=0-26 ev 
and 7'~500°K, A/k7'=6, and hence the first term in the square bracket 
in (12) is about 36. The other two terms in this case cannot be expected 
to be negligible in comparison ; and in fact for spherical energy surfaces 
and lattice scattering they are such that the whole bracket comes to 100. 
However, if we adopt the mobility ratio p4/u,%85 found by Welker 
(1954), for the intrinsic state the first factor o,0,/c? has a value of about 
1/87, so that Q~1. (The maximum possible value of this first factor is 4, 
giving Q~25: we cannot reach the experimental value (15), from (12), 
even when the energy transport is not reduced by this factor arising from 
the thermoelectric field equalizing the two fluxes.) Thus it appears that 
Busch and Schneider’s effect cannot be accounted for by the energy 
transported by the ambipolar thermodiffusion. 

It seems desirable that the effect be restudied by measuring both 
thermal and electrical conductivities, as a function of temperature, on 
the same single crystal specimen, and further observing the effect of strong 
magnetic fields on both. 


§ 3. Nernst Errect 

We now consider the effect of a magnetic field H, small enough not 
to change the mobilities, crossed with the temperature gradient. Since 
the electrons and holes are drifting in the same direction they will be 
deflected in opposite directions and thus give rise to reinforcing transverse 
currents. The magnitudes of these currents are given each by the product 
of (e/c)N x H with an appropriate mobility. We now show, subject to 
assumptions (b), (c) and (g) of § 1, that this is just the Hall mobility for 
the band. We need only consider a single band in isolation, so drop the 
subscript i. The Boltzmann equation for the distribution function f(r, p) 
is, in our case, 


a oer or enor 
Of are Rad if bel IETS 


where Yf is the ‘ collision derivative’, the rate of relaxation of f, and 
F is the Lorentz force 


F—g[E+(vxH)/c] ie) 


for charge q=-+e and group velocity v(p)=d«/dp. A Lorentz—Gans 
solution (Seitz 1940 b) of (16) is obtained by writing f as a series in orders 
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of the transport ‘ forces ’, f=f,)+/,+..., with the choice of fy such that, 
Bie the operator Y depending appropriately on the local temperature 

i 2i=0 : 
fo=mr)A exp(—<¢kT), . . . ~. . (18) 


where A(7’) is a normalizing factor (proportional to 7/2 on assumption ( f) 
of § 1). We have 


Of Cej,, a d(log n) d(log A) ce Me 
dp ape” =| +(e + Em) 7 | (18°) 
The term on the right of (16) involving f, and H vanishes, so that 
—dropping the terms explicitly proportional to dT /dr,* which we are 
not interested in here—we find f, is a solution of 


—Pf—(gle)(vx A). fy/p=fyv [Vlog n—E(gk7)]. (19) 


Thus an electric field E and a concentration gradient —nE(q/kT’) lead, 
even in the presence of a magnetic field, to the same disturbance 
(expressed by f,(p)) of the velocity distribution. For the component 
of drift velocity parallel to V/ log n—E(q/kT), this result merely verifies 
the Einstein relation (3): for the perpendicular component it shows that 
the ‘ Hall’ mobilities for drift currents due to an electric field and to 
concentration diffusion (our ambipolar effect is a combination of the two) 
are the same. 
If B is the Nernst coefficient, defined by 


Ee = Py xe (20) 


as measured by the transverse field E, for virtually infinite surface 
recombination rate (Landauer and Swanson 1953), then from the 
foregoing results we find 


B= (p47 +19")(o409/0?)(klec)[A/kT +3+Y1+7el+Bo ete) 


where p,;* are the Hall mobilities and by is the ‘ genuine’ Nernst 
coefficient contributed by each band taken separately. To indicate 
the magnitude of the effect we estimate B from (21) for intrinsic 
germanium at 300°K, by dropping B, and setting y;+y,~1. (Strictly 
speaking, supposing y,;=} is consistent only with also supposing 
uf=37p,/8, which is not true; but the y; make a small contribution 
to the square bracket of (21), so we give them this somewhat arbitrary 
value while using experimental values for the p,”.) Substituting 
A=0-75 ev and the values for w; and »,” found by Morin (1954), we 


obtain 


B=4-3 x 10-8 volts/gauss degree, theoretical estimate. (22) 


* These are responsible for the Soret term (the term in D®) of eqn. (2). 
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This is a considerably larger value than, to my knowledge, has been 
found by measurement on any substance. For comparison, a table of 
selected published experimental values follows : 


Substance Temperature, and B in volts gauss degree Source 
Bi =—119°o, 9x10 = 6276, 4-6 4 10-@ 
Te 33°c, 3-6 10-°; 204°C, 2:7x 10-° 
Pb 57°o, —bx10-* International 
Ir 18°o, —5x 10-44 Critical Tables 
Sn 57°a, ? —4x10-4 1929, 6, 420 
Si 47-4°o, 3-3 10- (Gottstein 1914) 
50°o, 1-1 x 10-!° and 2-0 x 10~!° (Buckley 1914) 
PbTe 50°c, 1-4 10-® (maximum value) 
PbSe 420°c, 4 10-1 (maximum value) lee 


For intrinsic silicon, (21) would predict a value for B of order 
10-* volts/gauss degree. However, appreciable doping could reduce the 
value considerably: if we drop the B, term from (21), and consider 
various dopings, at fixed temperature, not so large as to change the 
mobilities, then the constancy of n,n, leads to the rule* 


Bo?=constant=(Bo*), ntrinsic- Se USN ee eae 


In Buckley’s paper he mentions that his Si sample had several per cent 
of impurities (mostly Fe, some Al, and smaller fractions of other elements). 
This amount is, of course, far beyond the range of validity of (23), because 
of the effect of impurities on mobility and on band structure and crystal 
perfection, apart from the fact that B, will no longer be the negligible 
term in (21). It would evidently be of interest to make Nernst Effect 
measurements on Ge and Si erystals of present-day purities. 
Putley (1955) quotes a formula which reduces to 


377 a,0,fk\[ A 37 [ k\nypy?+nops” 
B= — vg) ge (| ee ed ed ee 
8 (Ha tHe) a (=) a+ | vs (se) Nip +Nofo ee 
in the present notation and units. It is evident by comparison with (21) 
that 


(a) the first term of (24) corresponds to the first term of (21) and the 
second term of (24) to the second term (B,) of (21); and 

(b) (24) has been derived by assuming spherical energy surfaces and 
lattice scattering. 
Thus (21) is a generalization of (24) ; but (24) gives an explicit expression 
for By for a special case. In the conditions of Putley’s observations on the 
intermetallic semiconductors PbTe and PbSe, the two terms of (24) were 


— 


* This result might be used as a test to determine when By is negligible 
compared with B. 
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of comparable magnitude. Because they are of opposite sign, the Nernst 
coefficient as a function of temperature passed through a zero in each case. 
He Since in practice semiconductor crystals even of highest purity become 
non-intrinsic at a low enough temperature, the first (thermodiffusion) 
term of (21) will have a maximum as a function of temperature, falling 
off to low values at both ends of the scale.* For a reasonably pure 
erystal of a semiconductor with a large energy gap and not too unequal 
mobilities (such as Ge and Si), however, there should be a range of 
temperature where we may both neglect B, and approximate this first 
term by 

Be (uy4+12")(o109/0?)(d/cT’) pe eee (21) 
where 4=e¢. We consider now what information could be obtained by 
determining B for a fixed temperature and different dopings, in conditions 
where (21’) is a good approximation. It is convenient to define a ‘ Nernst 
mobility ’ 


bY = BET /$)=(us7 +49")039,/0°, ea) 
just as the Hall mobility is defined by 
ee = oR=(u,"0,—p5",)/c. ey eae 20) 


Using the fact that o,+o0,=0c, we may eliminate the o’s from (25) and 
(26). Then 
p= (oy — pF) (eo +) (uy +9"). 5 Spee OE8) 

Over the doping range for which the band mobilities are constant, the 
plot of u* against ~” should thus lie on a parabola. The vertex of the 
parabola (maximum of «%) should be for the same doping as that for 
which the conductivity is a minimum, and have the coordinates 

| PE —po™), w= d(uy tye"). . -  ~ (28) 
The determination of the point (28) thus would give the values of wu,” 
and yp.” separately, 6 being known. This is the same information as is 
obtainable from the Dunlap ellipse (Dunlap 1950)—though in the present 
method it is obtained from a single portion of the doping curve, that 
around the resistivity maximum. However, if there were difficulty in 
getting absolute values of B, or uncertainty over the value of 4, the 
present method would determine only ».,%—p,": this is just what can 
be determined from the neighbourhood of the resistivity maximum 
by the Dunlap ellipse. Furthermore (21’) is only an approximation, 
with a limited area of usefulness. 


* There will in consequence be a change in sign of B (to negative at the 
lowest temperatures) where extrinsic conditions set in. There can also be a 
second change of sign (to negative at the highest temperatures), if (21) becomes 
negative (By, dominant) when we set n;=n, and A/kT=0. (In Putley’s 
formula (24)—that is, for spherical energy surfaces and lattice scattering— 
the condition for this is that either ~,>7:87u, or w2>7-87u,: this condition 
is satisfied for InSb.) It might be, however, that sufficiently heavy doping 
could dépress the first term of (21) far enough for there to be no changes of 
sign, the value of B remaining negative at all temperatures. 
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There is an analogy between the phenomenon, contributing to the 
Nernst Effect, discussed above and the phenomenon underlying the 
photomagnetoelectric effect (Moss, Pincherle and Woodward 1953), in 
which electrons and holes generated by light falling on a semiconductor 
surface diffuse inward together and are deflected by a magnetic field 
parallel to the surface, so that they maintain an electric field perpendicular 
to the magnetic field and parallel to the surface. 
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CXXXVI. CORRESPONDENCE 


Observations on Etched Crystals of Barium Titanate 


By D. 8. CampBein 
The Plessey Co. Ltd., Caswell Research Laboratories, Towcester, Northants.* 


[Received August 19, 1955] 


OBSERVATIONS of the etching effect of cold concentrated hydrochloric 
acid on barium titanate single crystals have been reported by Hooton 
and Merz (1955). We have obtained similar results using hydrochloric 
acid and phosphoric acid on crystals direct from the melt, with no 
polarizing treatment. Figures 1 and 2 (figs. 1-4, Plate 28) show the effect 
of etching a crystal containing two wedges of optically active material 
(‘a axis normal to the surface). 

Figure 1 shows the crystal examined with crossed nicols and mono- 
chromatic light under which conditions the interference fringes produced 
by the wedges are clearly visible. Figure 2 shows the surface of the 
crystal after etching in hydrochloric acid for 1 minute. The positive 
ends of the ‘c’ domain areas etch most heavily, the ‘a’ areas less 
heavily and the negative ends of the ‘c’ areas least of all. It is to be 
noted that fig. 2 shows the ‘c’ domains running parallel to the wedges. 
That the ‘c’ domains run vertically through the crystal is shown in 
figs. 3 and 4, where the correspondence between the top and bottom 
surfaces is very good. 

Fig. 5 


Diagram of coupling through a wedge. 


Examination of a wedge shows coupling through the wedge as illustrated 
in fig. 5. Dark areas, unetched, at one side of the wedge, correspond to 
light areas, etched, on the other. Figures 6 and 7 (figs. 6-9, Plate 29) 
show the type of micrograph obtained. 


* Communicated by H. F. Kay. 


1262 Correspondence 


Observation of several wedges has shown that coupling does not hold 
after a certain distance (approximately 0-07 mm) indicating some break- 
down in the head to tail domain arrangement over relatively large 
distances. 

Thicker crystals have shown similar results. Figure 8 shows the top 
and fig. 9 the bottom surface of a crystal 0-16 mm thick and it can be 
seen that coupling does exist but is very poor, the vast majority of the 
domains present at the top surface not appearing at the bottom. 

The etching technique has proved useful for the examination of barium 
titanate crystals. Failure of head to tail coupling after a certain minimum 
distance is observed, so that thin crystals should respond more easily than 
thick ones to applied fields. An electron microscopy study of the etched 
surface is being undertaken and will be reported elsewhere. 


The author would like to thank the Plessey Co. for permission to 
publish this note. 
REFERENCE 
Hooton, J. A., and Merz, W. J., 1955, Phys. Rev., 98, 409. 
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Fig. 1 


Hh 


Interference pattern produced Effect of etching in cone. HCl 
by ‘a’ wedges. x 100. for lmin. 100. 
Fig. 3 Fig. 4 


Top surface of ‘ ¢’ region crystal Bottom surface of ‘ ¢’ region crystal 
0-048 mm thick. Ultropak. 530. 0-048 mm thick. Ultropak. 350. 


D. S. CAMPBELL Phil. Mag. Ser. 7, Vol. 46, Pl. 29. 


Top surface of wedge area. Bottom surface of wedge area showing 
Ultropak. x 350. reflection. Ultropak. x 350. 


Top surface of “¢’ region crystal Bottom surface of ‘ ¢’ region crystal 
0-16 mm thick. Ultropak. 200. 0-16 mm thick. Ultropak. x 200. 


CXXXVIT. REVIEWS OF BOOKS| 


Solubilization and Related Phenomena. By M. KE. L. McBartn and E. Hurcutiy- 
son. [Pp. xv+259.] (Academic Books Ltd.) Price 56s. 

SOLUBILIZATION is an awkward word coined to express the increased take-up 
of a wide variety of substances in certain colloidal solutions. The widespread 
technological importance of the phenomenon is familiar from the enormous 
growth of the detergent industry ; less generally familiar technical applications 
of solubilization are found in practices such as emulsion polymerization, the 
solubilization of dyestuffs for treating textiles and food, and the solubilization 
of proteins in certain processes for the manufacture of fibres. 

This monograph aims to reduce a wealth of empirical knowledge and 
application to some kind of fundamental principles. After a brief survey of 
the historical background and of the properties of colloidal electrolytes some 
95 pages are devoted to a very useful factual survey of experiments on 
solubilization. 

Some kind of aggregate or micelle formation plays an essential part in the 
process. But discussions of the mechanism of solubilization as presented in 
this book make it clear that experimental evidence on the precise structure 
of such micelles still is very scanty. The present review of the experimental 
findings and of the current controversies of interpretation will be found 
valuable both by specialists in this field, and by those who are interested in 
the progress of physico-chemical theory generally. A. R. UBBELOHDE. 


Optical Properties of Thin Solid Films. By O. 8. Heavens. [Pp. 261.] 
(London: Butterworths Scientific Publications, 1955.) Price 35s. 
Tuts book, by Dr. O. 8. Heavens, dealing with the optical properties of thin 
solid films is a valuable and timely contribution to a subject which, within 
recent years, has grown to be of extreme importance both scientifically and 
technologically. Written by an expert the book covers (1) the formation of 
thin films, (2) the structure of thin films, (3) thin film optics, (4) the measurement 
of thin-film thickness and optical constants, and (5) practical applications of 
thin films in optics. The experimental production of solid thin films is 
admirably described by an author who himself has clearly had much practical 
experience in the production of such films. Of particular value is the discussion 
of the electron-microscopic studies of the structures of thin films. Thin-film 
optics, which is a formidable subject, is treated sympathetically and with 
the minimum of mathematical analysis needed for the study of such a complex 
subject. This section is unavoidably the heaviest in the book, yet at the same 
time will probably prove to be the most valuable in that it will save the reader 
rom the difficult task of reading many papers, mostly in foreign languages, and 
at times obscure. It is particularly in this chapter, that dealing with film 
optics, that the author renders a considerable service to those who are working 
in this field, especially to the practically minded man. The various experi- 
mental methods available for measuring thin-film thicknesses are discussed. 
in. detail. These include mechanical methods, electrical methods, optical 
methods and interferometric methods. The use of radioactive tracers is 
briefly indicated. The last chapter, dealing with practical applications, covers 
a very wide field such as that of anti-reflection systems, interference films, 
interference filters, multilayers, etc. and although consisting of a bare 50 pages 
is in fact a summary of a very considerable number of substantial publications. 
The book is practically free from blemish. The reviewer has only noted in 
two minor instances incorrect allocation of priority for some idea or experiment. 
This is of little consequence. 
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Production is excellent and the cost not excessive as far as modern prices go. 

There is no doubt that this will become a most valuable source of reference 
to all workers in the field of thin-film optics. The book is hardly to be 
recommended for the undergraduate, and was not intended for readers at this 
level. It was intended for the research worker, and as such, we welcome this 
as a valuable contribution to a growing and difficult subject. Sg 


Quantum Theory of Solids. By R. E. Prrerts. [Pp. vi+229.] (Oxford : 
Clarendon Press, 1955.) Price 30s. 

Tue chief danger to a student of solid state physics is that he will acquire a 
glib facility for manipulating various phenomenological concepts (e.g. phonon, 
hole, collision, dislocation, relaxation time) without really understanding their 
limitations. Professor Peierls begins with the atomic lattice and deduces its 
behaviour—thermal and electrical conduction, cohesion, magnetism, inter- 
action with light, superconductivity—directly from first principles. The book 
exhibits a rare rigour, not of pedantic mathematical formalism but in its faith 
to the physical model. The omission of numbers, graphs, detailed experimental 
results and complicated analytical formulae leaves the argument unobscured, 
whilst real fundamental difficulties are emphasized. The lucid exposition of 
early work not previously available in English is especially valuable, although 
there are places (e.g. in discussing electron—electron interaction) where the 
theory is not quite up to date. Students will find it simple and clear, making 
an admirable introduction to the basic theory ; all who work in the field will 
benefit from the criticism of many fallacious over-simplifications which are 
widely current. J. M. Z. 


High-Energy Accelerators. By M. 8S. Livinaston. (New York and London: 
Interscience Publ. Inc., 1954.) [Pp. viii+157.] Price 26s. 


THERE could hardly be an author more suitable for this kind of review than 
Professor Livingston, who has been constructing high energy accelerators from 
their very beginning, and as a general common-sense introduction to the subject 
the book could not be better. Most of it is concerned with the actual physics of 
the various accelerating systems ; only the first tenth of the book discusses their 
applications and the reasons for the drive towards ever higher energies. After 
the second chapter, in which the general physical principles are explained 
(relativistic motion, oscillations about stable orbits etc.), the next four chapters 
are concerned with the four main types of high energy machines: electron 
synchrotron, synchrocyclotron, linear accelerator and proton synchrotron. The 
last chapter describes the alternating-gradient focusing (‘ strong focusing ’) 
which has made possible the latest push towards energies above 101° electron 
volts. Mathematical formulae are given where needed but their derivation 
is only briefly indicated. Points of engineering are mentioned but not described 
in detail, and there is only brief mention of the problems of shielding. 

The construction of high energy machines has become an extremely complex 
and highly specialized subject which nobody could expect to learn from a book, 
but as a quick survey or as an introduction to the more specialized reports, this 
little book will serve admirably. O. Rak 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. ] 


